*H FLOATING POINT 

BINARY FLOATING POINT REPRESENTATION 

24 BIT SIGN-MAGNITUbE NORMALIZED FRACTION. 

8 BIT EXCESS-128 EXPONENT (BASE 2) 

SIGN BIT REPLACES MOST SIGNIFICANT FRACTION BIT (IMPLIED 11 
EXPONENT FIELD OF ZERO IMPLIES VALUE IS ZERO 


7F FF FF FF 
00 00 00 81 
00 00 00 01 
XX XX XX 00 
FF FF FF FF 


= 16777215/16777216*2*127 ■ 1.7014117*10*38 
= 1 / 2 * 2*1 =: 1.0000000 
= 1/2*2*(-127) = 2.9387359*10+(-39) 

= 0 

- “16777215/16777216*2+127 - -1.7014117*10+38 


ALL RESULTS ARE ROUNDED TO 24 BIT FRACTIONS. 


DECIMAL FLOATING POINT REPRESENTATION? 


13 BCD DIGIT SIGN-MAGNITUDE NORMALIZED FRACTION 
8 BIT EXCESS-128 EXPONENT (BASE 10) 

SIGN BIT IS MOST SIGNIFICANT DIGIT 
EXPONENT FIELD OF ZERO IMPLIES VALUE IS ZERO 


09 99 99 99 99 99 99 FF = 

01 00 00 00 00 00 00 81 = 1 

01 00 00 00 00 00 00 01 = 

XX XX XX XX XX XX XX 00 = 0 

89 99 99 99 99 99 99 FF 


99999 99999 999 * 10+127 


1 * 10+(-127) 


-.99999 99999 999 * 10+127 
ALL RESULTS ARE TRUNCATED TO 13 DIGIT FRACTIONS. 
CALLING SEQUENCE: 


IX POINTS TO FIRST BYTE OF FIRST ARGUMENT 
IY POINTS TO SECOND ARGUMENT (IF ANY) 

RESULT REPLACES FIRST ARGUMENT, IX IS UNCHANGED 

VALUE OF SECOND ARGUMENT IS UNCHANGED 

NO OTHER REGISTERS PRESERVED (INCL. ALTERNATE BANK) 

STANDARD FUNCTIONS: 


SUB 

X'=X-Y 



ADD 

X - X+Y 



MULT 

X 

1! 

X 

* 

-< 



DIVI 

X - X / Y 



RND 

X=RND PSUEDO-RANDOM NUMBER 

0< = X 

RANDMZ 

RANDOMIZES RND SEQUENCE 

I NT 

H=INT(X) 

FLOOR FUNCTION 


ABS 

X - A B S ( X ) 

ABSOLUTE VALUE 


SGN 

X=IF X = 0 

THEN 0 ELSE IF X>0 

THEN 

NEG 

X--X 



ROW 

X=X + Y 



SQR 

X-SQR(X ) 



EXP 

X=EXP(X) 



LOG 

X s3 LOG ( X ) 

NATURAL LOG 


TAN 

X=T AN(X) 

X IN RADIANS 


COS 

X-COS(X) 

X IN RADIANS 


SIN 

X=SIN(X) 

X IN RADIANS 


ATN 

X-ATN(X) 

-PI/2 <- RESULT < = 

PI/ 





RELATIONAL FUNCTIONS (RETURN 1 FOR TRUE, 0 FOR FALSE) 


EQ 

> 

If 

X 

11 

X 

GE 

X = ( X > = Y ) 

GT 

X=(X>Y) 

LE 

X-(X<=Y) 

LT 

X = ( X < Y ) 

NE 

X=(XOY) 


BOOLEAN FUNCTIONS (READ NONZERO AS TRUE, ZERO AS FALSE); 


AND X = X.AND.Y 

OR X=X.OR.Y 

NOT X=„NOT.X 

SPECIAL FUNCTIONS; 

RSIJB X = Y ) - X 
REC X=l/X 

BRK X=X-INT(X) 

ONE X=1 

ZERO X*0 

INF X=IF X>=0 THEN (LARGEST POS. #) ELSE (LARGEST NEG 

*H INTEGERS 

INTEGER REPRESENTATION: 

16 BIT 2'S COMPLEMENT 

REVERSED BYTES (LOW ORDER BYTE AT LOWER ADDRESS) 

FF 7F = 32767 
01 00 = 1 
00 00 = 0 
FF FF = -1 
00 80 = -32768 

CALLING SEQUENCE; 

IX POINTS TO FIRST BYTE OF FIRST ARGUMENT 
IY POINTS TO SECOND ARGUMENT (IF ANY) 

RESULT REPLACES FIRST ARGUMENT 

IX IS UNCHANGED 

NO OTHER REGISTERS PRESERVED 


STANDARD FUNCTIONS: 


IADD X = X + Y 

ISUB X*X-Y 

IMUL X = X * Y 

IDIV X=X/Y 

IABS X=ABS(X) 

ISGN X=IF X = 0 THEN 0 ELSE IF X<0 THEN -1 ELSE 1 
RELATIONAL AND BOOLEAN FUNCTIONS (SEE FLOATING POINT DEFINITIONS) 


IEQ,IGE,IGT,ILE,ILT,INE 
IAND,IOR,INOT 

SPECIAL FUNCTIONS (SEE FLOATING POINT DEFINITIONS): 


IONE,IZERO,INEG 


*H CONVERSION 
CONVERSION FUNCTIONS : 

FIX CONVERTS X FROM FLOATING POINT TO INTEGER FORMAT 

FLOAT CONVERTS X FROM INTEGER TO FLOATING POINT FORMAT 

(NOTE THAT 4 OR 8 BYTE RESULT REPLACES 2 BYTE ARG) 
DTOB CONVERTS ASCII CHARS TO FLOATING POINT NUMBER 

(NOTE SEE LISTING FOR CALLING SEQUENCE) 

BTOD CONVERTS FLOATING POINT NUMBER TO ASCII STRING 

(NOTE SEE LISTING FOR CALLING SEQUENCE) 

MODULES; 

BMATHO,BMATH1,BMATH2,BMATHM 
DMATHOrDMATHl,DMATH2,DMATHM 


FILE RECORD RECORD FILE STARTING DATE OF DATE OF 


FILENAME 

DRIVE 

TYPE 

COUNT 

LENGTH PROPS ADDRESS 

CREATION 

LAST MOD . 

DMATHO.S 

2 

A 

182 

0080 

770801 

780718 

DMATH2.S 

2 

A 

147 

0080 

770801 

780830 

DMATHI.S 

2 

A 

127 

0080 

770919 

780718 

DMATHO.OBJ 

2 

B 

28 

0080 

770801 

780718 

DMATH1.OBJ 

2 

B 

18 

0080 

770801 

780718 

DMATH2.OBJ 

9 

B 

19 

0080 

770801 

780830 

DMATHM.S 

2 

A 

8 

0080 

770914 

780718 

DMATHM.OBJ 

9 

B 

3 

0080 

770914 

780718 

BMATHO.OBJ 


B 

22 

0080 

770802 

780718 

BMATH.1. .OBJ 

2 

B 

19 

0080 

770802 

780718 

BMATH2.OBJ 

9 

B 

22 

0080 

770802 

780718 

MATHLIB.S 

9 

A 

35 

0080 

770802 

780718 

BNSRCH.S 

9 

A 

8 

0080 

770912 

780801 

BMATH1.S 

2 

A 

156 

0080 

770802 

780718 

BMATHO.S 

9 

A 

147 

0080 

770802 

780718 

BMATH2.S 

2 

A 

134 

0080 

770802 

780718 

MATH.DOC 

2 

A 

26 

0080 

* * « s r. t 

780718 

DIV3.S 

2 

A 

14 

0080 

770815 

780828 

BMATHM.S 

9 

A 

8 

0080 

770914 

780719 

BMATHM.OBJ 

9 

B 

3 

0080 

770920 

780719 

CAL . S 

2 

A 

120 

0080 

780718 

780828 

ADM.S 

9 

A 

6 

0080 

780719 

780719 

BEEHIVE, S 

9 

A 

6 

0080 

780719 

780719 

CAL.CMD 

9 

A 

1 

0080 

770920 

780719 

INFOTON.S 

26 FILES 

2 

EXAMINED 

A 

6 

0080 

780719 

780719 


25 FILES LISTED 

1290 TOTAL SECTORS FOR LISTED FILES 
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*H ADD AND SUBTRACT 


FLOATING POINT REVERSE SUBTRACT ROUTINE 
NEGATE X AND ADO 

rryryyrryyyyyyyryyyyyyyrrryyyyyyyyyyyyyyyryryyyyyrryryyyyyyryyyryyyyyyy 


L.D 

A , (IX ) 

XOR 

SDH 

LD 

(IX) ,A 

LD 

H, A 

LD 

D,(IY) 

JP 

ADOS 


I SET NEW SIGN 


FLOATING POINT SUBTRACT ROUTINE 
NEGATE Y AND ADD 


SUB LD A , ( IY ) ? GET SIGN OF Y 

XOR 80H ; CHANGE IT 

LD D , A 
L D H, ( IX ) 

JP ADOS 


FLOATING POINT ADD ROUTINE 

INPUT: IX - FIRST ARC & RESULT 

IY - SECOND ARG 
OUTPUT: IX - SUM 

ERRORS: OVERFLOW, UNDERFLOW 

REGISTERS FOR ADDITION: 

H'L'HL. ARG WITH LARGER EXP 
D’E'DE 2ND ARG 

B ’ BIT 7 IS SIGN OF RESULT 

C' LARGER EXPONENT 


t it* 


7777777777777777777777777777777777777777777777777 


77777777777777777 


1 

« 

« 

< 

< 

« 

« 

« 

I 

i 

1 

« 

I 

i 


; GET ME 

T 

: BYTES 



r 

ADD 

LD 

H , < IX ) 


1 


LD 

Of(IY) 



r 

: CHECK 

FOR 

ZERO ARGUMENTS 

1 

T 

ADDS 

XOR 

A 

J CLEAR A 



CP 

( IX + 3 > 


< 


JP 

NZ, ADOS 

; X IS NONZERO 



CP 

( IY + 3) 




RET 

Z 

; Y IS ALSO ZERO 

1 


LD 

A, D 

; GET Y MS BYTE 



PUSH 

IY 




POP 

HI¬ 


1 


CALL 

LO 

; COPY Y TO X 



LD 

< IX) f A 

? STORE MS BYTE 



RET 



4 

ADOS 

LD 

(IX ) f H 

; STORE NEW X SIGN 



CP 

(IY+3) 

f IN CASE Y IS ZERO 



i 


i 





RET 


; Y IS ZERO, RETURN X 


r 

; LOAD ARGUMENTS, TEST SIGNS 
; H r L'H IS X FRACTION 

J D'E’D IS Y FRACTION 


7 

L 

IS X EXPONENT 

1 

E 

IS Y EXPONENT 

r 

B" 

IS X SIGN 


i 

c ’ 

IS Y SIGN 


r 

ADD? 

LD 

A, H 

} EXTRACT SIGN 


AND 

80H 



LD 

B, A 



SET 

7 , H 

; RESTORE IMPLIED 1 


LD 

A, D 



AND 

80H 



LD 

C y A 



SET 

7, D 



XOR 

B 



EX 

AF,AF ' 

? SAVE XOR OF SIGNS 


LD 

L,(IX + i ) 



LD 

E,(IY+1) 



EXX 


f MAIN BANK 


LD 

H, < IX + 2) 



LO 

Dr«IY+2) 



LD 

L,(IX+3) 



LD 

E,(IY+3) 


1 

; put 

LARGER 

EXP IN HL. 


i SHIFT de right to make 

! EXP'S EQUAL 

r 

A DIO 

LD 

A,L. 



CP 

E 



JR 

NO, AD20 

5 NO SWAP NEEDED 


EX 

DE,HL 

; SWAP X AND Y 


EXX 


5 ALT BANK 


EX 

DE, HL 



LD 

B, C 

; B’ IS SIGN OF LARGER 


EXX 


} MAIN BANK 

AD20 

CALL 

ALIGN 

5 ALIGN FRACTIONS 

; DEC 
* 

IDE WHETHER TO ADD 

OR SUBTRACT 

y 

EX 

AF, AF' 

} RETRIEVE XOR OF SIGN 


JP 

M r SB20 

1 SIGNS ARE DIFFERENT 

r 

I ADD 




f 

ADD 

HL , DE 

? ADD FRACTIONS 


LD 

DE,80H 

; ROUNDING CONSTANT 


EXX 


; ALT BANK 


ADC 

HL , DE 


7 

; NORMALIZE 

(AT MOST i S 

HIFT RIGHT) 

r 

JP 

NC,FROUND 

i ALREADY NORMALIZED 


INC 

c 

? INCR EXPONENT 


JP 

Z,OVERB 

; OVERFLOW 


RR 

H 

,’ SHIFT RESULT RIGHT 


RR 

L. 







MAIN BANK 



EXX 

RR 

H 

; MAIN BANK 


RR 

L 



r w 
tZ i\f\ 


; ALT BANK 


UP 

FROUND 


T 

I SUBTRACT 



T 

SB20 

AND 

A 



SBC 

HL , DE 

i SUBTRACT FRACTIONS 


LD 

DE , BOH 

PREPARE FOR ROUNDING 


EXX 


; ALT BANK 


SBC 

HL 7 DE 



JP 

NC.SB30 

; X WAS > Y 

r 

; IF RE 

SUIT I 

S NEGATIVE, NEGATE HL 

7 

LD 

A , 80H 

} CHANGE SIGN 


XOR 

B 



LD 

B, A 



EXX 


; MAIN BANK 


EX 

DE, HL 



LD 

HL , 0 



AND 

A 



S BC 

HL ,DE 

; FORM NEGATIVE 


LD 

DE,80H 

? PREPARE FOR ROUND 


EXX 


? ALT BANK 


EX 

DE , HL 



LD 

HL , 0 


t 

SBC 

HL, DE 


7 

; post 

NORMAL 

.IZE BY SHIFTING LEFT 

SB30 

JP 

M,FROUND 

; ALREADY NORMALIZED 


EXX 


? MAIN BANK 


LD 

B, 2 4 

J PREPARE FOR POST-SHIFT 

SB 40 

ADD 

HL ,HL 

r SHIFT RESULT LEFT 


EXX 


? ALT BANK 


DEC 

c 

? DECR EXPONENT 


•JP 

Z,UNDER 

; UNDERFLOW 


ADC 

HL , HL 



JP 

M,FROUND 

; NORMALIZED 


EXX 


? MAIN BANK 


D JNZ 

SB 40 

? LOOP NOT DONE 


JP 

ZERO 


r 

; ALIGN 
* 

ARGUMENT IN DE 


7 

1 SHIFT RIG 

;HT UNTIL 

EXPONENTS EQUAL 

? RETURN X 

EXP IN C' 


1 

DE IS 

; ALIGNED 

FRACTION 

T 

L IS 

ZERO 


7 

ALIGN 

LD 

A , L 

? GET X EXP 


LD 

L , 0 

f CLEAR LSBYTE OF X 


EXX 


} ALT BANK 


LD 

C r A 

f PUT IN C' 


EXX 


MAIN BANK 


SUB 

E 

? CALC SHIFT COUNT 


LD 

E, L 

J CLEAR LSBYTE OF Y 



ALIO 


RET Z j SHIFT COUNT IS ZERO 

ID B,25 ; MAX COUNT 


ALIO EXX 


r 

ALT BANK 

SRL 

D 

t 

SHIFT RIGHT 

RR 

E 



EXX 


r 

MAIN BANK 

RR 

D 



RR 

E 



DEC 

A 

r 

DECR SHIFT COUNT 

RET 

Z 

t 

r 

COUNT IS ZERO 

D JNZ 

ALIO 

T 

CONTINUE SHIFTING 

RET 




*H MULTIPLY 




7 J M M M M M 

T*T*%**¥* 

777777777 

J ? y r 

771777777777777777 


r 

; FLOATING POINT MULTIPLY ROUTINE 
1 PRODUCES ROUNDED RESULT 
J INPUT? IX POINTS TO MULTIPLIER 

? IY POINTS TO MULTIPLICAND 

; OUTPUT; IX POINTS TO NORMALIZED PRODUCT 
5 ERRORS; EXPONENT OVERFLOW OR UNDERFLOW 
! MAIN LOOP REGISTER ASSIGNMENTS; 



AHL 

- 24 BIT ACCUM, INITIALLY ZERO 

r 

H ’ L r 

D’ - 24 BIT RIGHT 

EXTENSION 

r 


OF ACCUM, INIT 

MULTIPLIER 

t 

y 

CDE 

- 24 BIT MULTIPLIC 

AND 

r 

B 

- LOOP COUNTER 


r 

A’ 

-■ BIT 7 IS SIGN OF 

RESULT 

T 

TIMING; 

APPROX 2531 (1943 

INNER LOOP) 

y : 

tfTT77TTTrr 

**»***%»»*¥******» 

yyyyyyyyyyyy 

r 

r 

CHECK FOR 

ZERO OPERANDS 


MULT XOR 

A 



CP 

(IX+3) 



JP 

Z,ZERO ; X IS 

ZERO 


CP 

(IY+3) 



JP 

Z,ZERO ; Y IS 

ZERO 

? 

y 

LOAD OPERANDS> INITIALIZE 


} 

LD 

H,(IX) ; X IN 

H'L’D’ 


LD 

L,(IX+1> 



LD 

D,(IX+2) 



LD 

A,H i SAVE 

SIGN 


SET 

7,H ; RESTORE IMPLIED 


EXX 

; MAIN 

BANK 


CALL 

L.OADY ; Y IN 

CDE & FINAL 


EX 

AF,AF' r SAVE 

SIGN 


XOR 

A ; ZERO 

ACCUM AHL 


LD 

H, A 



LD 

L, A 



LD 

B,24 ; LOOP 

COUNTER 


JP 

ML1 


i 

y 

MAIN MULTIPLY LOOP 


T 

PRODUCES 

48 BIT RESULT IN 

AHLH'L'D' 


; SHIFT ACCUM RIGHT 


ML OOP 


RRA 


RR 

EXX 

L. 

RR 

H 

RR 


RR 

EXX 

D 

JR 

NC , ML.2 

ADD 

Hi ,DE 

ADC 

A» C 

DJNZ 

ML OOP 


r 

5 NORMALIZE (AT MOST .1 


i ALT BANK 

; SHIFT MPL.IER RIGHT 


; MAIN BANK 
; MPLIER BIT IS 0 
; ADD MPL.ICAND 
? MAY CAUSE CARRY 
; COUNTER NOT DONE 

SHIFT RIGHT) 


JR NCr MROUND 
RRA 

R R H 

R R L 

JP MROUN1 


ALREADY NORMALIZED 
SHIFT RESULT RIGHT 

ROUNDING BIT IN CY 


; ROUND 


MROUND 

EXX 


J 

ALT BANK 


RL 

H 

T 

GET ROUNDING BIT 


E XX 


r 

MAIN BANK 


DEC 

B 

r 

SET B TO -1 

MROlINi 

CALL 

RSTO 

t 

T 

ROUND AND STORE RESULT 

r 

; ADD 

EXP'S 

AND SHIFT 

COUNT 

; o 
* 

OR 2 OVERFLOWS 

OK, 

1 IS ERROR 

T 

ADD 

A r C 

f 

ADD OTHER EXPONENT 

MCHK 

JP 

PE r M10 

T 

OVERFLOW 


ADD 

A , B 

r 

ADD SHIFT COUNT 


JP 

PE » MERR 

T 

OVERFLOW 


JP 

M20 



M 10 

ADD 

A i B 

T 

ADD SHIFT COUNT 


JP 

PO , MERR 

f 

NO OVERFLOW 

M 20 

ADD 

A t BOH 

T 

ADD BIAS 


JP 

Z,UNDER 

} 

UNDERFLOW 


LD 

(IX + 3 ) > A 

y 

STORE EXP 


RET 





*H DIVIDE 


! FLOATING POINT DIVIDE ROUTINE 
, PRODUCES ROUNDED RESULT 

! INPUT: IX POINTS TO DIVIDEND 

} IY POINTS TO DIVISOR 

! OUTPUT: IX POINTS TO QUOTIENT 

; ERRORS: DIVISOR IS ZERO 

: EXPONENT OVERFLOW OR UNDERFLOW 

t REGISTER ASSIGNMENT FOR MAIN LOOP: 

J AHL - 24 BIT ACCLiM, I NIT DIVIDEND 

i CDE - 24 BIT DIVISOR 

5 D'E'H'L' - RESULT REGISTER 

i B -- LOOP COUNTER 

J A ' BIT 7 IS SIGN OF RESULT 

J TIMING: APPROX 3595 (3022 INNER LOOP) 




; CHECK FOR ZERO OPERANDS 


D IVI 

XOR 

A 



CP 

(IY+3) 



JP 

Z.DIVO ; Y 

IS ZERO 


CP 

(IX+3) 



JP 

Z r ZERO 5 X 

IS ZERO 

r 

; LOAD 

OPERANDS r CALCULATE S 

iIGN 


LD 

Br(IX) 

T 

LD 

H,(IX+1) 


LD 

L., (IX + 2) 


LD 

A,B 

T 

CALL 

LOADY 

r 

EX 

AFrAF' 

r 

LD 

ArB 

i 

OR 

8 QH 

F 

LD 

B, 26 

r 

JP 

DVSIJB 



LOAD X INTO BHL 


GET SIGN IN A 
LOAD Y INTO CDE 
SIGN BIT IN A 7 
X IN AHL 

RESTORE IMPLIED 1 
LOOP COUNTER 


MAIN DIVIDE LOOP 

PRODUCES 26 BIT RESULT IN (BIT 0 OF)D’E'H'L'CY 


DVOV 

AND 

A 


SBC 

HI. i DE 


SBC 

A r C 

DVRESi 

SCF 



DEC 

B 


JR 

Z,DVNORM 

DVLOOP 

EXX 



ADC 

HL,HL 


RL 

E 


RL 

D 


EXX 



ADD 

HLrHL 


ADC 

A r A 


JR 

C t DVOV 

DVSIJB 

SBC 

HL, DE 


SBC 

A r C 


JR 

NC,DVRESI 


ADD 

HL.. DE 


ADC 

A, C 


AND 

A 


DJNZ 

DVLOOP 

r 

; NORMAL 

. IZE 

RESULT 

y 

DVNORM 

EXX 



L.D 

B,0 


LD 

A r E 


BIT 

0 r D 


JR 

ZiDROUND 


INC 

B 


SCF 



RRA 



RR 

H 


RR 

L. 


SUBTRACT DIVISOR 

SET CARRY 
DEC LOOP COUNTER 
LOOP DONE 
ALT BANK 

SHIFT RESULT INTO D’E'H’L 


MAIN BANK 
SHIFT ACCUM LEFT 


OVERFLOW ON SHIFT 
SUBTRACT DIVISOR 

POSITIVE RESULT 
RESTORE ACCUM 

TURN OFF CARRY 
LOOP NOT DONE 


ALT BANK 

CLEAR SHIFT COUNT 
RESULT NOW IN D’AH’L' 
CHECK FIRST RESULT BIT 
MSBIT IS ZERO 
INC SHIFT COUNT 
SHIFT RESULT RIGHT 


ROUNDING BIT IN CY 




0ROUND CALL RSTO 


; ROUND AND STORE FRACTION 


? 


5 SUBTRACT EX 

:ponents 


T 

SUB 

c 


JP 

MCHK 


f 

: RECIPROCAL 



T 

RFC CALL 

ST 01 

SAVE X 

LD 

HI- , Cl 


CALL 

LD 

: GET A 1 

I.D 

IYiTEMPI 


CALL 

DIVI 


RET 



*H MULT , DIVI 

UTILITIES 


T 

: RSTO ROUNC 

i, STORE FRACTION, LOAD EXPONENTS 

: INPUT: 

ROUNDING BIT IN CARRY 

J 

NORMALIZED 

FRACTION IN AHL 

r 

FINAL SIGN 

BIT IN A' 

OUTPUT: 

UPDATED SHIFT COUNT IN B 

f 

UNBIASED EXP'S IN A AND C 

TIMING: 

APPROX 186 


T 

RSTO LD 

DE , 0 


ADC 

HL , DE 

; ADD ROUNDING BIT 

ADC 

A, 0 


JR 

NC , $ + 3 

? NO OVERFLOW 

INC 

B 

J INC SHIFT COUNT 

; NOTE: DUE 

TO IMPLIED 

1, ACTUAL. SHIFTING UNNECE 

LD 

(IX+i),H 

; STORE FRACTION 

LD 

(IX+2),L 


AND 

7FH 

; TURN OFF IMPLIED 1 

LD 

H, A 

J SAVE REST OF BYTE 

EX 

AF,AF' 

? GET FINAL SIGN 

OR 

H 

? COMBINE 

LD 

(IX ) , A 

; STORE MSBYTE 

LD 

A,(IY+3) 

? GET EXP 

SUB 

BOH 

; REMOVE BIAS 

LD 

C , A 


LD 

A,(IX+3) 

; GET OTHER EXP 

SUB 

8 OH 

; REMOVE BIAS 

RET 




LOADY LOADS FRACTION PART OF (IY) INTO CDE 
INPUT: M8BYTE OF X IN A 

OUTPUT: SIGN OF RESULT IN BIT 7 OF A 

TIMING: 86 


LOADY 

LD 

C, (IY > 

? LOAD FRACTION 


LD 

D,(IY+1) 



LD 

E,(IY + 2) 



XOR 

C 

? CALC SIGN BIT 


AND 

80H 



SET 

7, C 

? RESTORE IMPLIED 1 

? 

RET 




5 OVER OR UNDERFLOW ERRORS 



M ERR 


MERR JP 

P ,UNDER ? 

UNDERFLOW 

EX 

AF , AF' ? 

GET SIGN 

LD 

( I X ) r A ; 

SET SIGN OF RESULT 

JP 

OVER 


1 

; DIVISION 

BY ZERO ERROR 


DIVO LD 

HL,ERROR 


SET 

4, (HD 


JP 

INF 


*H FIX 



TTrTTTTTTrr 

yyyyyyyyyyyyyy 

yyyyyyyyyyyyyyyyyyyyyyyyyyyyy 


FIX CONVERTS FLOATING POINT NUMBER TO 
16 BIT 2'S COMPLEMENT 
RESULT IS ROUNDED 

INPUT; IX PTS TO FL PT NUMBER 
OUTPUT; IX PTS TO 2'S COMP INTEGER 

HI. CONTAINS 2'S COMP INTEGER 
ERRORS: OVERFLOW IF INPUT CANNOT BE REPRESENTED IN 
16 BITS. MAX OUTPUT IS 32.767, MIN IS -32768 
ALGORITHM; SHIFT COUNT IS DERIVED FROM EXPONENT 
FRACTION PART IS SHIFTED RIGHT AND CONVERTED 
TO 2'S COMPLEMENT. 


} y 

yyyyyyyyyy 

********** 

7777777777 

*************************** 
777777777777777777 7 77777777 

FI 

X LD 

IY i C _ 5 



CALL 

ADD 

ADD 1/2 TO ROUND 


CALL 

FIXA 

? CONVERT NUMBER 


JR 

NC , FX60 

; NO OVERFLOW 

r 

? 

INTEGER OVERFLOW 


r 

LD 

HL,ERROR 



SET 

6, (HL) 

; FLAG OVERFLOW 


L.D 

HL. , 7FFFH 

, LARGEST INTEGER 


BIT 

7, (IX) 



JR 

Z i F X 6 0 

; SIGN POSITIVE 

* 

INC 

HL 

f SMALLEST INTEGER 

r 

y 

STORE RESULT 


y 

FX60 I..D 

( IX > , L. 



LD 

(IX + 1) , H 


t 

RET 



y 

FIXA MAIN 

16 BIT FIX 

OPERATION 

y 

INPUT; 

IX PTS TO 

FL PT NUMBER 

y 

OUTPUT; 

HL IS 16 BIT 2'S COMPLEMENT 

y 


CY IS SET 

IF HL IS INVALID (OVERFLOW) 

y 

FIXA XOR 

A 



LD 

H , A 

; CLEAR ACCUM 


LD 

L» A 



C P 

(IX+3) 



RET 

Z 

5 INPUT WAS ZERO 


LD 

H 7 (IX) 

; LOAD MS 16 BITS 


SET 

7 , H 

; RESTORE IMPLIED 1 


LD 

1.7 (IX+1) 




« 

i 

f 

t 

« 

f 

« 

« 

t 

« 

t 

€ 

« 

t 

i 

« 

f 



LD 

A r 16 + 8 0 H 



SUB 

(IX + 3) 

; CALC SHIFT COUNT 


RET 

c 

; X >= 2+16 


JR 

Z.FX25 

; COUNT IS ZERO 

y 

J INITIALIZE 

FOR LOOP 


T 

LD 

Br 16 

; SET MAX COUNT 


CP 

B 



JR 

NCrFXlO 

; COUNT >= 16 


LD 

B r A 

} SET COUNT 

F XI 0 

XOR 

A 

! CLEAR A 


CP 

( IX + 2 ) 



RLA 


; A IS 1 IF ( IX+2 )'> 

r 

; SHIFT 

t 

RIGHT 

' UNTIL C 

0UNT = 0 

r 

FX20 

SRL 

H 



RR 

L 



ADC 

A r 0 

; ACCUM BITS SHFTD 


DJNZ 

FX20 


r 

; CONVERT TO 

2’S COMPLEMENT 

T 

F X 2 5 

BIT 

7,(IX) 



JR 

Z ? FX40 

POSITIVE 


AND 

A 



JR 

Z ! FX30 

; EXACT INTEGER 


INC 

HL, 


FX30 

EX 

DE,HL 



LD 

HL ? 0 



AND 

A 



SBC 

HL,DE 

? NEGATE 


JP 

PrFXSO 

} OVERFLOW 


AND 

A 

? CLEAR CY 


RET 



FX40 

AND 

A 

; CLEAR CY 


BIT 

7 ? H 



RET 

Z 

? NO OVERFLOW 


; OVERFLOW 

T 

FK50 


SCF 

RET 


SET CY 


*H FLOAT 

rrrry?'r?yy*yrryfrrrrrr*?rryr*rrr*rr*rr*y'FrrF9’rr?rrrr?r7)’r)’a'rrrrr?irr*r)'r 

* 

* 

? FLOAT CONVERTS 16 BIT 2’S COMPLEMENT INTEGER 
J TO FLOATING POINT FORMAT 

; (FLOAT1 ASSUMES INPUT IS IN HL) 


INPUT 


IK 


PTS TO 16 BIT 2'S COMP INTEGER 
; OUTPUT; IK PTS TO FLOATING PT NUMBER 

; ERRORS; NONE 

; ALGORITHM; CONVERT TO SIGN MAGNITUDE AND NORMAL 
; REGISTERS; HL IS 16 BIT ACCUM 

; A IS SIGN 

; C IS EXPONENT 

T 

t Ml? j It t t ? r ? r t ? ? I yyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyy 


IZE 


yyyyyyyyyyyyyy 


€ 

i 

< 

€ 

€ 

I 

1 

I 

i 

i 

4 

4 

4 

4 

I 

4 

I 

4 

4 

4 

4 


t 


4 





f 

r 




1 


; TEST INPUT 





II 

r 

FLOAT ID 

L » (IX) 

; LOAD 16 BIT INTEGER 


1 

LD 

Hr(IX + 1) 





FLOAT! LD 

ArH 

5 TEST FOR ZERO 



ft 

OR 

L 



1 


jp 

Z, ZERO 

; INPUT WAS ZERO 



• 

r 

; CONVERT TO 

SIGN MAGNITUDE, SIGN IN A 


1 


r 

LD 

ArH 

? GET SIGN BIT 



• 

AND 

BOH 



1 


JR 

Z , FL..10 

; POSITIVE 




EX 

OErHL 




« 

LD 

HL , 0 



i 


SBC 

HI., DE 

; NEGATE 



ft 

r 

} NORMALIZE 
? 




i 


r 

FL..10 LD 

Cr16+SOH 

; INITIAL EXPONENT 



ft 

JP 

M , FL30 

} ALREADY NORMALIZED 


i 


AND 

A 

J CLEAR CY 



ft 

FL2Q DEC 

c 

; DEC EXP 



ADC 

HL., HL 

5 SHIFT LEFT 


i 


JP 

P , FL20 

; NOT NORMALIZED YET 



ft 

T 

J STORE RESULT 



i 

ft 

r 

F 1.30 RES 

7 , H 

J CLEAR IMPLIED 1 



OR 

H 

; ADD SIGN 


i 


LD 

(IX ) , A 





LD 

(IX+1) , L 




ft 

LD 

< IX+2),0 



i 


LD 

(IX + 3 ) r C 




ft 

RET 





BRK 




i 


; BRK SEPARATES INTEGER 

AND FRACTION PARTS 


i 

i 

? INPUT? 

IX PTS TO 

X 


m 


r OUTPUT; 

IX PTS TO 

FRACTION 0 <= X < 1 



ft 

r 

HL IS LOW 

ORDER 16 BITS OF 2'S COMP 

FIX32 ( X > 

1 


T 

CY IS SET 

IF INPUT > 2+24 




! ALGORITHM; USES 32 

BIT FIX AND FLOAT, AND Ft 

PT SUB 


ft 

f ALWAYS RETURNS 

VALID FRACTION PART 


1 


r 

BRK CALL 

ST 01 

J SAVE COPY OF X 



ft 

CALL 

FIX32 

J GET 32 BIT INTEGER 


1 


JR 

C , BRK20 

t v \ “7 .f. n L 




EXX 


ALT BANK 




BIT 

*7 p 



1 

IT 





F XX 


; MAIN BANK 




PUSH 

HL 

i SAVE HL AND FLAGS 



ft 

PUSH 

AF 



1 


CALL 

FL.032 

5 CONVERT TO FLOATING 




LD 

IY,TEMPI 




• 

CALL 

RSUB 

,’ CALC X - INT ( X ) 


1 


POP 

AF 

; RESTORE REGS 



ft 

POP 

HL 



1 

f 





i 





POSITIVE SIGN 


€ 


« 

« 

t 

€ 

€ 

€ 

f 

i 

i 

« 

i 

f 

t 

i 

€ 

C 

i 

i 



JR 

Z, BRK 


EX 

DE,HL 


LD 

HL , 0 


AND 

A 


SBC 

HL , DE 

BRKIO 

AND 

RET 

A 

BRK 20 

*H CP 

CALL 

SCF 

RET 

ZERO 


; NEGATE HL 
; CLEAR CY 

; FORCE ZERO FRACTION 
i SET CY 


? COMPARE ROUTINE 


r 

RETURNS 

CY SET IF 

X < Y 

T 

RETURNS 

Z SFT IF 

X = Y 

r 

CP 

LD 

A f (IX) 

? GET MS BYTE OF X 


AND 

80H 



JR 

NZ , CPIO 

; X IS NEGATIVE 


BIT 

7 , (IY) 



RET 

NZ 

i Y IS NEGATIVE 


LD 

A, ( IX + 3) 

; TEST EXP'S 


CP 

(IY+3) 



RET 

NZ 



LD 

A,(IX) 

J TEST MS BYTE 


CP 

(IY) 



RET 

NZ 



L.D 

A, ( IX + 1 ) 



CP 

( IY+1 ) 



RET 

NZ 



LD 

A, < IX+2) 



CP 

( IY + 2 ) 



RET 



CPiO XOR 

( IY) 

; TEST SIGN AND SET NZ 


RLA 


; PUT SIGN IN CY 


RET 

c 

; SIGNS DIFFERENT 


LD 

A t (IY + 3) 



CP 

( IX + 3 ) 



RET 

NZ 



LD 

A i ( IY ) 



CP 

(IX) 



RET 

NZ 



LD 

A r (IY + 1 ) 



CP 

(IX+1 ) 



RET 

NZ 



LD 

A,(IY+2) 



CP 

(IX+2) 



RET 



& H 

FLO 3? 



T 

t 

FL032 CONVERTS 32 

BIT INTEGER TO FLOATING POINT 

} 

BY NORMALIZING AMD 

ROUNDING TO 24 BITS 

? 

ENTRY TO 

FL.0321 WITH DE ZERO AVOIDS ROUNDING 

r 

INPUT 5 

H'L'HL IS 

32 BIT INTEGER 

t 


DE IS ROUNDING VALUE FOR FL.0321 

r 


B' BIT 7 

IS SIGN 


OUTPUT; 

IX POINTS 

TO FL PT NUMBER 


« 

t 

« 

< 

< 

i 

i 

i 

t 

« 

i 

« 

i 

l 

i 

i 

1 


I 


; TEST FOR ZERO 

« 


i 


€ 


€ 




FL032 

LD 

DE.80H 

FL.0321 

LD 

A, L 


OR 

EXX 

H 


OR 

L 


OR 

H 

9 

JP 

Z i ZERO 

r 

! NORMAL 

. IZE 

ACCUM. COU 

r 

F NORM 

LD 

C» OAOH 


BIT 

7 r H 


JR 

NZ,FROUND 

FNORL 

DEC 

EXX 

c 


ADD 

EXX 

HL r HI.. 


ADC 

HL.,HL 


JP 

P,FNORL 


? ROUNDING BIT 
? OR ALL 4 BYTES 

; ALT BANK 


; INPUT IS ZERO 

EXP IN C' 

; I NIT EXP TO +32 

? ALREADY NORMALIZED 
J DEC EXP 
i MAIN BANK 

SHIFT ACC LEFT 
; ALT BANK 

; NOT NORMALIZED YET 


ROUND FRACTION TO 24 BITS 

NOTE? ROUNDING OVERFLOW IMPLIES FRACTION 

IS A POWER OF 2, SO HIGHEST BIT NEED 
NOT BE PUT IN BEFORE STORING 


FROUND 

EXX 



ADD 

HL. r DE 


EXX 



LD 

DE i 0 


ADC 

HL , DE 


JR 

NC r FSTO 


INC 

c 

9 

JP 

Z.OVERB 

r 

; STORE 

RESULTS (ENTER 

t 

FSTO 

LD 

(IX + 3),C 


LD 

A i H 


AND 

7FH 


OR 

B 


LD 

(IX)iA 


LD 

(IX+1) r L 


EXX 



LD 

(IX + 2 ) t H 


RET 


7 

; LOAD 

FRACTION PART (C 

7 

F LD 

LD 

H. 0 


LD 

L r ( IX) 


LD 

ArL 


AND 

boh 


LD 

B , A 


SET 

7 t i. 


EXX 



LD 

H,(IX+1) 


LD 

L r (IX + 2) 


RET 


*H F1X32 



f MAIN BANK 
; ADD ROUNDING BIT 
! ALT BANK 


; NO OVERFLOW 
? INC EXP 
; OVERFLOW 

IN ALT BANK) 

; BINARY EXPONENT 
J PUT IN SIGN BIT 


; MSBYTE 
} 2ND BYTE 
? MAIN BANK 
J 2RD BYTE 


(CALLED IN ALT BANK) 

; H'L’HL IS ACCUM 

5 GET SIGN BIT 

l PUT IN B’ 

) RESTORE IMPLIED 1 
I MAIN BANK 



! FIK32 CONVERTS Ft PT NUMBER TO 32 BIT INTEGER 
i INPUT: IX POINTS TO Ft PT NUMBER 

i OUTPUT: H ’ L ’ HI... - 32 BIT INTEGER 

f B' ~ BIT 7 IS SIGN BIT 

; A -- ZERO IF INPUT WAS INTEGER 


T 


CY - SET 

IF INPUT >= 2+32 

7 

FIX32 

XOR 

A 

7 

TEST FOR ZERO 


CP 

( IX + 3 ) 




JP 

NZrFXA 

7 

NOT ZERO 


LD 

H ; A 

r 

CLEAR H'L' 


LD 

L 7 A 




L D 

B,A 

% 

r 

CLEAR SIGN 


EXX 


r 

MAIN BANK 


L D 

Hr A 

7 

CLEAR HL 


LD 

L r A 




RET 


7 

A IS ZEROr CY OFF 

; CALCULATE 

SHIFT COUNT 



r 

F XA 

CALL 

FLD 

y 

LOAD FRACTION TO L ’ HI. 


LD 

A r 98H 

T 

CALC SHIFT COUNT 


SUB 

(IX+3) 




RET 

c 

J 

INPUT WAS TOO BIG 


RET 

z 

7 

SHIFT COUNT IS ZERO 


LD 

B. 24 

7 

MAX COUNT IS 24 


CP 

B 




JR 

NC.FX1 

r 

COUNT >= 24 


LD 

B r A 

7 

PUT COUNT IN B 

F XI 

XOR 

A 

T 

CLEAR A 

f 

; SHIFT 

LOOP 




F XL 

EXX 


T 

ALT BANK 


SRL 

H 

T 

SHIFT ACCUM RIGHT 


RR 

L 




EXX 


7 

MAIN BANK 


RR 

H 




RR 

L 




ADC 

ArO 

y 

COUNT i’S SHFTD OUT 


DJNZ 

FXL 

r 

COUNT NOT DONE 

r 

: USE c 

EI.LI.NG FUNCTION 

FOR NEGATIVE NUMBERS 

T 

RET 

z 

y 

NO I’S SHIFTED OUT 


EXX 


1 

T 

ALT BANK 


BIT 

7 y B 




EXX 


r 

MAIN BANK 


RET 

z 

r 

POSITIVE SIGN 


L D 

DE» 1 

7 

INCREMENT ACCUM 


ADD 

HL t DE 




RET 

NC 

* 

t 

NO CY FROM INCR 


EXX 


r 

ALT BANK 


INC 

HL 

7 

ADD CY 


EXX 


T 

r 

MAIN BANK 


AND 

A 

T 

TURN OFF CY 


RET 




*H UTILITIES 








; UTILITIES 


r 

ST 0.1. 

I..D 

DE,TEMPI 


JP 

STO 

ST 02 

I..D 

DE i TEMF'2 


JP 

STO 

ST 03 

LD 

DE » TEMP3 

STO 

PUSH 

IX 


POP 

HL. 

MOV 

LD 

LOIR 

RET 

BC, 4 

ID 

PUSH 

IX 


POP 

DE 


JP 

MOV 


*H BOOLEAN 
*H BOOLEAN 


r 


? LOG 

ICAL OPERATORS 

r 

I OR 

LD 

A, (IX) 


OR 

(IX + 1 ) 


OR 

(IY) 


OR 

(IY+1 ) 


JP 

INE1 

I AND 

LD 

A,(IX ) 


OR 

(IX+1) 


JP 

Z,IZERO 


L D 

A,(IY) 


OR 

(IY + 1) 


JP 

INEi 

I NOT 

LD 

A,(IX ) 


OR 

(IX+1 ) 


JR 

IEQ1 

OR 

LD 

A , ( IX + 3 ) 


OR 

(IY+3) 


JP 

NE1 

AND 

LD 

A , (IX+3) 


AND 

A 


J P 

Z»ZERO 


LD 

A,(IY+3) 


AND 

A 


JR 

NE1 

NOT 

LD 

A,(IX+3) 


AND 

A 


JR 

E 01 


T 


} RELATIONAL 

OPERATOR 

SGT 

CALL 

SCP 


JR 

IGT 1 

SLE 

CALL. 

SCP 


JR 

ILE1 

SL.T 

CALL 

o L F 


JR 

ILT1 

SGE 

CALL 

SCP 


JR 

IGE 1 

SEQ 

CALL. 

SCP 


JR 

IE Q 1 

SNE 

CALL 

S C P 


; DEST IS TEMPI 

; DEST IS TEMP2 

; DEST IS TEMP3 
f SOURCE IS ACCUM 

; COUNT IS 4 

; DEST IS ACCUM 


; X IS ZERO 




JR 

INE 1 

IGT 

CALL 

I CP 

IGT1 

JR 

Z,IZERO 


JR 

Ci IZERO 


JR 

I ONE 

I L E 

CALL 

I CP 

I LEI 

JR 

Z t ION E 


JR 

C.ZONE 


JR 

I ZERO 

ILT 

CALL 

I CP 

ILT1 

JR 

Cr IONE 


JR 

I ZERO 

IGE 

CALL 

TCP 

IGEI 

JR 

C , IZERO 


JR 

IONE 

IEQ 

CALL 

I CP 

IE Q1 

JR 

ZrIONE 


JR 

IZERO 

INE 

CALL 

I CP 

I NE1 

JR 

Z, IZERO 


JR 

IONE 

GT 

CALL 

CP 

GT1 

JR 

Z t ZERO 


JR 

C t ZERO 


JR 

ONE 

LE 

CALL 

CP 

LEI 

JR 

Z ,ONE 


JR 

CrONE 


JR 

ZERO 

LT 

CALL 

CP 

L T1 

JR 

CrONE 


JR 

ZERO 

GE 

CALL 

CP 

G E1 

JR 

C r ZERO 


JR 

ONE 

EG! 

CALL 

CP 

EQ1 

JR 

ZrONE 


JR 

ZERO 

NE 

CALL 

CP 

NE.1, 

JR 

Z r ZERO 


JR 

ONE 

* H CONS 

TANT 

FUNCTIONS 

t 

l PRODUCE UNITY RESULT 

t 

I ONE 

LD 

( I X ) r 1 


LD 

( I X + 1 ) r 0 


RET 


ONE 

LD 

( I X ) r 0 


LD 

(IX+i) ,0 


LD 

( I X + 2 ) y 0 


LD 

(IX + 3)r81H 

T 

T 

RET 



; PRODUCE ZERO RESULT 


r 

IZERO 

LD 

o 

X 

i~j 


LD 

(IX + i) y 0 


RET 


ZERO 

LD 

( I X > r 0 


I.D 

l.D 

LD 

RET 


<I X +1 ) ,0 
(IK+2) ,0 
(IX+3) .0 


PRODUCE INFINITE RESULT 


INF 

RL 

(IX) ; GET OLD SIGN 


LD 

(IX)rOFEH 


RR 

(IX) 5 REPLACE SIGN 


LD 

(IX+1) i OFFH 


L D 

(IX+2)iOFFH 


LD 

(IX + 3)* OFF H 


RET 


t 

! FLAG 

OVERFLOW ERROR , STORE INFINITY 

r 

OVERB 

EXX 

J ALT BANK 


LD 

(IX)rB i SET SIGN BIT 

OVER 

LD 

HL,ERROR 


SET 

hr (HI.) J FLAG OVERFLOW 


JP 

INF 

T 

; FLAG 

UNDERFLOW ERROR> STORE ZERO 

LINDER 

LD 

HL,ERROR 


SET 

5,(HL) ; FLAG UNDERFLOW 


JP 

ZERO 

*H CONSTANTS 

AND WORK AREAS 

Cl 

DEFW 

0 ; 1,0 


DEFW 

8.100H 

C... 5 

DEFW 

o ; 0 . 5 


DEFW 

8000H 

SEED 

DEFW 

0 


DEFW 

0 

TEMPI 

DEFS 

4 

TEMP2 

DEFS 

4 

TEMP3 

DEFS 

4 


REFERENCES 


EXTERNAL 

GLOBAL 

GLOBAL 

GLOBAL 

GLOBAL 

GLOBAL 

GLOBAL 

GLOBAL 

GLOBAL 

GLOBAL 


ERROR,ICP,SCP 

ADD,SUB,RSUB»MULT,DIVI,RE C 

FIX,FIXA,FL OAT,BRK,FIX32,FL032,FL 0321 

CP,STO1,ST02 ,ST03,TEMP1,TEMP2 

TEMP3 , SEED , LD , FL.D, FLOAT 1 

OVER,UNDER,ZERO,ONE,INF 

GT,LT f GE » LE » EQ,NE »AND,OR,NOT 

IGT,ILT,IGE»ILE,IEQ,INE,IAND/IOR,INOT 

SGT rSLT,SGE,SLE,SEQ,SNE 

Cl,C_5 







# *H DECIMAL TO BINARY 


y y t y y 


€ 


€ 

« 

i 

€ 

i 

i 

f 

f 

< 

€ 

i 

i 

f 

t 

< 

i 

« 

i 

« 


DECIMAL TO BINARY CONVERSION 

INPUT; IK POINTS TO OUTPUT AREA 
OUTPUT; IX POINTS TO FL PT NUMBER 
ERRORS; NO DIGITS IN FRACTION PART 
NO DIGITS IN EXPONENT 
EXPONENT OVERFLOW 
EXPONENT UNDERFLOW 
REGISTER ASSIGNMENT IN MAIN LOOP: 

H’L'HL - 32 BIT ACCOM, I NIT ZERO 
B - DECIMAL EXPONENT, INIT 0 

C - FLAGS: BIT 6 - DIGITS FOUND 

BIT 7 - FRACTION OVRFLW 


f 

A 

- NEW DIGIT 

r 

t 

B' 

- SIGN 

IN 

BIT 7 

* f V t 1 V V 

T 7 7 7 7 7 7 

t r J t t ? 

TTrtTfTTTJ 

t y r 

yyyyyyyyyyyyyyyyyyy 

DTOB 

CALL 

LMARK 

y 

MARK SCAN HERE 


CALL 

DBINIT 

; 

INITIALIZE REGISTER 


CALL 

LIN 

y 

GET FIRST CHAR 


CP 

r * » 




JR 

Z, DIO 

y 

SIGN IS + 


CP 

T T 




JR 

N Z , D 2 0 

r 

NO SIGN 


EXX 


/ 

ALT BANK 


SET 

7 , B 

y 

FLAG NEGATIVE SIGN 


EXX 


y 

MAIN BANK 

T 

; LOOP 

TO PROCESS CHARS 

BEFORE DECIMAL POINT 

r 

DIO 

CALL 

LIN 

y 

GET NEXT CHAR 

D20 

CP 

7 T 

xs 




JR 

Z , D 4 0 

IT 

7 

DECML PT FOUND 


CALL 

DIGIT? 




JR 

C, D50 

t 

NOT A DIGIT 


SET 

6 , C 

y 

FLAG DIGIT FOUND 


BIT 

7 r C 




JR 

NZ,D30 

y 

OVERFLOW 


CALL 

MAD 

y 

CONVERT THIS DIGIT 


JP 

DIO 



D30 

INC 

B 

t 

INC DECML EXP 


JP 

DIO 



T 

; loop 

TO PROCESS CHARS 

AFTER DECIMAL POINT 

D 40 

i • A L L 

LIN 

y 

GET NEXT CHAR 


CALL 

DIGIT? 




JR 

C, D50 

y 

NOT A DIGIT 


SET 

6 r C 

y 

FLAG DIGIT FOUND 


BIT 

7 , C 




JR 

NZ.D40 

n 

r 

OVERFLOW 


DEC 

B 

r 

DEC DECML. EXP 


CALL 

MAD 

7 

CONVERT THIS DIGIT 


JP 

D40 




i CHANGE FRACTION TO FLOATING POINT FORMAT 

* 

T 

D50 BIT 6, C 


**%»*%»**«» 


77777777777 


t 

€ 

t 

€ 

i 

€ 

€ 


f 


i 

t 

i 

« 

« 

€ 

« 

i 

< 

i 


f 




NO DIGITS FOUND 
SAVE THIS CHAR 


JP 2 ,!REST 
PUSH AF 
CALL FL.032 
POP AF 


PROCESS EXPONENT 
REGISTERS? 


1 

C - 

BIT 7 IS SIGN 


T 


BIT 6 IS DIGIT 

-FOUND FLAG 

1 


BIT 5 IS DIGIT 

-OVERFLOW FLAG 

r 

D - 

IS ACCUM FOR D 

i-TO-B CONVERSION 

i 

LD 

D,0 ; 

INITIALIZE REGISTERS 


LD 

C j D 



CP 

’ E ' - ’ 0 ’ 



JR 

NZ»D90 5 

NOT AN ’E’ 


CALL 

LIN 



CP 

' + ' 



JR 

r* p 

Z t D60 ; 

r r 

SIGN IS + 


JR 

NZ r D70 ; 

NO SIGN 


SET 

7 i C i 

FLAG MEG SIGN 

D 60 

CALL 

LIN 


D70 

CALL 

DIGIT? 



JR 

C > D80 ? 

NOT A DIGIT 


SET 

6 t C r 

FLAG DIGIT FOUND 


B IT 

5»C 



JR 

NZ > Dioo ; 

OVER OR UNDERFLOW 


CALL 

MAD! t 

CONVERT THIS DIGIT 


JP 

D60 


D80 

BIT 

6 , C 



JR 

Z t D12 0 ; 

NO DIGITS FOUND 


BIT 

7,0 



JR 

NZ,EXPERR ? 

OVER OR UNDERFLOW 

D90 

CALL 

DSCALE ? 

SCALE BY DECML EXP 


JP 

RTN 


T 

; F'lNl 

SH SCANNING DIGITS 


r 

D 100 

CALL 

LIN 

FINISH SCANNING DIGIT 


CALL 

DIGIT? 



JR 

NCfDlOO ? 

DIGIT IS FOUND 

r 

? DETERMINE 

FROM SIGN BIT 

WHETHER OVER OR UNDER 

r 

EXPERR 

BIT 

7 , C 



JR 

NZ,Di1o ; 

SIGN IS NEGATIVE 


CALL 

OVER ; 

FLAG ERROR 


JP 

RTN 


D .1.10 

CALL 

UNDER ? 

FLAG ERROR 


JP 

RTN 


7 

? ILLEGAL NUMBER 


r 

D 120 

LD 

HL,ERROR 



SET 

7,(HL) 



RET 



*EJECT 




RTN 

CALL 

LBACK ? 

BACKSPACE SCAN PTR 



; DIGIT? ROUTINE CONVERTS ASCII CHAR TO 0-9 VALUE 


RETURNS CY ON IF NOT A DIGIT 


T 

DIGIT? 

SUB 

'O’ 



RET 

c 

; WAS < 30H 


CP 

10 



CCF 


5 CY SET IF > 9 


RET 



t 





; DBINIT CLEARS REGISTERS H'L',B’,HL,BC,A 


/ 

DBINIT XOR 

A 

i CLEAR A 

LD 

Hr A 


L D 

L. A 


LD 

Br A 


EXX 


; MAIN BANK 

LD 

Hr A 


LD 

Lr A 


L.D 

Br A 


LD 

C I A 


RET 




I 

; MAD ROUTINE TO MULTIPLY ACCUM BY 10 AND ADD NEW DIGIT 

T 


MAD CALL 

MUL10 

} MULTIPLY BY 10 

EXX 


? MAIN BANK 

LD 

D r 0 

; ADD A TO ACCUM 

L.D 

E t A 


ADD 

HL , DE 


EXX 


! ALT BANK 

JR 

NCr*+3 

? NO CARRY FROM LOW ORDER 

INC 

HL 

; ADD CARRY 

LD 

A r H 

? CHECK OVERFLOW 

AND 

OFOH 


EXX 


; MAIN BANK 

RET 

“7 

4- 

J NO OVERFLOW 

SET 

RET 

7rC 

; FLAG OVERFLOW 

r 

MUL10 MULTIPLIES 32 

BIT ACCUM BY 10 

■r ACCUM IS 

H’L'HL 


; USES D'E' 

DE 


; CALLED IN 

t 

1 MAIN BANK 

r RETURNS ALT BANK 

f 

MUL10 ADD 

HL r HL 

J SHIFT ACC LEFT 

LD 

E, L 

; COPY INTO DE 

LD 

DrH 


EXX 


! ALT BANK 

ADC 

HL , HL 


LD 

E,L 


LD 

DrH 


EXX 


; MAIN BANK 

ADD 

HL rHL 

? SHIFT ACC LEFT 

EXX 


; ALT BANK 

ADC 

HL , HL 


EXX 


; MAIN BANK 

ADD 

HL , HL 

; SHIFT ACC LEFT 

EXX 


; ALT BANK 


LEFT 


ML , OE 


ADD 

EXX 

ADC 

RET 


HL. , DE 


; ADD DE TO ACC 
; ALT BANK 


; MADi MULTIPLY 8 BIT ACC BY 10, ADD A 
; SAVES RESULT IN D 


MADI 

PUSH 

AF 

7 

SAVE A 


LD 

A, D 

7 

GET OLD RESULT 


ADD 

A, A 

* 

y 

SHIFT LEFT 


ADD 

A,A 

» 

r 

SHIFT LEFT 


ADD 

A, D 

r 

ADD OLD VALUE 


ADD 

A,A 

r 

SHIFT LEFT 


POP 

ADD 

DE 

A, D 

7 

NEW DIGIT IN D 


LD 

AND 

D, A 

OFOH 

r 

STORE IN D 


RET 

Z 

* 

f 

NO DIGIT OVERFLOW 

*EJECT 

SET 

RET 

5, C 

7 

FLAG OVERFLOW 


1 

i DSCALE ROUTINE TO COMBINE DECIMAL EXPONENTS 
; AND MULTIPLY BY POWER OF 10 FROM TABLE 

! INPUT: C - BIT 7 IS SIGN OF DECML EXP 

: D - IS MAGNITUDE OF DECML EXP 

? B - IS 2'S COMPL EXP FROM FRACTION 

i OUTPUT: IX POINTS TO FL PT NUMBER 
; ERRORS: EXPONENT OVER OR UNDERFLOW 

; REGISTER ASSIGNMENT: 

: DE- POINTER INCREMENT 

; C - EXPONENT MAGNITUDE 

i B - LOOP COUNTER 

; IY- TABLE POINTER 


; INITIALIZE 

y 

DSCALE 


DSC 10 


LD 

IY,TOPTAB 

» 

y 

TABLE POINTER 

LD 

A , D 

7 

MOVE EXP MAGN TO A 

LD 

D E , - 5 

% 

y 

POINTER INCREMENT 

IT S 

IGN-MAGNITUDE 


TO BINARY AND ADD 

BIT 

7 r C 



JR 

Z, DSC 10 

y 

POSITIVE SIGN 

NEG 




LD 

C,0 

y 

CLEAR SIGN BIT 

ADD 

A , B 

r 

7 

ADD BOTH EXP'S 

RET 

Z 

y 

ZERO EXPONENT 

JP 

P, DSC40 

y 

RESULT IS POSITIVE 

JP 

PE, OVER 

y 

OVERFLOW 

NEG 


y 

GET MAGNITUDE 

CP 

100 



JP 

NC,UNDER 

y 

MAGNITUDE >= 100 

LD 

C , A 

y 

SAVE MAGN IN C 

LD 

B,LENTAB 

y 

LOOP COUNTER 



I LINEAR SEARCH THRU TABLE OF POWERS OF 10 


CP 

( IY + 4 ) 


JP 

C, DSC30 

5 EXP < TABLE EXP 

PUSH 

BC 

; SAVE REGS 

CALL 

DIVI 

; DIVIDE BY TABLE ENTRY 

POP 

BC 


LD 

DE i - 5 

} RESTORE DE 

LD 

ArC 


SUB 

(IY + 4 ) 

J SUB TABLE EXP 

L.D 

C , A 


ADD 

lYrDE 

; DECR POINTER 

DJNZ 

RET 

DSC 20 

; LOOP NOT DONE 


r 

; EXPONENT SIGN IS POSITIVE 


t 

r 

DSC 40 

JP 

PE,UNDER 

; UNDERFLOW 


CP 

100 



J p 

NC,OVER 

} EXP >= 100 


LD 

C, A 

J INIT FOR LOOP 


LD 

B,LENT AB 


T 

5 LINEAR SEARCH THRU 

TABLE OF POWERS OF 10 

r 

DSC 50 

CP 

(IY+4) 



JP 

C , DSC60 

? EXP < TABLE EXP 


PUSH 

BC 

; SAVE REGS 


CALL 

MULT 



POP 

LD 

BC 

DE , - 5 

? RESTORE DE 


LD 

A, C 



SUB 

(IY+4) 

; SUB TABLE EXP 


LD 

C, A 


DSC60 

ADD 

IY, DE 

i DECR POINTER 


DJNZ 

RET 

DSC50 

? LOOP NOT DONE 

*H BINARY TO 

DECIMAL 


i J r ? 7 f 

1f 

y f r f t r t 

V » t t 1 K t H 

7T*7T?rr 

*»*»*»»***»** **»»**»* 


i FLOATING POINT BINARY TO DECIMAL CONVERSION 
I RESULT IS ROUNDED TO 7 DECIMAL DIGITS 
J INPUT; HI.. POINTS TO OUTPUT AREA 

} IX POINTS TO FL PT NUMBER 

? OUTPUT; ASCII STRING IN OUTPUT AREA 

i OUTPUT; FL PT NUMBER IS MODIFIED 

; A IS 2’S COMPLEMENT EXPONENT C-128<=A<=126) 

i OUTPUT AREA IS 12 BYTES; 

! 1 SIGN OF FRACTION (SPACE OR -•) 

5 7 DIGITS OF FRACTION 

; 1 " E" 

? 1 SIGN OF EXPONENT < + OR ~) 

! 2 DIGITS OF EXPONENT 

; EXPONENT IS FOR DECIMAL PT AFTER FIRST DIGIT 

; ERRORS; NONE 

f REGISTER ASSIGNMENT FOR INNER LOOP: 

5 H'L'HL - 32 BIT ACCUM, INIT FRACTION PART 

; WITH 4 BITS TO LEFT OF POINT 

i C - DECIMAL. EXPONENT 

? 8 - LOOP COUNTER 


IY - OUTPUT POINTER 


7777777777 


7777777T7T777777777 7 7777777777 


; CONVERT SIGN, TEST FOR ZERO 


7 

BTOD 

BIT 

7, (IX) 

7 

TEST SIGN BIT 


LD 

(HL> , ' ' 

7 

DEFAULT SIGN 


JR 

Z,BD02 

7 

POSITIVE SIGN 


LD 

( HL. ) r ' - ' 



BP 02 

INC 

HL 




PUSH 

HL 

x 

7 

SAVE OUTPUT POINTER 


LD 

A,(IX + 3) 

x 

7 

CHECK FOR ZERO 


0 R 

A 




JP 

NZ,BD05 

? 

X IS NOT ZERO 


LD 

HL , 0 

t 

7 

CLEAR ACCOM 


E XX 


7 

MAIN BANK 


LD 

HL, 0 




LD 

C,0 

? 

CLEAR DECML EXP 


POP 

IY 


OUTPUT PTR IN IY 


JP 

BD50 



7 

i SCALE 

X, LOAD FRACTION 



7 

BD05 

CALL 

BSCALE 

* 

7 

SCALE BY POWER OF 10 


EXX 


7 

ALT BANK 


CALL 

FLD 

7 

LOAD FRACTION 


POP 

IY 

7 

PTR IN IY 


7 

; SHIFT 

TO GET ZERO BINARY EXP 

? AND 

ALIGN BIN POINT BETWEEN BITS 3 AND 

; OF 

H’ 


; INITIALLY, 7D <> EXP <= 80 

7 

LD 

A,(IX+3) J GET BINARY EXP 


SUB 

7CH ; REMOVE BIAS, ADD 

BD20 

ADD 

HL,HL J SHIFT LEFT 


EXX 

l ALT BANK 


ADC 

HL , HL 


EXX 

; MAIN BANK 


DEC 

A ; DECR EXPONENT 

* 

JR 

NZ.BD20 ; EXP NOT ZERO YET 

7 

l NORMALIZE 

TO? .1 < ACCUM < 1 

; . l 

IS .199999 (HEX) 

? 

EXX 

J ALT BANK 


LD 

A,1 ? FIRST HEX DIGIT 


CP 

H 


JP 

NZ.BD40 ; ACCUM <> .1 


LD 

A,99H ? ALL OTHER HEX DIG 


CP 

L 


JP 

NZ,BD40 ; ACCUM <> .1 


EXX 

; MAIN BANK 


CP 

H 


JP 

NZ,BD45 ; ACCUM <> .1 


CP 

L 


JP 

C,BD50 ; ACCUM > .1 

BD30 

CALL 

. MUL10 r MPY BY 10 


EXX 

; MAIN BANK 


t 


DEC 

c 

r 

DECR DECHL EXP 

1 



JP 

BD50 





BD40 

EXX 


r 

MAIN BANK 


« 

B D 4 5 

JR 

NC,BD30 

T 

ACCUM < .1 

« 


7 

; ROUND 

BY ADDING 1/2 

LAST 

DECIMAL POSITION 


« 

r 





< 


BD50 

LD 

B t 7 

? 

LOOP COUNTER 




LD 

DE,ODH 

r 

.5E-7 

< 

1 


ADD 

HL , DE 

7 

ADD ROUNDING VALUE 



EXX 


* 

r 

ALT BANK 




J R 

NC,$+3 

T 

NO CY FROM LOW BYTES 

« 

« 


INC 

HL 

V 

t 

ADD IN CARRY 



L.D 

A r H 

7 

GET MSBYTE 




AND 

OFOH 




c 


JR 

Z,BD55 

7 

NO ROUNDING OVERFLW 

€ 



INC 

c 

V 

T 

INCR DECML EXP 




JP 

BD65 

7 

USE THIS AS 1ST DIGIT 

i 

« 

BD55 

C V V 

c: /\ 


t 

MAIN BANK 

i 

7 

; loop 

TO CONVERT FRAC 

TION 

m 

m. 

7 

BP 60 

CALL 

MUL10 

7 

MPY BY 10 

m 


B D 6 5 

LD 

A, H 

7 

GET MSBYTE 


f 


LD 

D, H 



i 



AND 

OFH 

7 

CLEAR TOP 4 BITS 




LD 

H, A 




« 


L D 

A, D 

7 

GET MSBYTE AGAIN 

i 



EXX 


7 

MAIN BANK 




RRA 


7 

SHIFT RIGHT 4 

< 

1 


RRA 






RRA 





f 


RRA 






AND 

OFH 

7 

CONVERT TO ASCII 

« 



OR 

'O' 






LD 

(IY),A 

7 

OUTPUT CHAR 

i 

f 


INC 

IY 





D JN2 

BD60 

7 

LOOP NOT DONE 


f 

7 

1 CONVERT EXPONENT 



i 


r 

LD 

( I Y ) , ’ E ' 




« 


LD 

( IY + 1 ) , ' + 

7 


« 



LD 

A , C 

7 

GET 2 ’ S COMP EXP 




PUSH 

AF 

t 

T 

SAVE IT 


< 


DEC 

A 

7 

SUB ONE 

l 



JP 

P , BD70 

f 

r 

SIGN IS + 




LD 

< IY + 1 ) , ' - 

7 



< 


NEG 




i 


BD70 

LD 

B , OFFH 

1 

? 

QUOTIENT 




LD 

C, 10 

r 

T 

DIVISOR 


i 

BD80 

SUB 

c 

* 

f 

LOOP DIVIDES BY 10 

i 



INC 

B 

7 

INCR QUOTIENT 




JR 

NC.BD80 

7 

STILL POSITIVE 


i 


ADD 

A,C 

7 

RESTORE REMAINDER 

i 



OR 

’0’ 

7 

CONVERT TO ASCII 




LD 

(IY+3) r A 

r 

STORE 2ND DIGIT 


« 


LD 

A , B 



i 



OR 

! O r 

7 

CONVERT TO ASCII 


« 


LD 

(IY+2) , A 

7 

STORE 1ST DIGIT 

i 


f € 




POP AF ! GET BINARY EXP AGAIN 

RET 

* EJECT 

7 

; BSCALE TRY TO REDUCE BINARY EXP TO ZERO 
J SCALING BY POWER OF 10 

! INPUT: IX POINTS TO FL PT NUMBER 

; OUTPUT: IX POINTS TO SCALED Fl. PT NUMBER 

: .0625 < X < 1 

J C CONTAINS DECIMAL EXP OF SCALE FACTOR 

: ERRORS: NONE 

? REGISTER ASSIGNMENT FOR LOOP: 

: A - IS COMPARISON EXPONENT 

J 8 - IS LOOP COUNTER 

: C - IS DECIMAL. EXPONENT 

? DE - IS POINTER INCREMENT 

! IY - IS TABLE POINTER 


BSCALE 

LD 

IY r TOPTAB 

IT 

r 


LD 

DE r -5 

* 

T 


LD 

B r LENTAB 

t 

T 


LD 

C, 0 

T 


LD 

A 7 (1X43) 

r 


CP 

81H 



JR 

NC,BSC30 

T 


NEG 


r 

7 

; SCAN 

T 

TABLE 

OF POWERS 

OF 

r 

BSC 10 

CP 

( IY 4 3 ) 



JP 

C,BSC20 

t 


PUSH 

BC 



CALL 

MULT 

If 

1 


POP 

BC 



L. D 

DE r -5 

« 

7 


LD 

A r C 

T 


SUB 

(IY+4) 

7 


LD 

C t A 

T 


X OR 

A 



SUB 

(1X43) 

7 

BSC 20 

ADD 

IY,DE 

r 


DJNZ 

BSC 10 

7 


RET 



7 

; SCAN 

TABLE 

OF POWERS 

U p 

T 

BSC30 

ADD 

A r 3 

t 


JR 

C,BSC50 

r 

B8C4D 

CF* 

(IY+3 ) 



JP 

C,BSC60 

7 

BSC50 

PUSH 

BC 



CALL 

DIVI 

7 


POP 

BC 



LD 

DE r -5 

t 

7 


LD 

A rL 

♦ 

r 


ADD 

A r (IY + 4) 

7 


L.D 

C r A 



LD 

A, (IX + 3) 

7 


ADD 

A r 3 

7 


PT TO TOP OF TABLE 
5 BYTES PER ENTRY 
LOOP COUNTER 
CLEAR DEC EXP 
GET BINARY EXP 

EXP IS POSITIVE 
NEGATE EXP 

10 


EXP < TABLE EXP 

MULTIPLY BY TABLE ENTRY 

POINTER INCREMENT 
GET DECML EXP 
SUB LOG(TABLE ENTRY) 
RPLACE DECML EXP 

GET NEW BINRY EXP 
DECR POINTER 
LOOP NOT DONE 


10 

ADD 3 TO BINRY EXP 
OVERFLW f USE LARGEST SCALE 

EXP+3 < TABLE EXP 

DIVIDE BY TABLE ENTRY 

POINTER INCREMENT 

GET DECML EXP 

ADD LOG(TABLE ENTRY) 

GET NEW BINRY EXP 
ADD 3 








D JNZ 

BSC 40 

5 LOOP NOT DONE 

RET 



*H RND 



i 77777777777 

77777777777 

777777777777777777777777 

r 

; RANDOM NUMBER GENERATOR 

J PRODUCE 

S 32 BIT RANDOM NUMBERS: 

; X(I+1> = 3141592621 * X(I) + 2718281829 

; NORMALIZED AND TRUNCATED TO 24 BIT 

; FLOATING POINT, 

0 <> X < 1 

; INPUT; 

IX POINTS 

TO OUTPUT AREA 

t 

t 

SEED,SEED+ 

2 CONTAIN LAST 32 BIT RE 

; output; 

IX POINTS 

TO RANDOM Ft. PT NUMBER 

i 

SEED,SEED+2 CONTAIN NEW RESULT 

; ERRORS; 

NONE 


; REGISTER 

ASSIGNMENT 

FOR MAIN LOOP: 

; H ' L ' 

HL. - IS 32 

BIT ACCUM 

; D'E' 

DE - IS 32 

BIT REGISTER 

; 8 ' C ' 

CA - IS 32 

BIT RIGHT EXTENSION 

; & 

- IS LOOP COUNTER 

t irrirririrr 

rrrrirrrrrr 

vvvtvvvvvvvvtvvvvvtvvtt* 
777777 7 77777777777777777 

r 

i INITIALIZE 

REGISTERS 


RND LD 

HL.,0 

; CLEAR ACCUM 

LD 

DE,(SEED) 

f GET X(I> 

LD 

BC,0BB40H 

; MS HALF OF 3141592621 

EXX 


; MAIN BANK 

LD 

HL , 0 


LD 

DE,(SEED+2 

) 

LD 

C , 0E6H 

5 3RD BYTE OF 314159262 

LD 

A, 2DM 

? 4TH BYTE OF 314159262 

LD 

B, 33 

; LOOP COUNTER 

JP 

ML 2 



32 BIT MULTIPLY LOOP 


ML.OOP JR 

NC,ML 1 ; 

MPLIER BIT IS ZERO 

ADD 

HL,DE ? 

ADD MPLICAND TO ACCUM 

EXX 


ALT BANK 

ADC 

Hl., DE 


EXX 

* 

r 

MAIN BANK 

T 

! ROTATE ACCUM RIGHT 
* 


? 

ML1 EXX 

t 

ALT BANK 

RR 

H 


RR 

L 


EXX 

f 

MAIN BANK 

RR 

H 


RR 

L 


F 

; ROTATE RIGHT EXTENSION 


r 

HL2 EXX 

y 

ALT BANK 

RR 

B 


RR 

c 



E XX 
RR 


e 


MAIN BANK 



t 


RRA 


7 

MPLIER BIT IN CY 



DVNZ 

MLOOP 

7 

LOOP COUNT NOT DONE 

m 

7 

7 

MOVE LOW ORDER 32 BITS 


TO ACCUM 


7 

EXX 


7 

ALT BANK 

m 


LD 

H, B 





L. D 

L., C 





EXX 


7 

MAIN BANK 

i; 


LD 

H, C 





LD 

L. , A 



• 

7 

7 

ADD 27.18281829 




T 

LD 

DE,0B065H 

7 

LS HALF OF 27182818 

€ 


ADD 

HL, DE 





EXX 


7 

ALT BANK 



LD 

DE,0A205H 

7 

MS HALF OF 27182818 

« 


ADC 

HL , DE 




7 

J 

STORE NEW 

RESULT 



t 

7 

LD 

(SEED),HL 

7 

MS HALF 



LD 

B , 0 

7 

CLEAR SIGN BIT 

i 


EXX 


7 

MAIN BANK 



LD 

(SEED + 2) , HI¬ 

7 

LS HALF 



L.D 

DE , 0 

J 

SET FOR NO ROUNDING 

t 


CALL 

FL.0321 

* 

7 

FLOAT INTEGER 



LD 

A,(IX + 3) 

7 

DIVIDE BY 2+32 



SUB 

"7 O 
v3 jC« 



( 


L..D 

(IX + 3 ) , A 





RET 





i 

« 

f 

i 

« 

t 

€' 

t 


RANDOMIZE USES REFRESH REGISTER TO MODIFY SEED 


RANDMZ LD A,R ; GET REFRESH COUNTER 

LD (SEED),A ; MODIFY SEED 

RET 

* H SHORT FUNCTIONS 

777777777777777777777777777777777777777777777777777777777777777777777777 


? INTEGER FUNCTION --- FLOOR FUNCTION 
; INPUT,OUTPUT: IX POINTS TO FL PT NUMBER 


INT CALL FIX32 

RET C 
CALL FL032 
RET 


i CONVERT TO INTEGER 
: X >= 2+24 

,' CONVERT TO FLOATING POINT 


ABSOLUTE VALUE FUNCTION 

INPUT,OUTPUT: IX POINTS TO FL PT NUMBER 


ABS RES 7,(IX) 

RET 


? TURN OFF SIGN BIT 


I 

i 

i 

< 

I 

i 

i 

i 

I 

€ 

I 

I 

i 

i 

l 

I 

« 

i 

l 

I 

I 


i 



; SIGN FUNCTION 

; INPUT,OUTPUT? IX POINTS TO FL PT NUMBER 

r 

TrfrrrrjryTyjfffyrtrrttTrttyiirrtttrrfrrjrrrtittiJtrtfrtriTrrrTTirirjt 

SGN XOR A J CLEAR A 

CP (IX+3) 

; ZERO INPUT IS NOP 
J GET SIGN 

; OUTPUT i 
i SIGN OF INPUT 

. . , . . .. . . , , . 

7y777r7f?7f77777777777777tS7777tfrrrt7r77ry7r7t7777r7trr7f 7 7777f7fT7rr 


XOR 

A 

CP 

(IX+3) 

RET 

Z 

LD 

A,(IX) 

AND 

80H 

CALL 

ONE 

LD 

RET 

(IX),A 


i NEGATE FUNCTION 


NEG LO A, (IX) 

XOR 80H 
LD (IX),A 
RET 

*H INTEGER ADD AND SUB 


! IADD, ISUB INTEGER ADD AND SUBTRACT 


* 

INPUT: IX 

F'TS 

TO 

X 

7 

IY 

PTS 

TO 

Y 

1 

OUTPUT: IX 

F'TS 

TO 

SUM OR DIFFERENCE 

t 

T 

ERRORS: OVERFLOW 

(POSITIVE OR NEGATIVE) 


1 7 ? ? 7 7 7 7 

t t ? f 

7 7 7 r 

777777777777 

777777777777777 

I ADD 

CALL 

ASLOAD 

} LOAD ARGUMENT 


ADC 

HL,, DE 



JP 

F'O, ASSTO 

; NO OVERFLOW 

ASOV 

vJP 

PfNINF 

? OVFL.W SIGN PO 


JP 

F'INF 


I SUB 

CALL 

ASLOAD 

? LOAD ARGUMENT 


SBC 

HL, DE 



JP 

PO,ASSTO 

; NO OVERFLOW 


JR 

ASOV 


7 

: ASLOAD 

t 

' LOADS ARCS TO 

HL, DE, RESETS 

7 

ASL.OAD 

LD 

H,(IX+1) 



LD 

L , (IX ) 



LD 

D, (IY+1 ) 



LD 

E, (IY) 



AND 

A 



RET 




} ASSTO STORES HL AS RESULT 

r 

ASSTO LD (IX),L 

LD (IX+i),H 

RET 

*H INTEGER MULTIPLY 



IMUL INTEGER MULTIPLY 


INPUT; IX PTS TO X 

IY PTS TO Y 

OUTPUT; IX PTS TO PRODUCT 

ERRORS; OVERFLOW (POS OR NEG) 

REGISTERS; DE IS MULTIPLICAND 

AC IS MULTIPLIER AND ACCI.JM EXTENSION 
HL IS ACCUM 
B IS LOOP COUNTER 


777777 

1 » » » t t 

t r r * i t 

777777777777 

rrttttTtrt9t77J7777ri7T 

IMUL 

CALL 

MDLOAQ 

; LOAD ARCS 


LD 

A r B 

5 X IN AC 

r 

J 16 BIT BY 

16 BIT UNSIG 

NED MULTIPLY 

r 

LD 

HL , 0 

; Ct EAR AC:CUM- - 


LD 

B, 16 

} SET CcJUnTER 


R'RA 


; INITIALIZE CARRY 


RR 

c 


I MIO 

JR 

NC , IM20 

; MF'LIER BIT IS ZERO 


ADD 

HL , DE 

; ADD MULTIPLICAND 

I M20 

RR 

H 

; ROT ACC-MPLICAND INTO 


RR 

L 



RRA 




RR 

C 



DJNZ 

IM10 

? COUNTER NOT 0 

r 

! TEST 

F 0 R > 

16 BIT PRODUCT 

r 

LD 

B t A 

; PRODUCT IN BC 


LD 

A r H 



OR 

L 



JR 

NZ»INF? 

; OVERFLOW 

? 

; TEST 

FOR > 

15 BITS, AND 

STORE RESULTS 

r 

HDSTO 

EX 

AFr AF’ 

; GET SIGN OF RESULT 


AND 

8 OH 



JR 

Z, MDS20 

? POS RESULT 

; NEGATE BC 



r 

LD 

HL 7 0 



SBC 

HL, BC 



JR 

Z, MDSIO 

1 RESULT WAS ZERO 


JP 

P»NINF 

} NEGATIVE OVERFLOW 

MDSIO 

LD 

(IX) , L 

; STORE RESULT 


LD 

(IX+1),H 



RET 



7 

; TEST 

FOR > 

15 BITS 


MDS20 

BIT 

7 t B 



JR 

NZrF'INF 

; POSITIVE OVERFLOW 


LD 

( IX) ,C 



LD 

(IX+i),B 



RET 




777777777777777777777 








< 

ft 

ft 

f 


i 

« 

i 

€ 

< 

I 

ft 

f 

< 

« 

€ 

« 

i 

f 

ft, 

t 


5 INF? CHECKS SIGN BIT IN A’ 
f PINE PRODUCES POS OVERFLOW 
? NINE PRODUCES NEG OVERFLOW 


I DVD 

LD 

HL , ERROR 



SET 

4, (HL) 

; FLAG ERROR 

INF? 

EX 

AF, AF ’ 

; GET SIGN BIT 


AND 

8 OH 



JR 

NZ,NINF 

; NEG SIGN 

F'lNF 

LD 

(IX),OFFH 

J RESULT IS 7FFFH 


LD 

(IX+1 ) , 7FH 



JR 

FINF 


NINF 

LD 

( IX ) , 0 

; RESULT IS 8000H 


LD 

(IX+1) , 80H 


FINF 

LD 

HL , ERROR 



SET 

6 , (HL ) 

i FLAG ERROR 


RET 



*H INTEGER 

DIVIDE 


; ; ?; 

» » t » t T 1 

7 7 7 7 7 7 7 

************ 

777777717777 

yyyyyyyyyyyyyyyyyyy 

r 

: IDIV INTEGER DIVIDE 


r 

r 

INPUT: 

IX PTS TO 

X 

* 

r 


IY PTS TO 

Y 

T 

OUTPUT : 

IX PTS TO 

X/Y 

T 

ERRORS: 

Y IS ZERO 



ALGORITHM: NON-RESTORING DIVISION, 16 BY 15 
REGISTERS: AC IS DIVIDEND, ACCI.JM EXT 

DE IS DIVISOR 

HL IS ACCIJM, REMAINDER AT END 


t V Hr ft ft 

7 7 7 7 7 

% % % ft ft % ft 

7 7 7 7 7 7 7 

yyyyyrryyy 

yyyyyyyyryyyyyyyyyyyyyy 

IDIV 

CALL 

MDL.OAD 

; LOAD ARGS 


LD 

A, D 



OR 

E 



JR 

Z,IDVO 

J DIVISOR IS ZERO 


LD 

A , B 

? X IN AC 

y 

; 16 

BIT BY 

15 BIT UNS 

IGNED DIVIDE 

T 

LD 

HL ,0 

; CLEAR ACCUM 


LD 

B , 16 

1 SET COUNTER 

IDV10 

RL. 

C 

i 0 <> HL < DE 

J ROT ACC-RESULT LEFT 


RLA 




ADC 

HL, HL 

; NO CY POSSIBLE 


SBC 

HL , DE 

? SUB DIVISOR 

IDV20 

CCF 

JR 

NC , IDV50 

; CALC RESULT BIT 
ACCIJM WENT NEGATIVE 

IDV30 

DJNZ 

IDV10 

J COUNTER NOT 0 


JP 

IDV60 


IDV40 

RL 

c 

; -DE <= HL < D 

J RO ACC-RESULT LEFT 


RLA 




ADC 

HL , HL 



AND 

A 

: TURN OFF CY 


ADC 

HL. , DE 

; ADD DIVISOR 


JR 

C, IDV30 

; ACCUM WENT POSITIVE 


JR 

Z , IDV20 

? ACCUM BECAME ZERO 

I DV50 

DJNZ 

IDV40 

5 COUNTER NOT 0 


< 

ft 

ft 

ft 

ft 

€ 

i 

€ 


ft 

< 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 


t 


€ 


w* 

IDV60 

RL 

C 

? SHIFT IN LAST RESULT 

BIT 

m 



RLA 







LD 

B, A 

; QUOTIENT IN BC 



9 


JP 

MDSTO 

; STORE RESULT 


€ 


7 

; HDLOAD LOADS X INTO 

BC, Y INTO DE 



t 

r TAKES ABS OF EACH, 

STORES X0R OF SIGNS 


« 


i IN 

A' 



j 


t 

7 

HDLOAD 

L.D 

B.(IX+i> 

? LOAD X 

i 

t 



LD 

C,(IX) 


j 




LD 

A, B 


l 


« 


AND 

8QH 



i 



JR 

Z, HDLIO 

; X IS POSITIVE 

if 

i 


« 

r 

NEGATE X 



i 

t 


T 

LD 

HL ,0 




€ 


SBC 

HI., B C 



€ 



LD 

B, H 


/ 




LD 

C, L 


f 


€ 

T 




j 

« 


; LOAD 

Y 





t 

y 

HDLIO 

LD 

D,(IY+1) 


1 

? 

i 



LD 

E,(IY) 






X OR 

D 

? FORM SIGN OF RESULT 



« 


EX 

AF,AF r 

; SAVE IT 


i 



BIT 

7, D 


/ 




RET 

Z 

? Y IS POSITIVE 



i 

y 

— .-. 



/ 

« 


? NEGATE Y 





« 

7 

LD 

HL , 0 



« 



SBC 

HL, DE 






EX 

DE , HL. 




t 


RET 




i 


*H IABS 

AND 

ISDN 




f 

y ? 7 r ? r y 

T 

7 7 7 7 7 

yyyyyyyyyyy 

yyyyyyyyyyyyyyyyyyyyyyy 

yyyyyyyyyyyyyyyyyyyyyyyyyy 

i 


IABS 






t 

T 

» ****** 

7 7 7 1 7 7 7 

t « * » « 

7 17 7 7 

*********** 

yyyyyyyyyyy 

yyyyyyyyyyyyyyyyyyyyyyy 

yyyyyyyyyyyyyyyyyyyyyyyyyy 

i 


IABS 

BIT 

7,(IX+1) 






RET 

Z 

; X IS POSITIVE 



( 

INEG 

LD 

D, <IX + 1 ) 

; LOAD X 


4 



LD 

E , (IX ) 






LD 

HL , 0 

; NEGATE X 



• 


AND 

A 



4 



SBC 

HL,OE 






JP 

PE,NINF 

; NEG OVERFLOW 



« 


LD 

(IX+1) , H 

} STORE X 


i 



L.D 

(IX) ,L 






RET 





9 

V *«««** 

1 7 7 7 7 7 1 

»«*,** 
17 111 

*********** 

7 7 777 7 77777 

, f t t t , • t 1 t , t t t 1 1 1 t 1 T t « t 

yyyyyyyyyyyyyyyyyyyyyyy 

yyyyyyyyyyyyyyyyyyyyyyyyyy 

i 


7 

i I y b N 






• 

7 





( 


******* 

* ? r ? r r r 

y y i i y 

yyyyyyyyyyy 

yyyyyyyyyyyyyyyyyyyyyyy 

*t»»***»»**»t»tt*t*,*,,,,t 

yyyyyyyyyyyyyyyyyyyyyyyyyy 



IS G N 

LD 

A,(IX + 1) 

; GET SIGN 




• i 


« 



i 



BIT 

7, A 


JR 

Z,ISG10 J 

X IS NONNEG 

LD 

(IX)»OFFH ? 

STORE ~1 

LD 

(JX+1)iOFFH 


RET 



I SGI 0 OR 

(IX) 


RET 

2. ; 

X IS ZERO 

LD 

(IX+l)»i 

STORE +1 

L.D 

(IX) ,0 


RET 



*H ICP 



y yyyyyyyyyyy 


?7ttT7?7TT7?7frfTf 

; ICP INTEGER COMPARE 


7 

; RETURNS 

CY SET IF X < 

Y 

? RETURNS 

Z SET IF X = 

Y 

yyyyyyyyyyyy 

yyyyyyyyyyyyy 

77777777777777777 

ICP LD 

B» (IX + 1) 

GET SIGN OF X 

L.D 

A, B 


AND 

80 H 

TEST SIGN r CLEAR 

JR 

NZr. ICP10 r 

X IS NEG 

B IT 

7.(IY + i ) 


RET 

NZ ? 

Y IS MEG 

L D 

ArB 


CP 

(IY+i) 5 

SIGNS ARE BOTH P 

RET 

NZ 


LD 

Ar(IX) 


CP 

(IY) 


RET 



I CP10 XOR 

(IY+i) ; 

TEST SIGN > CLEAR 

RLA 

7 

SIGN BIT INTO CY 

RET 

C 

SIGNS DIFFERENT 

LD 

A t B 


CP 

(IY+1) J 

BOTH SIGNS NEG 

RET 

NZ 


LD 

A,(IX) 


CP 

(IY) 


RET 




#H TABLE 


f 

! TABLE OF POWERS OF 10 

* 

r 

LENTAB EQU 11 
LOTAB EQU $-5 

DEFW 0020H J FIRST 7 VALUES ARE INTEGERS 

DEFW 8400H 

DEFB 1 

DEFW 0048H 

DEFW 870OH 

DEFB 2 

DEFW 007AH 

DEFW 8AOOH 

DEFB 3 

DEFW 401CH 

DEFW 8E00H 

DEFB 4 

DEFW 5043H 

DEFW 910 0 H 






i 

€ 

t 


f 


« 


DEFB 

5 

DEFW 

2 47 4 H 

DEFW 

9400H 

DEFB 

6 

DEFW 

9 618 H 

DEFW 

9880 H 

DEFB 

7 

HI TAB EQIJ 

f-5 

DEFW 

0BC3EH 

DEFW 

9B20H 

DEFB 

8 

DEFW 

1B0EH 

DEFW 

OBACAH 

DEFB 

1A 

DEFW 

0C253H 

DEFW 

0001CH 

DEFB 

24 

TOPTAB DEFW 

0C51DH 

DEFW 

OEBAEH 

DEFB 

32 

7 

; REFERENCES 



THIS VALUE EXACT 

THESE 3 VALUES ROUNDED 
(1.000000027E1A) 

(1« QD0000014E24) 

(1.000000033E32) 


I 

I 

< 

I 

I 

i 

i 

i 


i 

€ 

i 


EXTERNAL MULT,DIVI»ERROR,OVER,UNDER>SEED 
EXTERNAL ZERO , ONE , FI..032 , FLO321, FI..D , FIX3 2 
E XTERNAL. RC2 - L.HARK , LREST , LBACK ,LIN 
GLOBAL DTOB,BTOD 

GLOBAL RND,RANDMZ,INT,ABS,SGN,NEG 

GLOBAL IADD,ISUB,IMUL,IOIV,TABS ? ISGN,ICR,INEG 

*H POWER 


I 

i 

I 


< 

i 

« 

« 


! X TO THE Y POWER 

r 

; INPUT; IX PTS TO X 

? IY PTS TO Y 

! OUTPUT: IX PTS TO X+Y 

? ERRORS; NONE DIRECTLY 

; SPECIAL CASES; IF Y IS 0, RETURN i 
J ALGORITHM; FOR INTEGER Yr A SEQUENCE OF 
J SQUARING AND MULTIPLYING BY X 

; PRODUCES THE RESULT 

J FOR NONINTEGER Y, RESULT IS EXP(Y*LOG(X)) 


I 

4 

i 

4 


€ 


DETERMINE IF Y IS INTEGER 




4 


f ROW PUSH IX ; SAVE IX 4 

PUSH IY ? COPY IY TO IX 

POP IX 

f CALL FIX32 ; GET INTEGER PART 4 

POP IX J RESTORE IX 

OR A 

gj JR Z,ROW10 ; Y WAS AN INTEGER | 

7 

; FOR Y A NONINTEGERr TAKE EXP(Y*LOG<X)> 

t ? I 

PUSH IY ; SAVE IY 

CALL LOG 5 TAKE LOG OF X 

t * 


i 


4 



t 


POP 

IY 


% 



CALL 

MULT 

: Y * L 0 G ( X ) 




CALL 

EXP 



t 

T 

RET 



1 


T 

; FOR Y 

AN INTEGER, MULTIPLY 


x 

r 




€ 


POWIO 

XOR 

A 

? CHECK FOR ZERO 




CP 

(IY+3) 



t 


JP 

Z, ONE 

; X TO 0 POWER 

« 



PUSH 

HL 

; SAVE HL 




EXX 


J ALT BANK 

i 

€ 


PUSH 

HL 

; SAVE H'L' 



BIT 

7, B 





CALL 

NZ , REC 

; NEGATIVE POWER 


C 


CALL 

ST 01 

} SAVE COPY OF X IN TEMPI 

i 



POP 

HL 

; RESTORE INTEGER 




EXX 


; MAIN BANK 


< 


POP 

HL 


« 



LD 

B, 32 

; LOOP COUNTER 


€ 

r 

; LOOK 

FOR FIRST 1 BIT 


i 


T 

P0W20 

PUSH 

BC 

? SAVE COUNTER 

X 

€ 


ADD 

HL, HL 

SHIFT LEFT 1 



EXX 


? ALT BANK 




ADC 

HL, HL 



t 


EXX 


; MAIN BANK 

i 



,JR 

C,POW40 

; 1 BIT FOUND 


X 


POP 

BC 

? GET COUNT 



DJNZ 

P0W20 

; LOOP NOT DONE 

< 


* 

JP 

ONE 

? THIS JUMP SHOULD NOT HAPPEN 


f 

r 

; SQUARE FOR 

EACH ZERO 

BIT 

< 


! SQUARE AND 

MULTIPLY 

BY X FOR EACH I 


1 

1 

P0W30 

PUSH 

BC 

} SAVE COUNT 

i 



ADD 

HL , HL 

1 SHIFT LEFT i 




PUSH 

HL 

? SAVE HL 


X 


EXX 


; ALT BANK 

« 



ADC 

HL , HL. 





PUSH 

HL 

; SAVE H’L’ 


i 


PUSH 

AF 

; SAVE CY 

I 



PUSH 

IX 

; COPY IX TO IY 




POP 

IY 



€ 


CALL 

MULT 

! SQUARE 

XI 



POP 

AF 

? GET CY BIT 




JR 

NC , P0W35 

; ZERO BIT WAS SHIFTED 


€ 


LD 

IY,TEMPI 


« 



CALL 

MULT 

; MULTIPLY BY X 



P0W35 

POP 

HL 

; RESTORE H’L’ 


i 


EXX 


MAIN BANK 

X 



POP 

HL 

; RESTORE HL 



P0W40 

POP 

BC 

; GET COUNT 


€ 


DJNZ 

P0W30 

} LOOP NOT DONE 

X 



RET 




« 

t 

; INTEGER VERSION, FLOATS ARCS FIRST, THEN 

X 


CAL 

LS POW, FIXING 

RESULT, USES TEMP2 , 



; TEMP3 SINCE ROW WILL FIND Y AN INTEGER 


« 





X 

f 





< 


I ROW 


PUSH IX ! SAVE ORIGINAL PTR 
LD Lr(IY) ? GET Y 

LD H , (IY+.1 ) 

LD IX ,TEMP3 

CALL. FLOAT! } FLOAT Y VALUE IN TEMP3 

POP IX ; GET ORIGINAL PTR 


PUSH 

IX 

LD 

L , (IX ) 

LD 

H, ( IX + 1) 

LD 

IX,TEMP2 

CALL 

FLOAT1 

LD 

IY,TEMP3 

CALL 

POW 

CALL 

FIX 

POP 

IX 

L.D 

(IX),L 

LD 

(IX+1),H 

RET 



* H SQUARE ROOT 


! GET X 


FLOAT X VALUE IN TEMP2 
PT TO Y 

USE REAL ROUTINE 
FIX RESULT 
GET ORIGINAL IX 
STORE RESULT 


1 


TTttTTTTrTTTTTrTtTtTTTTrtTlIfTf 


* 


; SQUARE ROOT FUNCTION 

; INPUT; IX PTS TO X 

5 OUTPUT; IX PTS TO SQR(X) 

; ERRORS; X IS NEGATIVE, RETURNS SQR(ABS(X)) 

J SPECIAL CASES; X IS 0 

; ALGORITHM; INPUT IS SCALED TO .25 <= X < 1 
; BY REMOVING AN EVEN EXPONENT. HALF THE EVEN 

? EXPONENT IS ADDED TO THE RESULT EXPONENT 

; AN INITIAL APPROXIMATION IS GIVEN BY 

! A 1ST ORDER POLYNOMIAL, AND IMPROVED TO 

; REQUIRED ACCURACY IN 3 ITERATIONS OF 

! NEWTON’S METHOD. 


; CHECK VALID ARGUMENTS 


LD 

A,(IX + 3) 

OR 

A 

JP 

Z , ZERO 

BIT 

7,(IX) 

CALL 

NZ,SQERR 


; GET EXP 
} ZERO INPUT 
; NEGATIVE ARG 


1 

! EXP OF RESULT IS ROUND(ORIGINAL EXP/2) 
; GIVE X EITHER 0 OR -! EXPONENT 


SUB 

80H 

9 

REMOVE BIAS 

LD 

B, A 

t 

SAVE UNBIASED EXP 

SRA 

A 

r 

DIVIDE BY 2 

ADC 

A , 0 

* 

r 

ROUND 

PUSH 

AF 

t 

SAVE EXP OF RESULT 

ADD 

A, A 

* 

r 

MULT BY 2 

SUB 

B 

t 

SUB ORIGINAL EXP 

NEG 

ADD 

A , 8 0 H 

y 

RESTORE BIAS 


CALC INITIAL APPROX 





DEC 

A 

T 

DIVIDE BY 2 

LD 

(IX + 3) , A 



CALL 

ST 01 

T 

SAVE X/2 IN TEMPI 

LD 

IY,PSQOO 



CALL 

ADD 



r 

i PERFORM 3 

ITERATIONS 



r 

LD 

B, 3 



S 010 PUSH 

BC 

t 

SAVE COUNTER 

CALL 

ST 02 

t 

SAVE LAST APPROX IN TEMP2 

LD 

HL, TEMPI 



CALL 

LD 

r 

GET X/2 

L.D 

IY,TEMP2 



CALL 

DIVI 

7 

DIVIDE BY LAST APPROX 

DEC 

(IY + 3) 

T 

DIVIDE LAST APPROX BY 2 

CALL 

ADD 

r 

NEW APPROX 

POP 

BC 

r 

GET COUNTER 

DJNZ 

♦ 

801 0 

7 

LOOP NOT DONE 

/ 

; STORE EXPONENT 



T 

POP 

AF 


GET SAVED EXP 

ADD 

A , 80H 

* 

r 

ADD BIAS 

L.D 

(IX+3) , A 



RET 




r 

; SQUARE ROOT OF NEGATIVE 

NUMBER 

T 

SQERR LD 

HL , ERROR 



SET 

3 , (HL) 

r 

FLAG ERROR 

RES 

7 , ( IX ) 

7 

TAKE ABS VALUE 

RET 




*H EXP 





TTtrtyttrtt 

t r 

77777777777777777777777777 


; EXPONENTIATION ROUTINE 

T 

; INPUT: IX PTS TO X 

f OUTPUT: IX PTS TO E+X 

ERRORS: OVER/UNDERFLOW IN HANDLING EXPONENT PART 
? SPECIAL CASES: INTERNALLY, FRACTION PART IS ZERO 
! ALGORITHM: MULTIPLY BY LOG BASE 2 OF E AT BEGINNING 
THEN CALCULATE 2 + X. INTEGER PART OF X IS ADDED 
! TO EXPONENT OF RESULT. FRACTION PART IS SCALED 

: TO 0 <= X < ,25 BY REMOVING HIGH ORDER 2 BITS 

I AND SCALING BY 1, 2+.75, 2*. 5, OR Z+.25 AT END. 

f APPROXIMATION TO 2 + X IS GIVEN BY: 

: <0<X+2) + X*P<X+2>> / (0<X+2) - X*P<X+2>> 

J WHERE 0 IS -1ST ORDER POLYNOMIAL, P IS OTH ORDER. 

? REFERENCE IS "COMPUTER APPROXIMATIONS", HART, JOHN 

; ET AL, FUNCTION NUMBER 1080 


; CONVERT E+X TO 2 + X 

r 

EXP LD IY, I..GE4 5 MULT BY 4*L0G2(E> 

CALL MULT 



; GET 8 BIT 2’S COMPl. INTEGER PART IN A 
} GET HIGH 2 BITS OF FRACTION PART IN OE 


CALL 

FIXA 

7 

GET 16 BIT INTEGER PART 

JP 

CfEXERR1 

r 

INTEGER OVERFLOW 

PUSH 

HL 

t 

SAVE INTEGER 

CALL 

BRK 

T 

CALC FRACTION PART 

P OP 

HL 

7 

RESTORE INTEGER 

SR A 

r< r;» 

H 

L 

y 

DIVIDE INTEGER BY 4 

ft ft. 

RRA 

y 

SAVE 2 BITS IN A 

SRA 

H 



RR 

L. 



RLA 




PL. A 


y 

2 BITS NOW IN A 

RLA 


* 

r 

MULT BY 4 

RLA 




AND 

OCH 



L D 

E , A 



LD 

DrO 

* 

r 

DF. IS NOW TABLE DISPLACEMENT 

LD 

A,H 

r 

CHECK SIZE OF INTEGERR 

RLC 

L. 

7 

SIGN BIT IN CY 

ADC 

A , 0 



JR 

NZ » EXERR2 

* 

r 

INTEGER WAS > 8 BIT 2’S COMP 

LD 

RRCA 

A , L 

7 

RESTORE 8 BIT INTEGER 

PUSH 

AF 

* 

t 

SAVE INTEGER PART 

PUSH 

DF: 

t 

7 

SAVE NEXT TWO BITS 

XOR 

A 



CP 

( IX + 3) 



JR 

Z i EXIO 

7 

FRACTION PART IS ZERO 

DEC 

(IX + 3 ) 

t 

r 

SCALE TO 0 <= X < 1/4 

DEC 

CIX+3) 

7 

UNDERFLOW CANNOT OCCUR HERE 


7 

! EVALUATE RATIONAL APPROX. TO 2+X 


CALL 

ST Oi 

7 

SAVE 

A COPY 

OF X 

IN TEMPI 

LD 

IY ,PEXQO 






CALL 

MULT 






CALL 

ST 02 

7 

SAVE 

X*P<X+ 

2) IN 

TEMP2 

LD 

HL ,TEMPI 






CALL 

LD 

7 

GET ANOTHER 

COPY 

OF X 

LD 

IY * TEMPI 






CALL 

MULT 

> 

CALC 

X + 2 



LD 

IYrQEXOO 






CALL 

ADD 






CALL 

ST 01 

J 

SAVE 

Q( X+2 ) 

IN TEMPI 

LD 

IY j TEMP2 






CALL. 

SUB 






CALL 

ST 03 

7 

SAVE 

Q-X#P 

IN TEMP3 

LD 

HL,TEMPI 






CALL 

LD 

7 

GET G 

j 



LD 

IY, TEMP2 






CALL 

ADD 

7 

CALC 

Q + X*P 



LD 

IY,TEMP3 






CALL, 

DIVI 

t 

T 

FINAL 

. RESULT 



; ADD IN FRACTIONAL AND INTEGER POWERS OF 2 




; GET TABLE DI8PL 


POP DE 
LD IY,EXTAB 
ADD IY ? DE J POINT TO FRACTIONAL POWER 

LD A , E 
0 R A 

CALL NZ 7 MULT ; MULTIPLY BY TAB ENTRY 
} ADD INTEGER PART TO EXPONENT OF RESULT 



POP 

BC 

7 

INTEGER PART IN B 


LD 

A,(IX + 3) 

r 

GET EXP 


SUB 

80H 

7 

REMOVE BIAS 


ADD 

A, B 




JP 

F‘E , EXERR3 

7 

OVER/UNDERFLOW 


ADD 

A, 80H 

7 

RESTORE BIAS 


JP 

Z , UNDER 

7 

UNDERFLOW 


LD 

CIX+3),A 




RET 




r 

; EXP 

ERRORS 




T 

EXERRi 

LD 

H , (IX > 

7 

GET ORIGINAL SIGN 

EXERR2 

BIT 

7 , H 

1 

r 

TEST SIGN 


JP 

NZ, UNDER 

7 

NEGATIVE SIGN 


JP 

OVER 



EXERR3 

JP 

P , UNDER 

» 

7 

POS SIGN AFTER OVERFLOW 


JP 

OVER 



*H LOGARITHM 

t *********** 


* » 

************************ 


; LOG ROUTINE 

7 

; INPUT; IX PTS TO X 

; OUTPUT; IX PTS TO NATURAL LOG OF X 

; ERRORS; X IS NEGATIVE, RETURN L.OG<ABS<X>> 

? X IS ZERO, RETURN NEGATIVE OVERFLOW 

; SPECIAL CASES; X IS 1, RETURN ZERO 
? ALGORITHM; EXPONENT IS STIPPED OFF, FLOATED, AND 
; ADDED TO RESULT AT END. REMAINING FRACTION 

; PART IS MULTIPLIED BY SQR<2> TO GIVE RANGE 

; OF SQR (2.) / 2 <= X < SQR < X) . LOG FUNCTION 

J IS APPROXIMATED BY; 

; z * p ( z2) 

; WHERE Z = (X-l) / (X+l) AND P IS 2ND ORDER. 

; REFERENCE IS FUNCTION 2661. 

LOG (SQR { 2 ) ) IS SUBTRACTED FROM RESULT TO CORRECT 
; FOR SCALING. 


7 7 7 7 7 7 7 

7 7 7 7 7 

77777777777777 

******** 

Ttrirrrr 

r 

; CHECK 

FOR 

VALID ARGUMENT 

8 

7 

L OG 

BIT 

7, (IX) 



CALL 

NZ,LGERR1 ; 

NEGATIVE 


LD 

IY,UNITY 



CALL 

CP 



JP 

Z,ZERO ; 

LOG ( 1 ) 


LD 

A,(IX + 3 5 



OR 

A 



JP 

Z,LGERR2 ; 

ZERO ARG 



MAKE EXP INTO 16 BIT 


COMPLEMENT 


SUB 

80H 

; REMOVE BIAS 

l.D 

L , A 

? PUT IN HL 

LD 

H i 0 


JR 

NC r L.010 

J SIGN IS POSITIVE 

LD 

H r OFFH 

; MAKE NEGATIVE 

r 

; CONVERT EXP 

t 

■ TO LOG BAS 

E E 7 SAVE IN TEMP3 

T 

LOIO PUSH 

IX 

; SAVE ACC PTR 

LD 

IX r TEMP3 


CALL 

FLOAT! 

; FLOAT HL 

LD 

IY r LN2 


CALL 

MULT 

; MULT BY LOG(2) 

POP 

IX 


r 

; SCALE TO 1/ 

S Q R(2 ) <= X 

< S0R(2) 

r 

LD 

(IX+3)» 80H 

; 1/2 <= X < 1 

LD 

IY f SQR2 


CALL 

MULT 

; MULT BY SQR(2 > 

? 

! APPROXIMATION FOR LOGARITHM 

r 

CALL 

ST 01 

; SAVE X IN TEMPI 

PUSH 

IX 

; SAVE PTR TO ACCIJM 

LD 

IX 7 TEMPI 


LD 

IY,UNITY 


CALL 

ADD 

? X+l IN TEMPI 

POP 

IX 


CALL 

SUB 

? X-l IN ACCIJM 

LD 

IY r TEMPI 


CALL 

DIVI 


CALL 

ST 01 

; (X-l)/(X+l) IN TEMPI 

CALL 

MULT 


CALL 

ST 02 

) Z+2 IN TEMP2 

LD 

IY t PL002 

? CALC P(Z+2) 

CALL 

MULT 


LD 

IY t PLOD 1 


CALL 

ADD 


LD 

IY p TEMF‘2 


L A L L 

MULT 


LD 

IY 7 PLOOO 


CALL 

ADD 


LD 

I.Y 7 TEMPI 

; CALC 2 * P(Z*2) 

CALL 

MULT 


LD 

IY 7 LNSQ2 


CALL 

SUB 

J SUB LOG(SQR(2)) 

LD 

IY 7 TEMP3 


CALL 

ADD 

? ADD EXP PART 

RET 

V 



T 

ERRORS 



r 

LGERRl LD 

HL 7 ERROR 


SET 

2, (HL.) 

; FLAG LOG OF NEG 

RES 

7 7 (HL) 

} TAKE ABS VALUE 

RET 






LGERR2 SET 7 , < IX > } SET NEGATIVE SIGN 

JP OVER ; LOG OF ZERO 

*H TANGENT 


! TANGENT ROUTINE 

t 

; INPUT: IX PTS TO X (IN RADIANS) 

J OUTPUT: IX PTS TO TAN(X) 

! ERRORS: NONE DIRECTLY 

J ALGORITHM: RANGE REDUCTION TO 0 < = X < = PI/4. 

J MULTIPLY BY 4/PI. LOW 3 BITS OF INTEGER 

? PART GIVE THE OCTANT NUMBER, AND IS USED 

5 AS A FLAG. TAKE FRACTIONAL PART, IGNORE 

: INTEGER PART (EXCEPT FLAG BITS). 

? OCTANTS 1,3,5,7 FORM 1 FRACTION PART. 

? APPROXIMATION IS GIVEN BY: 

! X*P(X+2) / Q(X+2) 

; WHERE P IS 1ST ORDER AND 0. IS 2ND ORDER POLY. 

J REFERENCE FUNCTION IS 4282. 

: OCTANTS 1,2,5,6 FORM RECIPROCAL OF RESULT. 

i OCTANTS 2,3,6,7 RETURN NEGATIVE SIGN. 


r 

; SCALE BY 4/PI 


TAN 

LD 

IY , F 0 V R PI 




CALL 

MULT 




CALL 

BRK 

r 

BREAK INTEGER AND FRACTION 


LD 

A,L 




PUSH 

AF 

? 

SAVE FLAG BITS 


AND 

1 

F 

TEST FLAG 


JR 

Z , T A10 

y 

THEY WERE 0 OR 2 


LD 

IY,UNITY 




CALL 

RSUB 

y 

FORM 1 - FRACTION PART 

F 

; TAN 

APPROX 




y 

T AID 

CALL 

ST 03 

y 

SAVE X IN TEMP3 


LD 

IY , T E M P 3 




CALL 

MULT 




C ALL 

ST 02 

T 

y 

SAVE X+2 IN TEMP2 


LD 

IY,QTA01 




CALL 

ADD 




L.D 

IY, TEMF‘2 




CALL 

MULT 

y 

MULT BY X+2 


LD 

IY , QTAOO 




CALL 

ADD 




CALL 

ST 01 

t 

y 

SAVE E)(X + 2) IN TEMPI 


LD 

HL,TEMP2 




CALL 

LD 

y 

GET X+2 


LD 

IY,PTA01 




CALL 

MULT 




LD 

IY, F'TAOD 




CALL 

ADD 




LD 

IY , TEMP3 




CALL 

MULT 

y 

X * P(X+2) IN ACCUM 


POP 

AF 

r 

GET 2 BIT FLAG 


PUSH 

AF 






I 


i 

€ 



AND 

3 


JP 

PE , TA20 


CALL 

ST02 


LD 

HL , TEMPI 


CALL 

L.D 


LD 

IY,TEMP2 


JP 

TA30 

TA20 

LD 

IY,TEMPI 

TA30 

CALL 

DIVI 


; GET SIGN OF RESULT 


; THEY WERE 0 OR 3 
? SAVE X*P IN TEMP2 


GET 

Q 


SET 

FOR 

Q/XP 

SET 

FOR 

XP/G 


r 


m 

POP 

AF 


AND 

O 

rf- 

« 

RET 

*7 

SET 

7 , ( IX ) 


RET 



*H COSINE, S 

INE 

4 " 

t **** ******* 

irryTTfirirr 

******* 
r r r r y y r 


GET 2 BIT FLAG 

? ON 0 OR 1 
? MAKE RESULT NEG 


y y y y 


* * * 
y y y 


y y y y y y 


yyyyyyyyyyyyy 


4 

4 

4 

4 

4 

4 


t 


COSINE ROUTINE 

EFFECTIVELY ADDS PI/2 AND CALLS SINE 


4 


€ 

i 


y y y 

COS 


? y y 


yyyyyyyyyyyyyyyyy 

LD lYrFOVRPI 
CALL MULT 
CALL BRK 
L. D A > L 
A DD A , 2 
JP SHQ 


yyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyy 


; SCALE BY 4/PI 

; SEPARATE INTEGER FROM FRACTION 
; ADD EFFECTIVE PI/2 


yyyyyyyy 


i * * » 

y y y y 


y y yyyyyyyyyyyyyyyyyyyyyyyy 


yyyyyyyyyyyyy 


yyyyyyyyyyyyy 


4 

4 

4 


4 

4 

C 

4 

C 

€ 

€ 

« 

4 


; SINE ROUTINE 

* 

r 

INPUT: IX PTS TO X (IN RADIANS) 

f OUTPUT: IX PTS TO SIN(X) 

ERRORS: NONE 

; ALGORITHM: RANGE REDUCTION TO D <= X <= PI/4 
: BY MULTIPLYING BY 4/PI AND TAKING FRACTION 

PART, LOW 3 BITS OF INTEGER PART GIVE 
; OCTANT NUMBER. 

; OCTANTS 1,3,5,7 FORM 1 - FRACTION PART 

i OCTANTS 0,3,4,7 USE SINE APPROXIMATION: 

! X*P(X+2> / Q(X*2) 

i WHERE P IS 2ND ORDER AND Q IS 1ST ORDER POLY. 

; REFERENCE FUNCTION IS 3060. 

: OCTANTS 1,2,5,6 USE COSINE APPROXIMATION: 

! P (X + 2 ) / Q ( X t ' 2 ) 

; WHERE P IS 2ND ORDER AND G IS 1ST ORDER POLY. 

! REFERENCE FUNCTION IS 3840. 

; OCTANTS 4,5,6,7 RETURN NEGATIVE SIGN. 


SIN 

LD 

IY,FOVRPI 



CALL 

MULT 

5 SCALE BY 4/PI 


CALL 

BRK 

? GET INTEGER AND FRACTION 


LD 

A r L 


SI 10 

PUSH 

AF 

! SAVE 3 BIT FLAG 


AMD 

1 



JR 

Z , S12 0 

: 0,2,4, OR 6 


yyyyyyyyy 


4 

4 

4 

4 

4 

4 

4 

4 

4 

4 




4 




i 

S X 2 0 

t 

« ; 

; SIN 

t 

€ 

€ 

i 

I 

€ 

« 

f 

* 

r 

; COS 

t ; 

S 130 

< 

f 

i 

€ 

t 

C 

S 140 

1 

i 

*H ARC 

t ?? 7 t T 

« 


LD 

IY , UNITY 



CALL 

RSIJB 

l 

CALC 1 -• FRACTION PART 

POP 

AF 

t 

GET 3 BIT FLAG 

PUSH 

AF 



AND 

3 



JP 

PO,8130 

1 

1 7 2 7 5 OR 6 

APPROX 




CALL 

ST 03 

7 

SAVE X IN TEMP3 

LD 

IY r TEMP3 



CALL 

MULT 

7 

SQUARE X 

CALL 

8 T 0 2 

T 

SAVE Xt-2 IN TEMP2 

LD 

IY.QSIOO 



CALL 

ADD 



CALL 

ST 01 

r 

SAVE Q(X+2) IN TEMPI 

LD 

HL , TEMP2 



CALL 

LD 

t 

GET X + 2 

LD 

IY/PSI02 



CALL 

MULT 



LD 

IY 7 F‘ S10 1 



CALL 

ADD 



LD 

IY 7 TEMP2 



CALL 

MULT 



LD 

IY 7 P S10 0 



CALL 

ADD 



LD 

IY 7 TEMPI 


PT TO 9. 

CALL 

DIVI 

r 

P/Q 

LD 

IY r TEMP3 



CALL 

MULT 



JP 

SI 40 



APPROX 




PUSH 

IX 

1 

COPY IX TO IY 

POP 

IY 



CALL 

MULT 



CALL 

ST 02 

7 

SAVE X+2 IN TEMP2 

LD 

IY'QCOOO 



CALL 

ADD 



CALL 

ST 01 

7 

SAVE Q(X+2> IN TEMPI 

LD 

HL 7 TEMP2 



CALL 

LD 

1 

GET X+2 

L.D 

IY 7 PCOO 2 



CALL 

MULT 



LD 

IY,PCOO 1 



CALL 

ADD 



LD 

IY 7 TEMP2 



CALL 

MULT 



LD 

IY 7 PCOOO 



CALL 

ADD 



LD 

IY 7 TEMP 1 



CALL 

DIVI 



POP 

AF 

1 

GET 3 BIT FLAG 

AMD 

4 



RET 

Z 

1 

4 , 5,5 OR 7 

SET 

7, ( IX) 

r 

NEGATE RESULT 

RET 




TANGENT 



9 * V V ? * 

f t ? 17 t 

fit t r l l y 1 l 

fit 

tTtttrtfftfftrttftttfft 


« 

i 

i 

€ 

i 

< 

< 

i 

< 

< 

« 

i 

i 

« 

i 

« 

« 

« 

« 

( 

« 


I 


f 



; ARCTANGENT ROUTINE 


INPUT; IX PTS TO X 

OUTPUT; IX PTS TO ATN(X) (IN RADIANS) 

ERRORS; NONE 

ALGORITHM; RESULT WILL HAVE SIGN OF INPUT. 
RANGE REDUCTION BY TAKING ABS VALUE AMD 
TESTING FOR 1 OF 3 REGIONS; 

1) X < T A N < PI/8) USE APPROX. DIRECTLY 

2) X < TAN(3*PI /8) REPLACE X BY 
<X~1> / (X+l) AND ADD PI/4 AT END 

3) REPLACE X BY i/X AND SUB FROM F'I/2 
AT END. 

APPROXIMATION FOR 0 <= X <= TAN(PI/8); 

X * P ( X + 2 ) / Q(X+2) 

WHERE. P IS 1ST ORDER AND G 2ND ORDER POLY. 
REFERENCE FUNCTION IS 5071. 


; TAKE ABS VALUE, DETERMINE RANGE 


ATN L..D A, (IX) 

AND 80H 
PUSH AF 
RES 7,(IX) 

LD IY f TAPIS 
L- A L L U P 
LD B,0 
JR C , A T 2 0 

LD IY r TA3PI8 
CALL CP 
JR CfATIO 

F 

; TAN(3*PI /8) <= X 

; TAKE 1/X NOW 


GET SIGN 

SAVE SIGN 
MAKE X POSITIVE 

COMPARE TO 1ST LIMIT 
SET FLAG TO 0 

< TAN(PI/8) 

COMPARE TO 2ND LIMIT 

< TAN(3*PI/8) 


IJB FROM PI/2 AT END 


CALL REC i TAKE RECIPROCAL 

LD B r 2 i SET FLAG TO 2 

JP AT20 

* 

r 

; TAN(PI/8) <= X < TAN(3*PI/8) 

? CALC (X-1)/(X + 1) NOW , 
i ADD PI/4 AT END 


AT 10 

CALL 

ST 01 


PUSH 

IX 


LD 

IXfTEMPI 


LD 

IYfUNITY 


CALI- 

ADD 


POP 

IX 


CALL 

SUB 


L.D 

IY,TEMPI 


CALL 

DIVI 


LD 

B, 1 


; SAVE X IN TEMPI 
J SAVE PTR TO ACCOM 

; X+i IN TEMPI 
; X-l IN ACCUM 

J SET FLAG TO 1 


? APPROX FOR ARCTAN 


AT20 PUSH BC 

, SAVE FLAG 





CALL 

ST 03 

LD 

IY , TEMF'3 

CALL 

MULT 

CALL 

ST 02 

LD 

IY , Q A T 01 

CALL 

ADD 

LD 

IY , TEMF'2 

CALL 

MULT 

LD 

IY ,QATOO 

CALL 

ADD 

CALL 

ST 01 

LD 

HL ,TEMP2 

CALL 

LD 

LD 

IY ,PAT01 

CALL 

MULT 

L.D 

IY,PATOO 

CALL 

ADD 

LD 

IY,TEMPI 

CALL. 

DIVI 

LD 

IY, TEMF'3 

CALL. 

MULT 

P OP 

DC 

BIT 

Q,B 

JR 

Z, AT30 


} SECOND INTERVAL 

L.D IY , P14 
CALL ADD } 

JP AT 40 

AT30 BIT 1,B 

JR Z,AT40 ? 

; THIRD INTERVAL 

LD IY > P12 
CALL RSUB J 

T 

; GIVE RESULT SIGN OF ARG 


AT 40 

POP 

AF 

t 


OR 

(IX) 



L D 

«IX ) , A 



RET 



* H CON 

ST ANTS 



; cons 

TANTS 



UNITY 

DEFW 

0 

r 


DEFW 

810 0 H 


LGE4 

DEFW 

0AA38H 

» 


DEFW 

833BH 

r 

PEXOO 

DEFW 

950AH 

? 


DEFW 

84DBH 

r 

QEXOO 

DEFW 

0EF47H 

* 

7 


DEFW 

85BFH 

T 

EXTAB 

EQU 

$-4 

7 

7 


DEFW 

3718H 

7 


DEFW 

81F0H 

r 

7 


DEFW 

0435H 

7 


DEFW 

81F3H 

t 

7 


DEFW 

4 4 5 7 H 

7 


DEFW 

81FDH 

r 


SAVE K IN TEMP3 


SAVE X+2 IN TEMP? 


SAVE Q(X+2) IN TEMPI 
GET X+2 


GET FLAG 

1ST OR 3RD INTERVAL 


ADD PI/4 


1ST INTERVAL 


SUB FROM PI/2 


GET ORIGINAL SIGN 


FL PT 1.0 

4 * LOG BASE 2 OF E 
5.770780164 
P COEFFICIENTS 
8.6615859416262 
Q COEFFICIENTS 
24.9920639458976 
QEX01 IS 1.0 
FRACTIONAL POWERS OF 
2 + 1/4 

1.189207115 
2 + 1/2 

1.414213562 
2 + 3/4 

1.681792831 


PSQOO 

DEFW 

0F655H 


DEFW 

7F1AH 

LN2 

DEFW 

7231H 


DEFW 

8018H 

LNSQ2 

DEFW 

7 2 31H 


DEFW 

7F18H 

SQR2 

EQU 

EXTAB+8 

PL 000 

DEFW 

OOOOH 


DEFW 

8204H 

PL 001 

DEFW 

9B2AH 


DEFW 

80D1H 

PL002 

DEFW 

925 4 H 


DEFW 

7F22H 

FOVRPI 

DEFW 

0F922H 


DEFW 

8183H 

PTAOO 

DEFW 

6C54H 


DEFW 

88A9H 

PTA01 

DEFW 

0DAC8H 


DEFW 

844EH 

QTAOO 

DEFW 

3B07H 


DEFW 

89CEH 

QTA01 

DEFW 

318FH 


DEFW 

872FH 

J QTA02 

IS 1 

.0 

PSIOO 

DEFW 

4653H 


DEFW 

863FH 

P S101 

DEFW 

0A294H 


DEFW 

8335H 

PS 102 

DEFW 

0AC31H 


DEFW 

7D59H 

QSIOO 

DEFW 

80Q6H 


DEFW 

8760H 

; qsioi 

IS 1 

.0 

PC 000 

DEFW 

0BF3EH 


DEFW 

86 CAH 

PC001 

DEFW 

53DBH 


DEFW 

84F8H 

PC002 

DEFW 

4865H 


DEFW 

7FF8H 

QCOOO 

DEFW 

0BF3EH 


DEFW 

86C9H 

! Q C 001 

IS 1 

. 0 

PATOO 

DEFW 

8F4AH 


DEFW 

844EH 

PAT01 

DEFW 

0D04BH 


DEFW 

8365H 

QATOO 

DEFW 

8F4AH 


DEFW 

844EH 

QAT01 

DEFW 

6D29H 


DEFW 

8 4 0 0 H 

? QAT02 

IS 1 

.0 

TAPIS 

DEFW 

1354H 


DEFW 

7FCDH 

TA3PI8 

DEFW 

821 AH 


DEFW 

827AH 

P 12 

DEFW 

0F49H 


DEFW 

81DBH 

P 14 

DEFW 

0F49H 


DEFW 

80DBH 


J SQR INIT APPROX 
; .4178932188 

J LOG BASE E OF 2 
; .6931471806 

; LOG BASE E OF SQR( 
; .3465735903 

; SQUARE ROOT OF 2 
? P COEFFICIENTS 
? 2.0000008368 

; .6664400777 

.41517739 
; 4/PI 

; i.273239545 
; P COEFFICIENTS 
? 212.42445758263 

5 -12.55329742424 
J Q COEFFICIENTS 
i 270.46722349399 

; -71.59606050466 

; P COEFFICIENTS 
J 52.81860134812 

? -4.644800481954 

? .0867545069521 

, Q COEFFICIENTS 
J 67.250731777791 

? P COEFFICIENTS 
J 47.68729218663 

; -13.70800004765 

; .447822343969 
J Q COEFFICIENTS 
i 47.68729082658 

} P COEFFICIENTS 
/ 12.65998609915 

? 6.3691887127 
? Q COEFFICIENTS 
? 12.65998646243 

5 10.5891113168 

} T AN(PI/8) 

; .4142135624 

; TAN(3*PI/8) 

; 2.414213562 
; PI / 2 
? 1,570796327 
; PI / 4 
i .7853981634 





REFERENCES 


GLOBAL. POW , IPOW , 8QR , EXP , LOG , SIN, COS , TAN , ATN 
EXTERNAL OVER,UNDER,ZERO,ERROR,ONE 

EXTFRNAL FIX32,FL032,ADD,SUB,RSUB 

EXTERNAL MULT,DIVI,REC,STOl,ST02 

EXTERNAL ST03,TEMPI,TEMP2,TEMP3,LD 

EXTERNAL CP,FIX,FIXA,FLOAT,BRK,FLOAT 1 


REQUIRED AS EXTERNALS FOR MATH PACKAGE 


1 ******* 
7 T 7 t T 7 T T 

MM 

M M J) I M M 

! MISCEI. 

.LANEOUS ROUTINE 

* »»»!*»» 
t t r t i y r i 

* * * * 

y t r r 

rrrrrrttttr 

LMARK 

LJ) 

(POINT),HL 


RET 


LREST 

RET 


LBACK 

RET 


RC2 

RET 


LIN 

PUSH 

HL 


L..D 

HL, (POINT) 

LINIO 

LD 

A,(HL) 


INC 

HL 


CP 

r r 


JR 

Z,LINIO 


L.D 

(POINT) , Hl. 


POP 

HL 


RET 


SCP 

RET 


POINT 

DEFS 

2 

ERROR 

DEFB 

0 

ILLNUM 

EQU 

7 

OVRFLW 

EQU 

6 

UNDFLW 

EQU 

5 

DVDZRO 

EQU 

4 

SQRNGN 

EQU 

3 

LOGNGN 

EQU 



REFERENCES 

GLOBAL 


f STORE PTR FOR DTOB 

5 ENTRY FOR UNSUCCESSFUL DTOB (NO DIGITS 
i DUMMY ROUTINE (BACK UP PTR FOR BASIC) 

? NORMAL RETURN, IF SOME KIND OF NUMBER 
} SAVE HL 
J GET PTR 
} GET NEXT CHAR 
5 ADVANCE PTR 

? IGNORE BLANKS 
; SAVE PTR 
i RESTORE HL. 

J STRING COMPARE 
; TEMP PTR STORAGE 
? ERROR FLAGS FOR MATH PACKAGE 
; ILLEGAL NUMBER (DTOB) 

5 OVERFLOW 
? UNDERFLOW 
J DIVISION BY ZERO 
; SQUARE ROOT OF NEGATIVE NUMBER 
} LOG OF NEGATIVE NUMBER 


E R R 0 R,SCP,LMA RK,LREST,LBA CK,R C2,LIN 


rtrtTtty 

FOUND) 

FOUND 


i 

i 

i 

i 

i 

i 

l 

l 

4 

4 

4 

4 

4 

4 


4 

4 

4 


4 


4 

I 

4 


4 






*H ADD AND SUBTRACT 


: FLOATING POINT REVERSE SUBTRACT ROUTINE 
J NEGATE X AND ADD 

t 

777777777777777777777777777777777777777777777777777777777777777777777777 


RSIJB 

LD 

fir( IX ) 

; GET MS 

BYTE 

OF 

X 


XOR 

80H 

? NEGATE 

SIGN 




LD 

D r A 

? SAVE IT 




LD 

Ef <IY> 

; GET MS 

BYTE 

OF 

Y 


JP A005 


; FLOATING POINT SUBTRACT ROUTINE 
} NEGATE Y AND ADD 


SUB 


LD 

A,<IY) 

r 

GET MSBYTE OF 

Y 

XOR 

80H 

•t 

7 

NEGATE SIGN 


LD 

E,A 

t 

SAVE IT 


LD 

D , (IX ) 

7 

GET MSBYTE OF 

X 

JP 

t t % t 

i i y r 

ADOS 

yyyyyyy? 

t t » t F 

7 7 7 7 7 

} M 7 1 } M 1 J 11 II 

1 t 

y r 


* * * If 


FLOATING POINT ADD ROUTINE 


7 

input; 

IX 

PTS 

TO 

7 


IY 

PTS 

TO 

« 

7 

output; 

IX 

PTS 

TO 

T 

ERRORS; 

EXPONENT 

O' 


SPECIAL.CASES; ZERO ARGUMENT(S) 


INPUT FOR ADOS s D CONTAINS MS BYTE OF X 

E CONTAINS MS BYTE OF Y 

ALGORITHM; ABSOLUTE VALUE MAGNITUDES ARE COMPARED, 


f THE 

LARGER IS 

LOADED TO AC 1 7 THE SMAL 

; AC 2 

IS SHIFTED 

TO ALIGN DECIMAL POINT 

T 

i yytTTTTTTiy 

** * ******* 
ytyyrrrtTf 

** *********** ********** 
yyyyyyyyyyyyyytyyyyyyyy 

ADD LD 

D,(IX) 

! GET MS BYTES 

LD 

E , < IY ) 


7 

; CHECK FOR 

ZERO ARGUMENTS 

7 

ADOS XOR 

A 

? CLEAR A 

CP 

( IX + 7 ) 


JP 

N Z 7 A D 0 7 

; X IS NONZERO 

CP 

( IY+7) 


RET 

Z 

J Y IS ALSO ZERO 

LD 

A r E 

; GET Y MSBYTE 

PUSH 

IY 


POP 

HL. 


CALL 

LD 

? COPY Y TO X 

LD 

<IX) , A 

J STORE MSBYTE 

RET 



AD07 LD 

(IX) ,D 

; STRE NEW SIGN IN CA 

CP 

( IY+7) 


RET 

*7 

; V IS ZERO, RETURN X 




OAD ARGUMENTS, LARGER MAGNITUDE IN ACi 


AD.1.0 


PUSH 

DE 

; SAVE MS BYTES 


RES 

7,0 

; CLEAR SIGN BITS 


RES 

7, E 



CALL 

CPI 

; COMPARE X TO Y IGNORING 

JR 

C,AD10 

; X < Y 


CALL 

ST A1 

; X TO ACI 


CALL 

STYA2 

; Y TO AC 2 


POP 

BC 

; GET SIGNS 


LD 

A, B 



X OR 

c 



LD 

A, B 

? USE SIGN OF X FOR 

RESULT 

JP 

AD 15 



CALL 

STYA1 

5 Y TO ACI 


CALL 

STA2 

i X TO AC2 


POP 

BC 

} GET SIGNS 


LD 

A, C 



X 0 R 

B 



LD 

A , C 

; USE SIGN OF Y FOR 

RESULT 

IMAL.IZE BY S 

HIFTING AC2 RIGHT 



A D15 

EX 

AF , AF' 

t 

7 

SAVE RESULT SIGN, XOR FLAGS 


LD 

A,(A C1E) 

y 

GET ACI EXPONENT 


LD 

HL.AC2E 




SUB 

<HL> 

7 

CALC SHIFT COUNT 


JR 

Z, AD30 

7 

FRACTIONS ALREADY ALIGNED 


CP 

13 




JP 

NC»FSTO 

7 

AC2 IS INSIGNIFICANT 


LD 

B, A 

r 

SET COUNTER 

AD 20 

X 0 R 

A 

7 

CLEAR A 


LD 

HL,AC2 

r 

PT TO ACCUM 


CALL 

SRDi 

* 

7 

SHFT RT 1 DIGIT 

* 

DJNZ 

AD20 

t 

r 

COUNTER NOT 0 

r 

; INITIALIZE 

FOR ADD 

OR 8 

UB 

y 

AD 30 

LD 

B, 7 


COUNTER 


LD 

HL,AC20 

y 

POINTERS 


LD 

DE,AGIO 




AND 

A 

r 

CLEAR CY 


EX 

AF , AF' 

? 

GET SIGN & FLAGS 


JP 

M , S B10 

y 

SIGNS WERE DIFFERENT 

; ADD 





f 

EX 

AF , AF' 

7 

SAVE SIGN OF RESULT 

AD40 

LD 

A,(DE) 




ADC 

A, (HL) 




DAA 





LD 

(DE) , A 




DEC- 

DE 




DEC 

HL 



* 

DJNZ 

AD 4 0 



t 

POST 

NORMALIZE, AT 

MOST 

1 SHFT RIGHT 

y 

LD 

A , (ACI) 

> 

CHECK FOR OVERFLOW 


AND 

OF.OH 




JP 

7., FSTO 

? ALREADY NORMALIZED 

LD 

HL,ACIE 


I NC 

(HI...) 

; INCREMENT EXPONENT 

JP 

Z, OVERA 

; OVERFLOW 

LD 

HL., AC1 

5 PT TO ACCUM 

CALL 

SRD1 

? SHIFT RIGHT 

JP 

FSTO 



? SUB 

t 

r 

SB10 

SB20 


EX 

LD 

SBC 

DAA 

LD 

DEC 

DEC 


AF,AF' 
A, (DE) 
A , ( HL. > 

(DE ) 7 A 

DE 

HL. 


SAVE SIGN OF RESULT 


DJNZ SB20 

POST-NORMAL.IZE BY SHIFTING LEFT 



LD 

A, (AC 1E ) 

; PUT EXP IN C 


LD 

C, A 



LD 

HL., AC1 



CALL 

N0RML.2 

5 NORMALIZE FRACTION 


JP 

Z,ZERO 

; ZERO FRACTION RESULT 


JP 

C,UNDER 

; EXPONENT UNDERFLOW 


LD 

A, C 

? STORE EXP 


LD 

( A C1E ) , A 



JP 

FSTO 


t 

» 

I 

COPY AC1 TO (IX) AND 

SET SIGN OF RESULT 

y 

FS 

TO LD 

HL,ACi 

r PT TO ACi 


CALL 

LD 



JP 

RSTO 

? SET SIGN FROM A’ 

*H 

MULTIPLY 



t r 

yyyrryyyyy 

yyyyyyyyyyy 

yyyyyyyyyyyyyyyyyyyyyyyy 

y 

i 

FLOATING POINT MULTIPLY ROUTINE 

7 

7 

INPUT; 

IX PTS TO 

MULTIPLICAND 

r 


IY PTS TO 

MULTIPLIER 

r 

OUTPUT; 

IX PTS TO 

PRODUCT (PRODUCT REPLACE 

r 

ERRORS; 

EXPONENT OVERFLOW OR UNDERFLOW 

? 

SPECIAL 

CASES; ARGUMENT(S) ZERO 

7 

ALGORITHM; X, 2X, 

4X, AND 8X ARE STORED IN 

7 

Y IS 

PLACED IN 

ACI . 


yrTtttrtrttrrjfrtrrr 


MULTIPLIER IS EXAMINED BIT AT A TIME WITHIN EACH DIGIT, 
CORRESPONDING TEMP REGS ARE ADDED TO (IX), AND (IX) 

IS SHIFTED RIGHT. 

TIMING? APPROX 19814 (16916 IN LOOP) 


; CHECK FOR ZERO OPERANDS 

MULT XOR A 

CP (IX+7) 





Jp 

CP 

JP 


Z.ZER'O 
(IY+7) 
Z,ZERO 


CALCULATE 

AND SAVE 

SIGN 

OF RESULT 

LD 

A , < IX ) 

V 

f 

MS BYTE OF X 

XOR 

( IY) 



EX 

AF,AF' 

t 

SAVE IN A' 

STORE X, 2 

X, 4Xr AND 8X 

IN TEMP REGISTERS 

PUSH 

IX 



POP 

DE 



LD 

HL , 6 



ADD 

HL.DE 



LD 

(SAVE ) , 

HL ; 

SAVE PTR TO LS BYTE 

EX 

DE 7 HL 



CALL 

MDTL 



SAVE Y IN 

AC1 (IN 

C A 8 E 

X AND Y ARE SAME) 

CLEAR (IX) 




CALL 

STYA1 

? 

STORE Y IN AC1 

LD 

IY 7 A C i 

r 

PT TO Y 

CALL 

ZERO 

r 

CLEAR X 

LD 

B r 7 

T 

LOOP COUNTER 

JP 

ML 20 

7 

SKIP SHIFT OF AC1 


J SHIFT PRODUCT RIGHT ONE DIGIT 
; GET NEXT MULTIPLIER DIGIT 


7 


M L 10 

XOR 

A 


CALL 

3RD 


CALL 

ML.DIG 


XOR 

A 


CALL 

SRD 

ML20 

LD 

C,(IY+6) 


DEC 

IY 


CALL 

MLDIG 


DJNZ 

ML 10 

; CHECK 

NORMALIZATION 

T 

DEC 

B 


L.D 

A,(IX) 


AND 

OF OH 


JR 

Z t ML60 


INC 

B 


XOR 

A 


CALL 

SRD 

7 

; CALCULATE 

EXPONENT OF 

r 

ML 60 

LD 

ft, (TEMPOE) 


SUB 

8 OH 


LD 

C, A 


L.D 

A,(IY+14) 


SUB 

80H 


ADD 

A, C 


7 CLEAR A 
i SHFT (IX) RIGHT 
; MULTIPLY BY 2ND DIGIT 

? SHFT (IX) RIGHT 
i GET 2 DIGITS 
; DECR PTR 

i MULTIPLY BY 1ST DIGIT 
? COUNTER NOT 0 


i DEFAULT SHIFT COUNT I 
} GET MS BYTE 

f ALREADY NORMALIZED 
r SET SHIFT COUNT TO 0 

? SHIFT 0 INTO (IX) 

RESULT 

J EXP OF ORIGINAL X 
; REMOVE BIAS 

EXP OF ORIGINAL Y 
J REMOVE BIAS 
»• ADD TWO EXPONENTS 





; CHECK EXPONENT CALCULATION FOR ERRORS 


HCHK 

JP 

PO, MiO 

? 

NO OVERFLOW 


ADD 

A, B 

r 

ADD SHIFT COUNT 


JP 

PO,MERR 

y 

NO SECOND OVERFLOW 


JP 

M20 



M10 

ADD 

A, B 

y 

ADD SHIFT COUNT 


JP 

PE, MERR 

7 

OVERFLOW 

M20 

ADD 

A, 80H 

* 

T 

RESTORE BIAS 


JP 

Z,UNDER 

1 

UNDERFLOW 


L.D 

(IX + 7 ) , A 

y 

STORE EXP 


JP 

RSTO 




*H DIVIDE 


I RECIPROCAL 

r 

REC CALL ST04 ! SAVE X 

LD HL,Ci 

CALL LD ; GET A 1 

LD IY , TEHP4 


? FLOATING POINT DIVISION ROUTINE 

r 

J INPUT! IX F'TS TO DIVIDEND 
IY PTS TO DIVISOR 

OUTPUT; IX PTS TO QUOTIENT (QUOT. REPLACES X, Y UNCHANGED) 
? ERRORS; DIVISOR ZERO 

; EXPONENT OVERFLOW OR UNDERFLOW 

? SPECIAL CASES? DIVIDEND ZERO 

; ALGORITHM; Y, 2Y,4Y,8Y ARE STORED IN TEMP. REGS. 

? X IS PLACED IN ACL. AC1 IS COMPARED TO THE TEMP. 

5 REGS. A BIT AT A TIME WITHIN EACH DIGIT, SUBTRACTING 

; WHEN TEMP < = ACL. RESULT BITS ARE ACCUMULATED IN C, 

; AND STORED A BYTE AT A TIME INTO (IX) 

,’ TIMING; APPROX. 22938 (20108 LOOP) 


yyyyyyyy 


yyyyyyyyyyyyyy 


y y y 


7 


; CHECK FOR ZERO OPERAND 


D IVI 

XOR 

A 


CP 

(IY+7) 


JP 

Z,DIVO 


CP 

(IX + 7 ) 

t 

JP 

Z , ZERO 

T 

; CAL 

CULATE 

AND SAVE 

r 

L.D 

A , ( IX ) 


XOR 

( IY) 


EX 

AF,AF » 


; STORE Y, 2Y, 4Y, AND 


; Y IS ZERO 
? X IS ZERO 
IGN OF RESULT 
; MS BYTE OF X 
5 SAVE IN A r 
8 Y IN TEMP REGISTER 



I 


; INITIALIZE ACl 


< 



PUSH 

T V 

JL f\ 

SAVE IX 

« 

CALL 

STA1 

? STORE X IN AC 1 


L D 

BC,701H 

; COUNTER IS 7, SHFT COUNT 1 


PUSH 

BC 

i SAVE COUNTERS 

i 

LD 

C. OFH 

r I NIT C 


LD ' 

HLrTEMPO 



CALL 

DVCP 

r COMPARE X TO Y (FRACTIONS) 

€ 

JR 

N C r D V 3 0 

? X >= Y 


POP 

BC 



DEC 

c 

; DECR SHIFT COUNT 

• 

PUSH 

BC 



JP 

DV20 

; BEGIN WITH SHIFT DIGIT 

€ 

7 

; DIVIDE LOOP 



r 

DV10 PUSH 

BC 

J SAVE COUNTERS 

i 

LD 

HLrAC1 



XOR 

A 



CALL 

SL.D1 

; SHIFT AC 1 LEFT 1 DIGIT 

1 

CALL 

DVDIG 

? GET 1 DIGIT RESULT 


DV20 LD 

HL rAC1 



KOR 

A 


€ 

CALL 

8LD1 

) SHIFT ACi LEFT A DIGIT 


DV30 CALL 

DVDIG 

} GET 2ND DIGIT RESULT 


LD 

A,C 

? COMPL RESULT IN A 

t 

CPI... 


? TRUE RESULT 


LD 

(IX > r A 

; STORE PAIR OF QUOTIENT DIG 


INC 

IX 


i 

POP 

BC: 



D JNZ 

DVi 0 

; COUNTER NOT 0 

€ 

r 

; CALCULATE EXPONENT OF 

RESULT 


r 

LD 

B,C 

? PUT SHIFT COUNT IN B 

f 

LD 

Ar(TEMPOE) 

? EXP OF ORIGINAL Y 


SUB 

8 OH 



LD 

L* r A 

J SAVE UNBIASED EXP 

1 

LD 

A,(A C1E) 

? EXP OF ORIGINAL X 


SUB 

80H 

} REMOVE BIAS 


SUB 

c 

} SUB TWO EXPONENTS 

€ 

POP 

IX 

J RESTORE IX 


JP 

MCHK 

? CHECK FOR ERRORS 


* H MULT f DIVI 

UTILITIES 



I * 

! MLDIG MULTIPLY DIGIT 

? INPUT; C CONTAINS DIGIT IN LOW 4 BITS 
i J (SAVE) IS PTR TO LS BYTE OF ACCUM 

; output; c is rotated 4 bits 

; PRODUCT OF mplicand with digit is 

I i ADDED TO ACCUM 

J TIMING 950 


« 

€ 

€ 

€ 

€ 

i 

i 

i 

1 

i 

« 

« 

i 

1 

i 

i 

« 


LD 

Dp 

HL. r TEMPOO 

; PT TO MPLICAND 

4 

CALL 

CrMLADD 

J ADD X 


LD 

D p 

HLrTEMP 10 

r* 


1 

r*. r*. 

CALL 

CrMLADD 

i ADD 2X 



t 4 


I 


I 







§ 

€' 

« 

€ 

« 

# 

« 

f 

« 

f 

f 

€ 


LD 

R R 

HL r TEMF‘20 

c 



CALL 

LD 

D £"» 

C , MLADD 

HL r TEMP30 

; ADD 

4K 

r*. r*. 

CALL 

RET 

L- 

C, MLADD 

; ADD 

SX 


MLADD ADD UTILITY FOR MULTIPLY 

INPUT! HL PTS TO NUMBER 

(SAVE) PTS TO LS BYTE OF ACCUM 
OUTPUT! 14 DIGIT NUMBER IS ADDED TO ACCUM 
TIMING 407 


LD 

DE, (SAVE ) 

; PT TO I..S BYTE 

AND 

A 

J CLEAR CY 

PUSH 

BC 

; SAVE BC 

LD 

B r 7 

f SET COUNTER 

LD 

A r (DE) 


ADC 

A,(HL) 


DAA 



LD 

(DE) r A 


DEC 

DE 


DEC 

HL 


DJNZ 

ML 40 


POP 

BC 


RET 



PRODUCES SIGN 

OF RESULT 


INPUT: IK PTS TO RESULT 

A’ BIT 7 IS SIGN BIT 

OUTPUT: IK PTS TO RESULT 

TIMING 74 


LD 

A » (IK ) 

J GET MS BYTE 

AND 

OFH 

? CLEAR UPPER DIGIT 

LD 

B, A 


EX 

AF.AF’ 

; GET NEW SIGN BIT 

AND 

80H 


OR 

B 


LD 

(IK) > A 

; REPLACE MS BYTE 

RET 




€ 

€ 

« 

€ 

« 


! DVDIG ROUTINE TO DIVIDE. 1 DIGIT 
f INPUT: AC 1 IS CURRENT REMAINDER (DIVIDEND) 

; TEMP REGS CONTAIN 1,2,4,8 *DIVISOR 

! OUTPUT: C CONTAINS QUOTIENT DIGIT IN LOW 4 
? (C ROTATED LEFT) 

? AC 1 CONTAINS REMAINDER << DIVISOR) 

! TIMING 1162 


DVDIG 

LD HL.TEMP3 

CALL DVCP 

; PT TO 

8Y 


JR C t DVD10 

LD HL t TEMP30 

! AC1 < 

8Y 


CALL DVSUB 

; 8Y < = 

ACi 

DVD 10 

RL C 

} COMPL 

RESULT BIT 


LD HL > TEMP2 

CALL DVCP 

; PT TO 

4Y 


f 


I 


€ 

i 

f 

f 

C 

< 

t 

< 

€ 

i 

i 

< 

i 

4 

4 

BITS 

4 

4 

i 

i 

i 






JR 

C r DVD20 

J ACl < 4Y 


LD 

HL t TEMP20 



CALL 

DVSUB 

; 4Y < ACl 

DVD20 

RL 

c 

; COMPl. RESULT BIT 


LD 

HL,TEMPI 

i PT TO 2Y 


CALL 

DVCP 



JR 

C r DVD3Q 

; ACl < 2Y 


LD 

HL»TEMPIO 



CALL 

DVSUB 

; 2Y < ACl 

DVD30 

RL 

C 

i COMPL RESULT BIT 


LD 

HL r TEMPO 

; PT TO Y 


CALL 

DVCP 

i COMPARE TO ACl 


JR 

C t DVD40 

; ACl < Y 


LD 

HL i TEMF'OO 



CALL 

DVSUB 

; Y < ACl 

DVD40 

EI¬ 

C 

COMPL RESULT BIT 


RE T 



T 

; dvcp 

ROUTINE TO COMPARE (DE) TO ACl 


f RETURNS 

CY SET IFF 

(DE) < ACl 

5 RETURNS 

2 SET IFF 

(DE) = ACl 

? TIMING APPROX 42 


7 

DVCP 

LD 

DErACl 

; PT TO ACl 


LD 

B r 7 

? SET COUNTER 

DVCIO 

LD 

A,(DE) 

5 GET BYTE 


CP 

(HL) 



RET 

NZ 

; < OR > 


INC 

HI¬ 

J INC PTRS 


INC 

DE 



DJN2 

DVCIO 

? COUNTER NOT 

T 

RET 




; DVSUB ROUTINE TO SUBTRACT (DE) FROM ACl 
J ASSUMES 14 DIGIT OPERANDS 

J TIMING 376 


DVSUB 

LD DE » AC10 

; PT TO ACCOM 


AND A 

LD B, 7 

? CLEAR CY 

DV40 

LD Ax(DE ) 

SBC Ax(HL) 

DAA 

LD (DE)r A 

DEC DE 

DEC HL 

DJNZ DV40 


r 

RET 


r 

i ROUTINE TO LOAD TEMP 

REGISTERS 

; TIMING 2014 


7 

HDTL 

LD DE,TEMPO 

) FIRST DEST 


CALL MOV 

} LOAD TEMPO 


LD HLxTEMPO 

; PT TO SRC 


RES 7 x(HL) 

x CLEAR SIGN 


L D B x 3 

5 SET COUNTER 


IN C 


IN C 


IN C 


; LOOP TO COPY TEMPN TO TEMPN+i & SHFT LET 1 BIT 





MDTL1 


LD 

A, B 

CALL 

MOV 

LD 

B, A 

AND 

A 

CALL 

SLB1 

INC 

HL 

D JNZ 

MDTL 

RET 



SAVE COUNTER 
HOVE TO NEXT TEMP 
RESTORE COUNTER 
CLEAR CY 

SHFT LEFT (HI. ALREADY OK) 
PT TO MSBYTE JUST SHIFTED 
COUNTER NOT 0 


J DETERMINE OVER OR UNDER 

7 

T 

HERR 


JP 

JP 


P t UNDER 
OVERA 


UNDERFLOW 


! DIVISION BY ZERO ERROR 

DIVO LD HL ,ERROR 

SET 4,(HL) 

JP INF 

*H FIX 

7 

; FIX CONVERTS FLOATING POINT NUMBER TO 
? 16 BIT 2'S COMPLEMENT INTEGER 

; (ROUNDING) 

f 

T 

i INPUT: ix PTS TO FL PT NUMBER 
OUTPUT: IX PTS TO INTEGER 
J ERRORS: OVERFLOW IF X > 32767 OR X 


•32768 


rrrrrrrrTrrTrrrrTfrrrrrrrrrrT 


FIX 

LD 

I Y r C „ 5 



CALL 

ADD 

? ADD .5 


CALL 

FIXA 

; CONVERT NUMBER 


JR 

N C » F X10 

J NO OVERFLOW 


LD 

HL t ERROR 



SET 

6 r ( HL. > 

; FLAG OVERFLOW 


LD 

HLr 7FFFH 

f 32767 


BIT 

7, (IX) 



JR 

ZrFXlO 

? POSITIVE SIGN 


L.0 

Hl. r 8000H 

? -32768 

F XID 

LD 

(IX) t L 

5 STORE INTEGER 


LD 

(IX+1) t H 



RET 



r 

; MAIN 

CONVERSION FOR 

FIX 

r 

} RETURNS 

CY SET IF 

OVERFLOW 

i RETURNS 

Z SET IF 

EXACT INTEGER 

; HL. 

IS LOW 16 BITS 

OF INT(X) FOR ALL : 

r 

FIXA 

XOR 

A 



LD 

H, A 

} CLEAR ACCUM 


LD 

L t A 



LD 

(SAVE) , A 

J CLEAR FLAG 

r 

; CALCULATE 

# DIGITS 1 

.EFT AND RIGHT OF Dl 

r 

LD 

A,(IX+7) 

; GET EXP 


SUB 

80H 

f SUBTRACT BIAS 





v 

/v 



JR 

CrFXAlOO 


JR 

Z t FXA100 


LD 

B * A 


LD 

A, 13 


SUB 

B 


JR 

C r FXA90 


LD 

C, A 

T 

J CONVERT LEFT DIGITS 

r 

PUSH 

IX 


POP 

IY 


LD 

A » ( IY ) 


INC 

IY 


JP 

FXA20 

F X A10 

CALL 

NZ/FXA110 


LD 

A r (IY ) 


INC 

IY 


PUSH 

RRA 

RRA 

RRA 

RRA 

AF 


CALL 

MAD 


DJNZ 

FXA15 


POP 

AF 


AND 

OFH 


JR 

FXA30 

FXA15 

CALL 

NZ/FXA110 


POP 

AF 

FXA20 

CALL 

MAD 


DJNZ 

FXA10 

? 

; CALCULATE 

OR OF ALL 

r 

XOR 

A 

FXA30 

LD 

B t C 


SRL 

B 


JR 

Z r F X A 5 0 

FXA40 

OR 

INC 

(IY) 

IY 


DJNZ 

FXA40 

: CONVERT TO 

2'S COMPL 

y 

FXA50 

BIT 

7 i (IX) 


JR 

Z t FXA80 


EX 

DE r HL 


LD 

HL, 0 


AND 

A 


JR 

Z > F X A 6 0 


INC 

DE 

FXA60 

SBC 

HL ( DE 


JP 

P r FXA90 

FXA70 

AND 

A 


L.D 

A,(SAVE) 


RLA 



RET 


FXA80 

BIT 

7 ? H 


JR 

Z > FXA70 


. i 

? X < 1 

; SAVE # LEFT DIGITS 


? X >= 10+13 
f SAVE # RIGHT DIGITS 

TO BINARY 

; SET POINTER 

; GET FIRST BYTE 


? OVERFLOW 

; GET NEXT TWO DIGITS 

5 SAVE THIS BYTE 
? GET LEFT DIGIT 


? CONVERT LEFT DIGIT 
? COUNTER NOT 0 
; GET RIGHT DIGIT 


; OVERFLOW 
; GET SAVED BYTE 
? CONVERT RIGHT DIGIT 
; COUNTER NOT 0 

DIGITS TO RIGHT OF DECML 

; initialize: a 

; USE OTHER COUNTER 
; DIVIDE BY 2 
? COUNT IS 0 
: OR 2 DIGITS 

: COUNTER NOT 0 

EMENT 


? POSITIVE SIGN 


; EXACT INTEGER 
? CORRECT MAGNITUDE 
; SUBTRACT FROM 0 
; OVERFLOW 
; SET Z OR NZ 
? GET FLAG 
; PUT IN CY 


; NO OVERFLOW 


€ 


i 

i 

< 

I 

i 

l 

i 

4 

I 

I 

PT I 

1 

I 

I 

I 

I 

I 

I 

I 

4 

1 


i 


V 

FXA90 

AND 

A 

T 



SCF 


1 



RET 



• 

F X A10 0 

LD 

Ai(IX+7) 

1 



JP 

FXA50 



F X A110 

L.D 

A,(SAVE) 

1 

€ 


OR 

8 OH 

r 



LD 

(. SAVE ) , A 


m 


RET 



% 

: MAD 
* 

ROUTINE 

! TO MULTIPLY 

m 

T 

MAD 

L.D 

DrH 

t 

i 



LD 

E t L 




ADD 

HL, HL. 

* 

y 

€ 


ADD 

HL , HL 




ADD 

HL. r DE 

J 



ADD 

Hl., HL 

y 

€ 


AND 

OFH 

y 



LD 

E, A 




LD 

D,0 


t 


ADD 

HL * DE 

y 



LD 

A r H 

7 

y 



AND 

OF OH 


C; 


RET 




SET Z OR NZ 
SET CY 

A IS 0 IFF FRACTION 0 

GET FLAG 
TURN IT ON 


ACCUM BY 10 AND ADD NEW DIGIT 

SAVE COPY OF HL. 

SHIFT LEFT 2 

ADD COPY 

SHIFT LEFT 1 

PUT NEW DIGIT IN DE 


ADD NEW DIGIT 
TEST HIGH BITS 


*H FLOAT 


? FLOAT CONVERTS 16 BIT 2’S COMPLEMENT INTEGER TO 
i FLOATING POINT FORMAT 

§ ; (FLOAT1 ASSUMES INPUT IS IN HL) 


« 

€ 

C 

i 

i 

i 

i 

< 

€ 

< 


5 INPUT: IK PTS TO 16 BIT INTEGER 

€ ; OUTPUT? IX PTS TO FLOATING PT NUMBER 

J ERRORS: NONE 



i 


7 7 7 * 

7 7 7 7 


€ 


I 

« 

€ 

« 

i 

t 

< 


; TEST INPUT 


FLOAT 

LD 

L r ( I X ) 

; LOAD 16 BIT INTEG 


LD 

H/(IX+1) 


FLOAT! 

PUSH 

HL 

; SAVE INTEGER 


CALL 

ZERO 

; CLEAR OUTPUT AREA 


POP 

HL 



LD 

A.H 

; TEST FOR ZERO 


OR 

L 



RET 

Z 

? INPUT WAS ZERO 

i CONVERT TO 

SIGN MAGNITUDE, SIGN IN A' 

y 

LD 

A i H 

5 GET SIGN BIT 


AND 

80H 



JR 

Z.FLIO 

; POSITIVE 


EX 

DE, HL 



LD 

HL r 0 



S B C 

HL , DE 

5 NEGATE 


y 

? CONVERT BINARY TO DECIMAL 


i 

€ 

< 

i 

l 

i 

i 

l 


« 


f 


PUSH 

IX 


POP 

HL 

5 PT TO OUTPUT AREA 

LD 

BC,iOOOO 


CALL 

FLDIV 

J DIVIDE BY 10000 

INC 

HL 

? ADVANCE OUTPUT PTR 

LD 

BC,1000 


CALL 

FLDIV 

J DIVIDE BY 1000 

LD 

BC,100 


CALL 

FLDIV 

; DIVIDE BY 100 

INC 

HL 

? ADVANCE OUTPUT PTR 

LD 

BC, 10 


CALL 

FLDIV 

; DIVIDE BY 10 

LD 

A, E 

; GET REMAINDER 

RL..D 


; ROT INTO OUTPUT 


T 

i NORMALIZE RESULT, STORE SIGN 

r 

LD C,S5H ; EXPONENT 

CALL NORML1 

JP RSTO J STORE SIGN 

r 

; ROUTINE TO DIVIDE DE BY BC 
J 4 BIT QUOTIENT IS ROTATED INTO (HL> 


FLDIV EX 

DE, HL 

; PUT ACCUM IN HL. 

XOR 

A 

i CLEAR A, CY 

DEC 

A 

? INITIALLY A IS -1 

FUJI VI SBC 

HL , BC 

J SUBTRACT DIVISOR 

INC 

A 

; INCR QUOTIENT 

JR 

NC, FL.DIV1 

J RESULT NOT NEGATIVE 

ADD 

HL , BC 

; RESTORE REMAINDER 

EX 

DE, HL. 

? SAVE ACCUM, GET PTR 

RLD 


5 STORE QUOTIENT 

RET 

*H CP 

f 

; CP FLOATING POINT COMPARE ROUTINE 
* 

T 

RETURNS 

CY SET IF X 

< Y 

} RETURNS 

Z SET IF X 

= Y 

; CPI ASS! 

f 

JMES MS BYTES 

FOR X,Y IN D,E 

T 

CP LD 

D , (IX > 

; LOAD MS BYTES 

LD 

E , < IY > 


CPI LD 

A, D 

5 GET SIGN OF X 

AND 

SOH 

(CY IS OFF) 

JR 

NZ,CPI 0 

; X IS NEGATIVE 

BIT 

7, E 


RET 

NZ 

; Y IS NEGATIVE 


; X AND Y ARE BOTH POSITIVE, COMPARE EXPONENT 

LD A,(IX + 7) 

CP (IY+7) 

RET NZ ; EXP'S NOT = 


COMPARE X FRACTION TO Y FRACTION 




LD 

A 7 D 

? MS BYTE. 


CP 

E 

? (SIGNS ARE EQUAL) 


RET 

NZ 



PUSH 

IX 

; SET UP F'TRS 


PUSH 

IY 



JR 

CF'20 


7 

; x I s 

NEGATIVE 7 CHECK 

SIGN OF Y 

7 

CF'10 

XOR 

E 

; TEST SIGN, SET NZ 


RLA 


J SIGN BIT IN CY 


RET 

c 

; Y WAS POSITIVE 

7 

; BOTH 

SIG N S 

NEGATIVE, 

TEST EXPONENTS 

T 

LD 

A, (IY + 7) 



CP 

(IX + 7 ) 


r 

RET 

NZ 


7 

; COMPARE Y ! 

t 

FRACTION TO 

X FRACTION 

7 

L.D 

A , E 



CP 

D 



RET 

NZ 



PUSH 

IY 

; SET UP F'TRS 


PUSH 

IX 


7 

; LOOP 

TO COMPARE 6 FRACTION BYTES 

7 

CF'20 

POP 

HL 

: GET F’TRS 


POP 

DE 



INC 

HL 



INC 

DE 



LD 

B , 6 

? COUNTER 

CP3G 

LD 

A, (DE) 



CP 

(HL ) 



RET 

NZ 



INC 

HL 



INC 

DE 



P JNZ 

CF'30 



RET 

*H NORMALIZE 

r 

; NORMLZ NORMALIZES AN ACCUM BY SHIFTING LEFT 

7 

! INPUT: IK PTS TO ACCUM (NORMLZ>NORML1) 

? HL. PTS TO ACCUM (NORMLZ) 

! C CONTAINS BIASED EXP (NORMLi,NORML 

; OUTPUT: NORMAL RETURN - NZ & NC 

J ZERO FRACTION - Z & NC 

f EXP UNDERFLOW - NZ & C 

? C CONTAINS UPDATED EXPONENT 

J (IX + 7 ) CONTAINS UPDATED EXPONENT 


NORMLZ 

NORMLI 


NORMLZ 


L.D C, (IX + 7 ) 
PUSH IX 
POP HL 
LD A 7 (HL > 
AND A 

JR NZ,NOR15 


; GET EXPONENT 
; F'T TO ACCUM 

? GET MS BYTE 


ALREADY NORMALIZED 




L.D 

B , 12 

; MAX COUNT 


X OR 

A 

5 CLEAR A 

NOR 10 

DEC 

c 

; DECR EXPONENT 


JR 

Z,N0R20 

? UNDERFLOW 


CALL¬ 

SLD1 

; SHIFT LEFT 


ING 

HL 

; PT TO 1ST BYTE 


LD 

A r ( HL ) 

? GET MS BYTE 


AND 

A 



JR 

NZ t NOR15 

; NORMALIZED 


DJNZ 

NOR 10 

5 COUNTER NOT 0 

NOR 15 

LD 

(IX + 7 ) ,C 

; STORE EXPONENT 


RET 



T 

UNDERFLOW, 

HL PTS TO 

MS BYTE» A CONTAINS MS 

7 

N0R20 

INC 

B 

; CALC BYTES LEFT 


SRI. 

B 


N0R30 

INC 

HL 

} ADVANCE PTR 


OR 

(HL) 

i COLLECT REMAINING BIT 


DJNZ 

N0R30 

? COUNT NOT 0 


RET 

C* f‘ p 1 

~r 

A. 

; FRACTION WAS ZERO 


o L- r 

RET 


; (UNDERFLOW) 


*H HALF 


r 

; HALF DIVIDES A FLOATING POINT NUMBER BY 2 
; INPUT MUST BE POSITIVE (0 IN SIGN DIGIT) 

r 

; INPUT: l>( PTS TO X (HALF) 

? OUTPUT: INPUT REPLACED BY X/2 
! FLAGS SET ACCORDING TO NORMLZ 

7 


AND 

A 

; CLEAR CY 

CALL 

SRB 

r DIVIDE FRACTION BY 2 

CALL 

NORMLZ 

; NORMALIZE RESULT 

RET 




*H SHIFT UTILITIES 

7 

r SLD SHIFT LEFT DIGIT 


INPUT: IX 

HL 
A 

OUTPUT: A 
HL 


PTS TO ACCUM (SLD) 

PTS TO ACCUM (SLD1 ) 
CONTAINS DIGIT TO SHIFT IN 
CONTAINS DIGIT SHIFTED OUT 
PTS TO i LESS THAN MS BYTE 


S i. D 

PUSH 

POP 

IX 

HL 

? 

SET PTR TO ACCUM 

SLD1 

PUSH 

LD 

BC 

BC i 6 

f 

SAVE BC 


ADD 

HL , BC 

7 

PT TO RIGHT END 


LD 

B r 7 

7 

LOOP COUNTER 

Si... DIO 

ELD 


7 

ROT LEFT 


DEC 

HL 

7 

PT TO NEXT BYTE 


DJNZ 

SLD10 

7 

COUNTER NOT 0 


POP 

RET 

BC 

* 

7 

RESTORE BC 

T 

; SRD 

SHIFT 

RIGHT 

DIGIT 






€ 

; INPUTS 

IK PTS 

TO 

ACCOM (SRD) 


r 


HL PTS 

TO 

ACCOM (SRDi ) 


t 


A CONTAIN 

S DIGIT TO SHIFT IN 

m 

l OUTPUTS 

A CONTAIN 

S DIGIT SHIFTED OUT 


t 


HL PTS 

TO 

1 MORE THAN LS BYTE 

€ 

r 

SRD 

PUSH 

IX 


S SET PTR TO ACCOM 



POP 

HL 




SRD1 

PUSH 

BC 


; SAVE BC 

€ 


L.D 

B , 7 


? COUNTER 


SRDIO 

RRD 



S ROT RIGHT 



INC 

HL 



i 


DUNZ 

SRDIO 


; COUNTER NOT 0 



POP 

BC 


J RESTORE BC 



RET 




€ 

r 






; SLB 

SHIFT 

LEFT BIT 



« 

1 

} INPUTS 

IX PTS 

TO 

ACCUM (SLB) 


f 


HL PTS 

TO 

ACCIJM (SLB1 ) 

t 

1 


CY CONTAIN 

S BIT TO SHIFT IN 

; OUTPUTS 

CY CONTAIN 

S BIT SHIFTED OUT 


r 


HL PTS 

TO 

1 LESS THAN MS BYTE 

« 

r 

SLB 

PUSH 

IX 


J SET PTR TO ACCUM 



POP 

HL 




SLB1 

PUSH 

BC 


; SAVE BC 

€ 


LD 

BC , 6 





ADD 

HL , BC 


? PT TO RIGHT END 



LD 

B ,7 


; LOOP COUNTER 

< 

S L BIO 

LD 

A » (HL. 5 


; GET BYTE 



ADC 

A, A 


; SHIFT LEFT 



DAA 




€ 


L.D 

(HL ) , A 


i REPLACE 



DEC 

Hl. 


} PT TO NEXT BYTE 



DJNZ 

SL BIO 


; COUNTER NOT 0 

« 


POP 

BC 


i RESTORE BC 



RET 




« 

? 

: SRB 

r 

SHIFT 

RIGHT BIT 




INPUT; 

IX 

PTS TO ACCUM 

(SRB) 

* 

r 


HL 

PTS TO ACCUM 

(SRB1 ) 

? 


CY 

CONTAINS BIT 

TO SHIFT IN 

i 

OUTPUT; 

CY 

CONTAINS BIT 

SHIFTED OUT 

f 

¥ 


HL 

PTS TO 1 MORE 

THAN LS BYTE 


RB 

PUSH 

IX 

7 

SET PTR TO ACCUM 


POP 

HL 



RBI 

PUSH 

BC 

» 

7 

SAVE BC 


LD 

8,7 

T 

COUNTER 

RBI 0 

LD 

A, (HL) 

f 

GET BYTE 


RRA 


r 

SHIFT 


RR 

C 

f 

SAVE CY IN C 


BIT 

7, A 




JR 

Z , SRB20 

r 

0 SHIFTED INTO UPPER 


SUB 

3 OH 

r 

CORRECTION 

RB20 

BIT 

3, A 




JR 

Z,SRB30 

7 

0 SHIFTED INTO LOWER 


SUB 

3 

r 

CORRECTION 


€ 



SRB30 


RL. C 
LD (HL) , A 
INC HI... 

DJNZ SRB10 
POP BC 
RET 

*H HOVE UTILITIES 


GET GY 
STORE BYTE 

COUNTER NOT 0 
RESTORE BC 


MOVE ROUTINES 


ST 04 

LD 

DE, TEMP4 




JP 

STO 



ST 05 

LD 

DE, TEMPS 




JP 

STO 



ST 06 

LD 

DE,TEMP6 



STO 

PUSH 

IX 




POP 

HL. 

5 SET SRC PTR 


MOV 

LD 

L.DIR 

RET 

BC , 8 



LD 

PUSH 

IX 




POP 

DE 

; SET DEST PTR 



JP 

MOV 



1 

) SPEC 

IALIZED LOAD ROUTINES, LOAD AC1 AND 

AC2 

T 

S T A1 

PUSH 

IX 

; SOURCE WILL BE 

(IX) 


JR 

ST 1 



STA2 

PUSH 

IX 

! SOURCE WILL BE 

(IX) 


JR 

ST 2 



STYA1 

PUSH 

IY 

J SOURCE WILL BE 

(IY) 

ST1 

POP 

HL 



M VI 

LD 

DE,AC1 

; DEST IS ACi 



CALL 

MOV 




LD 

HL,ACl 




LD 

D,(HL) 

} SAVE MS BYTE IN 

D 


RES 

7, (HL) 

J CLEAR SIGN BIT 



RET 




STYA2 

PUSH 

IY 

; SOURCE WILL BE 

(IY) 

ST2 

POP 

HL. 



MV2 

ID 

D E , A C 2 

; DEST IS AC2 



CALL 

MOV 




LD 

HL,AC2 




LD 

E,(HL) 

J SAVE MS BYTE IN 

E 


RES 

7, (HL) 

? CLEAR SIGN BIT 



RET 

*H BOOLEAN 
! LOGICAL OPERATORS 


I OR 

LD 

A , ( IX ) 


OR 

(IX +1 ) 


OR 

(IY) 


OR 

(IY+1) 


JP 

INF I 

I AND 

LD 

A, (IX 5 


OR 

(IX+1) 


JP 

Z,IZERO 


LD 

A,(IY) 


OR 

(IY+1) 


IS ZERO 


VALUE 



JP 

INEI 

INOT 

LD 

A t <IX ) 


OR 

( IX+1 ) 


JR 

IEQ1 

OR 

LD 

A, (IX + 7) 


OR 

( IY+7 ) 


JP. 

NE1 

AND 

L.D 

A r < IX+7) 


AND 

A 


JP 

Z,ZERO 


I..D 

A, ( IY+7) 


AND 

A 


JR 

NEi 

NOT 

LD 

A, (IX+7) 


AMD 

A 


JR 

E Q 3 

r 

RELATIONAL 

OPERATOR 

T 

SGT 

CALL 

SCP 


JR 

IGTi 

SLE 

CALL 

SCP 


JP. 

ILE 1 

SLT 

CALL 

SCP 


JR 

ILT 1 

SGE 

CALL 

SCP 


JR 

IGE 1 

SEQ 

CALL 

SCP 


JR 

I EG 1 

SNE 

CALL 

SCP 


JR 

INEi 

IGT 

CALL. 

I CP 

IGTi 

JR 

Zt IZERO 


JR 

C,I ZERO 


JR 

I ONE 

ILE 

CALL 

I CP 

I LEI 

JR 

Z, IONE 


JR 

C,IONE 


JR 

IZERO 

ILT 

CALL 

I CP 

ILT1 

JR 

C, IONE 


JR 

IZERO 

IGE 

CALL 

I CP 

I GE1 

JR 

C , IZERO 


JR 

IONE 

JEQ 

CALL 

I CP 

IEQ1 

JR 

Z , TONE 


JR 

IZERO 

I ME 

CALL 

I CP 

INEi 

JR 

Z,IZERO 


JR 

IONE 

GT 

CALL 

CP 

G T1 

JR 

Z t ZERO 


JR 

CfZERO 


JR 

ONE 

LE 

CALL 

CP 

LEI 

JR 

Z t ONE 


JR 

C i ONE 


JR 

ZERO 

LT 

CALL. 

CP 





i 


i 

i: 

€ 

€ 

€ 

« 

§ 

i 

i 

t 

< 

f 

i 

c 

< 

€ 


1. T1 

JR 

C, ONE 


JR 

ZERO 

GE 

CALL 

CP 

GE1 

JR 

C , ZERO 


JR 

ONE 

EQ 

CALL 

CP 

ECU 

JR 

Z r ONE 


JR 

ZERO 

NE 

CALL 

CP 

NE1 

JR 

Z r ZERO 


JR 

0 N E 

*H C 

0N8TANT 

FUNCTIONS 

r 

; PRODUCE UNITY RESULT 

r 

TONE 

LD 

(IX) , 1 


L.D 

(IX + 1),0 


RET 


ONE 

LD 

B» 7 


CALL 

ZEROl 


LD 

(HL)/81H 


LD 

(IX ) , 1 


RET 



? PRODUCE ZERO RESULT 


I ZERO 

LD 

(IX ) , 0 


LD 

RET 

CIX+1),0 

ZERO 

LD 

Br 8 

ZERO! 

PUSH 

IX 


POP 

HL 

ZER1D 

LD 

(HL ) rO 


INC 

HL. 


DJNZ 

ZER10 

r 

RET 



l PRODUCE INFINITE RES 


INF 

PUSH 

IX 


POP 

HL 


RL 

(HL ) 


LD 

(HL)r18 


RR 

(HL ) 


INC 

HL 


LD 

A r 99H 


LD 

B,6 

INF10 

LD 

(HL),A 


INC 

Hl. 


DJNZ 

INF 10 


LD 

(HL),OFFH 


RET 


T 

! FLAG 

OVERFL 

.OW ERROR» 

r 

OVERA 

EX 

AF»AF’ 


L.D 

(IX) ,A 

OVER 

LD 

HL,ERROR 


SET 

6,(HL) 


JP 

INF 


; PT TO RESULT 


T 

SET PTR 

; GET OLD SIGN 
? REPLACE SIGN 


; HIGHEST EXP 

STORE INFINITY 

; GET SIGN OF RESULT 
? STORE IT 

f FLAG OVERFLOW 


« 

« 

i 

« 

i 

« 

i 

« 

< 

« 

i 

i 

€ 

€ 

i 

i 

i 

l 

I 

i 

i 


i 



; FLAG UNDERFLOW ERROR, STORE ZERO 

UNDER LD HL,ERROR 

SET 5, (Hi. ) ; FLAG UNDERFLOW 

JP ZERO 

*H CONSTANTS AND WORK AREAS 
Cl DEFW 1 

DEFW 0 
DEFW 0 
DEFW 810QH 
CO DEFW 0 

DEFW 0 
DEFW 0 
DEFW 0 

C...5 DEFW 5 ; . 5 

DEFW 0 
DEFW 0 
DEFW 8000H 

SEED DEFW 0 f 0 

DEFW 0 
DEFW 0 
DEFW 0 

SEED10 EQU SEED+6 

CRNDA DEFW 1403H ? 31415 92853 581 

DEFW 9215H 
DEFW 3565H 
DEFW 8081H 

CRNDC DEFW 7102H i 27182 81828 459 

DEFW 8182H 
DEFW 8482H 
DEFW 8059H 

CRNDCO EQU CRNDC+8 
SAVE DEFS 2 
AC1 DEFS 8 

AGIO EQU AC1+6 

ACIE EQU AC1+7 ? EXPONENT BYTE 

AC2 DEFS 8 

A C20 EQU AC2 +6 

AC2E EQU AC2+7 f EXPONENT BYTE 

TEMPO DEFS 8 
TEMPOO EQiJ TEMPO + 6 
TEMPOE EQU TEMPO+7 
TEMPI DEFS 8 
TEMP10 EQU TEMPi+6 
TEMP2 DEFS 8 
TEMP20 EQU TEMP2+8 
TEMP3 DEFS 8 
TEMP30 EQU TEMP3+6 
TEMP4 DEFS 8 
TEMPS DEFS 8 
TEMP6 DEFS 8 
1 

; REFERENCES 

* 

T 

EXTERNAL ERROR,ICP,SCP 

GLOBAL ADD , SIJB, RSUB , MULT , REC , DIVI 

GLOBAL MLDIG,MDTL 

GLOBAL FIX,FIXA,FLOAT,FLOAT1,CP,MAD 

G LOBAL GT,LE,LT,GE,EQ,NE,OR,AND,NOT 





GLOBAL. 

GLOBAL 

GLOBAL 

GLOBAL 

GLOBAL 

GLOBAL 

GLOBAL 

#H DECIMAL TO BINARY 


IGT,ILErILT,IGE,IEQ»INE,IOR,IAND,INOT 
SGT , SLE , S'LT , SGE , SEQ , SNE 

SEED, SEED 10 , CRNDA, CRNDCO ,C_5, SAVE , C1, CO 

OVER,UNDER,ZERO,ONE,INF 

NORMLZ, NOR ML 1 , N0RML.2 , SLD , SRB , HAL F 

ST04 , ST05, ST06 , STO , L.D, SRD 

TEMP4,TEMPS,TEMP6,AC 1,AC 10,AC2,AC2 0 


! DECIMAL TO BINARY CONVERSION 
i INPUT: IK POINTS TO OUTPUT AREA 

; OUTPUT: IK POINTS TO FL PT NUMBER 

? ERRORS: NO DIGITS IN FRACTION PART 
5 NO DIGITS IN EXPONENT 

J EXPONENT OVERFLOW 

5 EXPONENT UNDERFLOW 

? REGISTER ASSIGNMENT IN MAIN LOOP: 

J C - DECIMAL EXPONENT 

; B - FLAGS: BIT 6 - DIGITS FOUND 

: BIT 7 - FRACTION OVRFLW 


T 

t 

t 

A 

A' 

- NEW DIGIT 

- SIGN IN BIT 7 

r r y y r r 

T 1 ? J 1 J 

rrrtrrTrrtr 

rrTifrTtrrrrrrrrrrTTT 

DTOB 

CALL 

LMARK 

} MARK BEGINNING OF 


CALL 

ZERO 

5 CLEAR RESULT 


X 0 R 

A 

} CLEAR A 


LD 

B, A 

J CLEAR FLAGS 


LD 

Cr80H+13 

J INIT EXPONENT 


EX 

AF , AF ' 

5 PUT AWAY SIGN 


CALL 

LIN 

,’ GET FIRST CHAR 


CP 

' + ’ 



JR 

ZrDIO 

; SIGN IS + 


CP 

r r 



JR 

NZ , D20 

; NO SIGN 


EX 

AF, AF' 



OR 

80H 

J FLAG NEGATIVE SIGN 


EX 

AF,AF' 


; loop 

TO PROCESS CHARS 

BEFORE DECIMAL POINT 

DIO 

CALL 

LIN 

J GET NEXT CHAR 

D20 

CP 

r r 

» 



JR 

Z , D 4 0 

; DECML PT FOUND 


CALL 

DIGIT? 



JR 

C, D50 

: NOT A DIGIT 


SET 

6, B 

? FLAG DIGIT FOUND 


BIT 

7, B 



JR 

NZ , D30 

? OVERFLOW 


CALL 

CNV 

; CONVERT THIS DIGIT 


JP 

DIO 


D 30 

INC 

C 

; INC DECML EXP 


JP 

DIO 


r 

: loop 

f 

TO PROCESS CHARS 

AFTER DECIMAL POINT 

t 

D 40 

CALL 

LIN 

; GET NEXT CHAR 


CALL 

DIGIT? 



JR 

C, D50 

: NOT A DIGIT 



SET 6. B 

; FLAG DIGIT FOUND 

BIT 7, B 


JR NZ.D40 

; OVERFLOW 

DEC C 

i DEC DECML EXP 

CALL CNV 

; CONVERT THIS DIGIT 

JP D40 



? CHANGE FRACTION TO FLOATING POINT FORMAT 


BIT 

6. B 


JP 

Z.LREST 

? NO DIGITS FOUND 

LD 

(SAVE).A 

? SAVE THIS CHAR 

CALL 

NORML 1 


JR 

Z.D130 

; ZERO FRACTION PART 

LD 

A.(SAVE) 

i GET CHAR AGAIN 


J PROCESS EXPONENT 
; REGISTERS; 

B - BIT 7 IS SIGN 

f BIT 6 IS DIGIT-FOUND FLAG 

I BIT 5 IS DIGIT-OVERFLOW FLAG 

? C - IS EXP FROM FRACTION PART 

; HL. - IS ACCLIM FOR D--TO-B CONVERSION 

r 

LD HL>0 ; INITIALIZE REGISTER 

L.D B r L. 

CP ' E'-’0’ 

JR NZ r D9D ; NOT AN ’ E’ 

! CONVERT SIGN 

f 

CALL LIN 

CP r + ' 

JR Z r D60 ; SIGN IS + 

CP 

JR NZ.D70 ; NO SIGN 

SET 7 r B ; FLAG NEC. SIGN 

T 

i CONVERT DIGITS 

D60 CALL LIN 

D70 CALL DIGIT? 

JR C.D80 ; NOT A DIGIT 

SET 6,B ; FLAG DIGIT FOUND 

BIT 5. B 

JR NZ.D100 5 OVER OR UNDERFLOW 

CALL MAD ; CONVERT THIS DIGIT 

LD ArH J CHECK FOR OVERFLOW 

AMD A 

JP Z.D60 ? HL < 256 

SET 5»B I SET FLAG 

JP D60 

t 

T 

; TEST FOR ILLEGAL EXPONENT PART 
D80 BIT 6. B 

JR Z.D120 J NO DIGITS FOUND 


; CONVERT HL. ACCUM TO 2' S COMPLEMENT 


D 90 BIT 7,B 

JR Z > D92 
EX DE , HL. 

L D HL y 0 

AND A 
SBC HL,DE 

! ADD PREVIOUS EXPONENT (C 
J CHECK FOR OVERFLOW 

I 

D92 LD B y 0 

AND A 

ADC HL y BC J 

JP Z»UNDER ; 

I. D A , H 

AND A 

JR NZ y D96 ; 

r 

; STORE SIGN AND EXPONENT 

* 

7 

EX AFiAF’ 5 

OR (IX) 

L D (IX > »A 5 

LD (IX + 7)iL t 

JP RTN 

T 

r 

; DETERMINE SIGN OF EXP 

7 

D96 BIT 7iH 

,JR D105 

V 

7 

; FINISH SCANNING DIGITS 


SIGN IS + 


ADD TWO EXPONENTS 
UNDERFLOW 


EXP < -128 OR > 127 


GET SIGN 

STORE SIGN 
STORE EXP 


D100 

CALL 

LIN 

? 

FINIS 

H SCANNING DIGITS 


CALL 

DIGIT? 





JR 

NCt D100 

y 

DIGIT 

IS FOUND 

T 

DETERMINE 

FROM SIGN 

BIT 

WHETHER OVER OR UNDER 

7 

BIT 

7, B 




D 105 

JR 

NZ,D.11.0 

* 

7 

SIGN 

IS NEGATIVE 


CALL 

OVER 

t 

7 

FLAG 

ERROR 


JP 

RTN 




DUO 

CALL 

UNDER 

7 

FLAG 

ERROR 


JP 

RTN 




7 

; ILLEGAL NUMBER 




7 

D 120 

LD 

HL,ERROR 





SET 

7 y (HL) 





JP 

RTN 




7 

; ZERO 

t 






T 

D130 

CALL. 

ZERO 





JP 

RTN 





*EJECT 
1 

} SCANNER ROUTINES 


i 


RTN 


CALL LBACK i BACK UP POINTER BY 1 

JP RC2 ; GO TO SPECIAL RETURN 


J DIGIT? ROUTINE CONVERTS 
i RETURNS CY ON IF NOT 


ASCII CHAR TO 0-9 VALUE 
A DIGIT 


DIGIT? SUB '0' 



RET C 

? WAS < 30H 


CP 10 

CCF 

! CY SET IF :: 

> 9 

RET 




i CNV ROUTINE TO MULTIPLY DECIMAL ACCUM i BY 10 AND ADD 


CNV CALL 

SL.D 

? SHIFT A INTO ACCUM 

LO 

Ar(IX) 

f GET HIGH BYTE 

AND 

A 


RET 

Z 

! < 13 DIGITS 

SET 

7»B 

J FLAG DIGIT OVERFLOW 

RET 



*H BINARY TO 

DECIMAL 



J BCD TO DECIMAL CONVERSION 


INPUT? IK PTS TO FL PT NUMBER 
HL PTS TO OUTPUT AREA 
OUTPUT? FL PT NUMBER IS MODIFIED 
OUTPUT AREA IS 19 BYTES? 

1 SIGN (SPACE OR -) 

13 DIGITS OF FRACTION 
1 " E " 

1 SIGN OF EXPONENT (+ OR -) 

3 DIGITS OF EXPONENT (LEADING ZEROS) 
EXPONENT IS FOR DECIMAL PT AFTER FIRST DIGIT 
A IS 2'S COMPLEMENT EXPONENT (-128<=A<=126) 


? CONVERT SIGN 


BTOD 

BIT 

7 , (IX > 

? TEST SIGN 


LD 

(HL),' ' 

? DEFAULT SIGN 


JR 

Z,BD10 

? + SIGN 


LD 

(HL), 

? -- SIGN 

BD10 

INC 

HL 


r 

; UNPACK 

: FRACTION PART 


r 

EX 

DE , HL 

; OUTPUT PTR IN DE 


CALL 

SID 

? SHIFT OUT SIGN BYTE 


LD 

B, 13 

» COUNTER 


LD 

A , ' 0 ' 


BD20 

CALL 

SLD 

? GET A DIGIT 


LD 

(DE),A 

; OUTPUT IT 


INC 

DE 



D JNZ 

BD20 

? COUNTER NOT 0 


T 


CONVERT EXPONENT SIGN 


GET SIGN MAGNITUDE EXP 





EX DE , HI... 

LD (HL),'E' 
INC HI 
LD (HL),'+’ 
L.D A, (IX + 7 ) 
AND A 
JR Z,BD30 
SUB 81H 
JR NC.BD30 
LD (HL ) r’~ 7 
NEC 

r 

i CONVERT HUNDREDS 


BD30 

INC 

HL 


L D 

(HL),'O’ 


CP 

100 


JP 

C , B D 4 0 


LD 

(HL) , ' 1' 


SUB 

100 

J DIVIDE 

: BY 

10 

r 

BD40 

INC 

HL 


LD 

B , 2 F H 


LD 

C, 10 

BD50 

SUB 

c 


INC 

B 


JR 

NC,BD50 


ADD 

A, C 


LD 

(HL) , B 


INC 

HL 


OR 

'0' 


LD 

(HL) , A 


LD 

A,(IX+7) 

*H RND 

SUB 

RET 

80H 

t i y ? ? 7 y y 

? r r y 

rtlrlr tr? 


OUTPUT PTR IN HL 


DEFAULT SIGN 
GET BIASED EXP 

ZERO EXPON 

REMOVE BIAS , DECR BY 1 
POSITIVE EXP 
- SIGN 

CALC MAGNITUDE 


DEFAULT DIGIT 
EXP < 100 


QUOTIENT COUNTER 
DIVISOR 
SUBTRACT 10 
INC QUOTIENT 
RESULT NOT NEGATIVE 
RESTORE REMAINDER 
STORE QUOTIENT DIGIT 


STORE REMAINDER DIGIT 
GET EXPONENT 
REMOVE BIAS 


; RANDOM NUMBER GENERATOR 

i INPUT: IX PTS TO OUTPUT/WORK AREA 
f SEED CONTAINS PREVIOUS RANDOM NUMBER 

l OUTPUT: IX PTS TO RANDOM NUMBER 
; SEED CONTAINS NEW NUMBER 

J ERRORS: NONE 

J ALGORITHM: ARITHMETIC DONE MODULO 10+13, 

? X :” 3141592653581 * X + 2718281828459 

? RESULT IS DIVIDED BY 10+13, AND NORMALIZED 


? INITIALIZE REGISTERS 

r 

RND LD IY, CRN DA ; PT TO MPLIER 

PUSH IX 
POP DE 
LD HL.,6 
ADD HL,DE 


PT TO LS BYTE OF (IX) 



LD (SAVE),HL 
LD HL i SEED 
CALL HDTL 
CALL ZERO 
LD B , 7 
JR RN20 


? SAVE FOR MLDIG 
J FT TO LAST NUMBER 
f SAVE MULTIPLES IN TEMPO-3 
: CLEAR WORK AREA 
; SET MULTIPLY LOOP COUNTER 


MULTIPLICATION, USING (IX) AS ACCUM, (IY) AS 
MULTIPLIER, (SEED) (SEED=IY-8> AS PRODUCT 


R N10 


RN20 


RN30 


MULTIPLY 

BY 31415 

92653 581 (CRNDA) 

PUSH 

AF 

i SAVE DIGIT SHIFTED OUT 

CALL 

MLDIG 

; MULTIPLY BY 2ND DIGIT 

XOR 

A 


CALL 

SRD 

; SHIFT ACCUM RIGHT 

RLA 


? PUT DIGIT IN HIGH 4 BITS 

RL.A 



RLA 



RLA 



AND 

OFOH 


LD 

C, A 

? SAVE UPPER DIGIT 

POP 

AF 

? GET LAST DIGIT 

OR 

c 

} COMBINE 

LD 

(IY-1),A 

} STORE IN SEED 

LD 

C,(IY+6) 

; GET NEXT 2 MPLIER DIGITS 

DEC 

IY 

? DECR PTR 

CALL 

MLDIG 

? MULTIPLY BY 1ST DIGIT 

XOR 

A 


CALL 

SRD 

f SHIFT ACCUM RIGHT 

DJNZ 

RNIO 


LD 

(IY-1),A 

; STORE LAST DIGIT 

D 27.182 

81828 459 

(CRNDC) TO (SEED) 

LD 

B, 7 

; LOOP COUNTER 

LD 

HL,CRNDCO 

} PT TO CONSTANT 

LD 

DE,SEED10 

; PT TO L.S BYTE OF SEED 

AND 

A 

! CLEAR CY 

iD LOOP 



i LD 

A,(DE) 

; GET 2 DIGITS 

ADC 

A, (HL.) 

} ADD 2 DIGITS 

DAA 



LD 

(DE) , A 

? STORE SUM 

DEC 

DE 


DEC 

HL 


DJNZ 

RN30 

? COUNTER NOT 0 

:AL.E AND 

NORMALIZE 

RESULT 

INC 

DE 

: PT TO MS BYTE 

AND 

OFH 

} CLEAR SIGN DIGIT 

LD 

(DE) , A 


EX 

DE, HL 

5 PUT PTR IN HL 

CALL 

LD 

; COPY SEED TO (IX) 

LD 

CrSOH 

; SET EXPONENT 

CALL 

N0RML1 

} NORMALIZE 


RET 


; RANDOMIZE 


USES REFRESH REGISTER TO MODIFY SEED 


y yyyyyyyyyy 

RANDMZ LD 
OR 
DAA 
L.D 
RET 

*H INT 


yyyyyyyyyyyyyyyyyyyyyyyyyrfyyyyyy 


A»R ; GET REFRESH REGISTER 

A 

; CONVERT TO DECIMAL 

(SEEDIO)r A 


* ¥ % * 
r r 7 r 


; ; ; ; 


; INTEGER FUNCTION -- FLOOR FUNCTION 
i INPUT,OUTPUTS IX POINTS TO FL PT NUMBER 


? y y y y r y 

y y r y i 

yyyyyyyyyyy 

TtJ717T17f7777T7?y77r 

y y y y y y y 

INT 

X 0 R 

A 




CP 

(IX + 7) 




RET 

Z 

? X WAS ZERO 



LD 

A, SDH 




SUB 

(IX + 7 ) 

; # DIGITS RIGHT OF 

PT 


RET 

c 

? X >= 10+13 



RET 

“7 

il. 

} X >= 10+12 



CP 

13 




JR 

NCrINT40 

; x < 1 



LD 

B, A 



? 

? SET POINTER TO 1ST DIGIT RIGHT OF PT 


T 

LD 

A, 14 




SUB 

B 

; i+# DIGITS LEFT OF 

PT 


LD 

E r A 




SRL 

E 

; # BYTES TO SKIP 



LD 

D r 0 




PUSH 

IX 

? GET F'TR 



POP 

HL 




ADD 

HL., DE 



T 

! CLEAR 

FIRS 

T DIGIT IF 

A IS ODD 


t 

AND 

1 




JR 

Z, INTIO 

? DEC PT IS BETWEEN 

BYTES 


LD 

A,(HL) 

; GET BYTE WITH DEC 

PT 


LD 

C r A 

; SAVE RIGHT HALF 



AND 

OFQH 

? CLEAR RIGHT HALF 



LD 

(HL) , A 

5 REPLACE IT 



INC 

HL 




LD 

A y C 

; GET RIGHT HALF 



AND 

OFH 



t 

} CLEAR 

REMAINING BYTES 

, OR-ING EACH INTO A 

FIRST 

INTiO 

SRL 

B 

J CALC # BYTES 



JR 

Z,INT30 

r NONE TO CLEAR 


INT20 

OR 

(HL ) 




LD 

(HL) r 0 




INC 

HL. 




DJNZ 

INT20 




r 

; POSSIBLE CORRECTION FOR NEGATIVE NUMBERS 



INT30 BIT 

7 , (IX ) 



RET 

"7 * 

a- r 

POSITIVE 


AND 

A 



RET 

Z ? 

INPUT WAS AN INTEGER 


LD 

1Y , C1 



CALI. 

. SUB ; 

SUBTRACT 1 TO CORRECT 


RET 



t 

i 

X < 1, NONZERO? OUTPUT 

0 IF +, -1 IF - 

r 

INT4Q BIT 

7, (IX) 



JP 

Z r ZERO ; 

POSITIVE FRACTION 


CALI. 

. ONE 



SET 

7 y (IX) ; 

NEGATIVE 1 


RET 




SHORT FUNCTIONS 


r 

i 

BRK SEPARATES INTEGER AND FRACTION PARTS 

r 

T 

INPUT? 

IX PTS TO X 


J 

OUTPUT? 

IX PTS TO FRACTIONAL PART OF X (PO 


BRK 


HI. CONTAINS 2' S COMPL. INTEGER PART OF X, 
LOW ORDER 16 BITS VALID IF X < 10+13 
CY SET IF INTEGER WON'T FIT IN 16 BITS 
Z SET IF EXACT INTEGER 
A' IS ORIGINAL SIGN 

; GET SIGN 

; COPY X 

i TAKE INTEGER PART 

} CALC. FRACTION PART 
i SAVE PTR 

? CALC. BINARY INTEGER PART 

; GET ORIGINAL SIGN 


LD 

A » (IX ) 

PUSH 

AF 

CALL 

ST 04 

CALL 

I NT 

LD 

IY.TEMP4 

CALL 

RSUB 

PUSH 

IX 

LD 

IX r TEMP4 

CALL 

FIXA 

POP 

IX 

EX 

AF,AF ' 

POP 

AF 

EX 

AF > AF' 

RET 

r r r r J 

1 M Jl 1 M 


A B S0 L UT E V A L UE F U NOTION 

INPUT t OUTPUT i IX POINTS TO FL PT NUMBER 


ABS RES 1, (IX) 

RET 


J TURN OFF SIGN BIT 


! SIGN FUNCTION 

5 INPUTrOUTPUT: IX POINTS TO FL PT NUMBER 


SGN X 0 R A 

CP (IX+7) 
RET Z 
LD A,(IX) 


i CLEAR A 

f ZERO INPUT IS NOP 
? GET SIGN 



I 


« 

t 

i 

« 

€ 


AND 

80H 


CALL 

ONE 

J OUTPUT 1 

OR 

« IK) 


LD 

(IK ) , A 

: SIGN OF INPUT 

RET 




! NEGATE FUNCTION 

; INPUT,OUTPUT; IK POINTS TO FL PT NUMBER 


NEG 

LD 

A,(IK) 

} GET SIGN BIT 


XOR 

80 H 

; COMPLEMENT IT 


LD 

(IK ) , A 

r RESTORE 


RET 




*H INTEGER ADD AND SUB 


5 IADD, I SUB INTEGER ADD AND SUBTRACT 


I 

< 

i 

< 

« 

i 


* 

f 

1 

t 

i 

€ 

< 

« 

« 

« 


? 

INPUT: IK 

PTS 

TO X 


r 

IY 

PTS 

TO Y 


7 

OUTPUT: IK 

PTS 

TO SUM OR 

DIFFERENCE 

7 

ERRORS: OVERFLOW 

(POSITIVE 

OR NEGATIVE) 


irrttyrt 

* Y » » 

7 r t t 

777777777 

7777777777777777777 

I ADD 

CALL 

ASLOAD 

J LOAD ARGUMENTS 


ADC 

HL , DE 



JP 

PO,ASSTO 

; NO OVERFLOW 

ASOV 

UP 

PrNINF 

J OVFLW SIGN POS 


JP 

PINF 


I SUB 

CALL 

ASLOAD 

J LOAD ARGUMENTS 


SBC 

HLrDE 



JP 

PO,ASSTO 

; NO OVERFLOW 

r 

JR 

ASOV 


7 

; ASLOAD 

LOADS ARCS 

TO HL, DE, RESETS C 

7 

ASLOAD 

LD 

Hr <IX+1 ) 



LD 

L , (IX ) 



LD 

D j (IY+1 ) 



LD 

E, (IY) 



AND 

A 



RET 



7 

; ASSTO 

f 

STORES HL. AS 

RESULT 

T 

ASSTO 

LD 

(IK), L 



LD 

(IX+1) , H 



RET 




*H INTEGER MULTIPLY 


77777777777777777 


77777 


7 r ? y 7 r 


« 

« 

€ 

< 

i 

< 

i 

€ 

i 

< 

i 


« 


! IMUL INTEGER MULTIPLY 

Y 

7 

} INPUT; IK PTS TO K 

J IY PTS TO Y 

; OUTPUT: IK PTS TO PRODUCT 

; ERRORS: OVERFLOW (POS OR NEG) 

» REGISTERS: DE IS MULTIPLICAND 


i 

< 

« 


t 


i 









AC IS MULTIPLIER AND ACCUM EXTENSION 
HL IS ACCUM 
B IS LOOP COUNTER 


JtTTTtrtrTTT 

yyyyyyyyy 

77777777777777777777777777 

IMUL CALL 

MDLOAD 

7 I- 0 A D A R G b 

LD 

A, B 

! X IN AC 

T 

i 16 BIT BY 

16 BIT UN 

SIGNED MULTIPLY 

y 

L.D 

HL.,0 

5 CLEAR ACCOM 

LD 

B, 16 

i SET COUNTER 

RRA 

C 

N C , IM 2 0 

; INITIALIZE CARRY 

K ft 

IMiO JR 

; MPLIER BIT IS ZERO 

ADO 

HL, DE 

) ADD MULTIPLICAND 

IM20 RR 

H 

; ROT ACC-MPLICAND INTO 

RR 

L 


RRA 



RR 

C 


D JNZ 

IMIO 

; COUNTER NOT 0 

r 

; TEST FOR > 

16 BIT PRODUCT 

' LD 

B, A 

; PRODUCT IN BC 

LD 

A, H 


OR 

L.. 


JR 

* 

NZ,INF? 

; OVERFLOW 

T 

TEST FOR >15 BITS, 

AND STORE RESULTS 

y 

MDSTO EX 

AF,AF ' 

5 GET SIGN OF RESULT 

AND 

80H 


JR 

Z, MDS20 

? POS RESULT 


! NEGATE BC 



LD 

HL , 0 



SBC 

H L., B C 



JR 

Z, MDSIO 

; RESULT WAS ZERO 


JP 

P, NINF 

; NEGATIVE OVERFLOW 

MDSIO 

L.D 

(IX),L 

; STORE RESULT 


LD 

(IX+1),H 



RET 



1 TEST 

FOR > 

15 BITS 


r 

MDS20 

BIT 

7 > B 



JR 

N Z , P1N F 

; POSITIVE OVERFLOW 


LD 

(IX),C 



LD 

(IX+1 ) ,B 



RET 

r 

; INF? CHECKS SIGN BIT IN A' 
? PINF PRODUCES POS OVERFLOW 
; NINF PRODUCES NEG OVERFLOW 

IDVO LD HL,ERROR 

SET 4,(HL) 

INF? EX AF,AF' 


FLAG ERROR 
GET SIGN BIT 







AND 

80H 



JR 

NZ,NINF 

; NEG SIGN 

PINF 

LD 

(IX),OFFH 

J RESULT IS 


L..D 

(IX+1),7FH 



JR 

FINF 


NINF 

LD 

(IX) ,0 

; RESULT IS 


LD 

(IX+1),80H 


FINF 

LD 

HL,ERROR 



SET 

6,(HL) 

? FLAG ERROR 


RET 

*H INTEGER DIVIDE 


lyytrrryryyyytrtrrtittrrryyyyryyyyyyyyyytfyytyrtyr 


IDIV INTEGER DIVIDE 


INPUT: IX PTS TO X 

IY PTS TO Y 

OUTPUT: IX PTS TO X/Y 

ERRORS: Y IS ZERO 


; ALGORITHM 

: NON-RES 

TORING DIVISION, 16 BY 1! 

: REG 

ISTERS 

: AC IS 

DIVIDEND, ACCUM EXT 

t 


DE IS 

DIVISOR 

t 


HL IS 

ACCUM, REMAINDER AT END 

7 

i »»»*«» 

t r r t r r y 

r t r r r r 

ttttlYttt 

7 r f r r r r t * 

yyyyyyyyyyyyyyyyyyyyyyyy 

IDI V 

CALL 

MDLOAD 

i LOAD ARCS 


LD 

AfD 



OR 

E - 



JR 

Z,IDVO 

; DIVISOR IS ZERO 

t 

LD 

A, B 

; X IN AC 

1 

5 16 BIT BY 1 

5 BIT UNSIGNED DIVIDE 

r 

LD 

HL , 0 

; CLEAR ACCUM 


LD 

B, 16 

; SET COUNTER 

IDV10 

RL 

C 

: 0 <= HL < DE 
; ROT ACC-RESULT LEFT 


RLA 




ADC 

HL ,HL 

? NO CY POSSIBLE 


SBC 

HL, DE 

? SUB DIVISOR 

IDV20 

CCF 


5 CALC RESULT BIT 


JR 

NC,IDV50 

; ACCUM WENT NEGATIVE 

IDV30 

DJNZ 

IDV10 

? COUNTER NOT 0 


JP 

IDV60 


IDV40 

RL 

C 

; -DE <= HL < 0 
? RO ACC-RESULT LEFT 


RLA 




ADC 

HL, H L 



AND 

A 

; TURN OFF CY 


ADC 

HL , DE 

; ADD DIVISOR 


JR 

C , I D V 3 0 

; ACCUM WENT POSITIVE 


JR 

Z,IDV20 

? ACCUM BECAME ZERO 

IDV50 

DJNZ 

ID V 4 0 

? COUNTER NOT 0 

IDV60 

RL 

c 

? SHIFT IN LAST RESULT 


RLA 




LD 

B, A 

: QUOTIENT IN BC 


JP 

MDSTO 

; STORE RESULT 


MDLOAD LOADS X INTO BC, Y INTO DE 

TAKES ABS OF EACH, STORES XOR OF SIGNS 





I 

; in 

A' 



1 

m 

? 

MDLOAD 

L..D 

b , (i x + 1 ) ; 

LOAD X 


m 


LO 

C, ( IX) 


m 



LD 

A, B 



m 


AND 

80H 


g 

I 


v.l R 

Z , M D L 1 0 ; 

X IS POSITIVE 

9 

< 

r 

; NEGATE X 



< 


7 

LD 

HL ,0 





SBC 

HI., BC 


M 

i 


LD 

BrH 


1 



LD 

C, L 



i 

7 

i LOAD 

Y 



€ 

.m 

7 

MDL..10 

LD 

D, ( IY+1 ) 


m 

1 


LD 

E, (IY) 


9 



X OR 

D J 

FORM SIGN OF RESULT 




EX 

AF,AF ' ; 

SAVE IT 

M 

• 


BIT 

7 , D 


1 



RET 

7 v 

7 

Y IS POSITIVE 


€ 

; NEGATE Y 



€ 


? 

LD 

HL , 0 


g 

I 


AND 

A 


1 



SBC 

HL,DE 



M 


EX 

DE, HL 


- m 

§ 


RET 



1 


*H IABS 

AND 

ISGN 



f 

* ****** 

? y ? 7 y 7 7 

r 

***** 
r y 7 7 7 

7777777777777 

777777777777777777777777777777777777777777 7 7777 

« 


5 IABS, 

INEG 




i 

F 




i 


y ? ? y ? ? y 

7 7 7 7 7 

7777777777777 

77777777777777777777777777777777777777777777777 



I ABS 

BIT 

7, < IX + i ) 





RET 

Z ; 

X IS POSITIVE 

m 

1 

INEG 

LD 

D , (IX +1) ; 

LOAD X 

1 



LD 

E. (IX) 





LD 

H L , 0 ? 

NEGATE X 


€ 


AND 

A 


i 



SBC 

HL , DE 



lS «-' 


JP 

PEfNINF ; 

NEG OVERFLOW 

M 

I 


LD 

(I x+1 ) , H ; 

STORE X 

f 



LD 

( IX ) , L. 



rf 


RET 





7 7 7 7 7 7 7 

* » t t i 

7 7 7 7 7 

************* 

tt?7fT7T7?I*T 

77777777777777777777777777777777777777777777777 

* 

< 

7 

; ISGN 

* 




< 


J > ? M M 

7 7 7 7 7 

*** ********** 
7777777777777 

77777777777777777 7 77777777777 7 77777777777 7 77777 



ISGN 

LD 

A , (IX +1 ) : 

GET SIGN 


i 


BIT 

7, A 


€ 



JR 

Z,IS G1o ; 

X IS NONNEG 


a 


LD 

(IX),OFFH i 

STORE -1 


49 


LD 

(IX+1),OFFH 


€ 



RET 





I SGI 0 

OR 

(IX) 



t 

f 





« 

f 




€ 


M 


i 


t 

« 

t 

t 

t 

t 

t 

< 

€ 

i 


RET Z 

U> (IX+1),1 
ID (IX) ,0 
RET 

*H ICR 


X IS 

STORE 


ZERO 

+ 1 


117111177117111111111171111171171111117171111771771111171111177117711711 


! ICP INTEGER COMPARE 

r 

: RETURNS CY SET IF X < Y 
; RETURNS Z SET IF X = Y 


7 7 7 7 

JfJJrrrrrj 

7 7 7 

♦ I t H t « 1 t » 1 V t Y » I 1 » 1 

7Tr7t7Tr777r77T7Tt 

LD 

Br(IX+1) 

T 

GET SIGN OF X 

LD 

A , B 



AND 

80 H 

7 

TEST SIGN r CLEAR 

JR 

NZ,ICP10 

7 

X IS NEG 

BIT 

7, (IY+1) 



RET 

NZ 

7 

Y IS NEG 

LD 

A > B 



CP 

(IY+1 ) 

7 

SIGNS ARE BOTH PO 

RET 

NZ 



LD 

A , (IX ) 



CP 

(IY) 



RET 




XOR 

(IY+1 ) 

7 

TEST SIGN, CLEAR 

RLA 


7 

SIGN BIT INTO CY 

RET 

c 

7 

SIGNS DIFFERENT 

LD 

A r B 



CP 

(IY+1 ) 

7 

BOTH SIGNS NEG 

RET 

NZ 



LD 

A,(IX) 



CP 

(IY) 



RET 





I 

I 

I 

« 

« 

i 

i 

< 

€ 

< 

i 


f 

f 

t 

€ 

C 

t 

f 


REFERENCES 


EXTERNAL MULT,DIVI,ERROR,OVER,UNDER,SAVE 
EXTERNAL SEED,SEED10,CRNDA,CRNDCO 

EXTERNAL SRD , LD, ML.DIG , MDTL 

EXTERNAL ZERO,ONE,NORML1,SLD,SUB,C1,MA D 

EXTERNAL ST04,TEMP4,RSUB,FIXA 

EXTERNAL LMARK,LREST,LBACK,RC2,LIN 

GLOBAL DTOB'BTOD 

GLOBAL RND,RANDMZ,INT,BRK,ABS,SGN,NEG 
GLOBAL I ADD , I SUB * IMUL » IDIV » I ABS » I SGN , ICF‘ > I NEG 


*H POWER 


1*777777777777777777777777777777777777 7 177777777777777777777 


J X TO THE Y POWER 

7 

; INPUT; IX PTS TO X 

; IY PTS TO Y 

; OUTPUT; IX PTS TO X+Y 

} ERRORS; NONE DIRECTLY 

; SPECIAL CASES; IF Y IS 0, RETURN i 

J ALGORITHM; IF Y IS AN EXACT INTEGER -32768 <= Y < 32768 

J THEN USE A SEQUENCE OF SQUARING AND 

l MULTIPLYING BY X TO PRODUCE RESULT 

! OTHERWISE » RESULT IS EXP(Y*LOG(X)) 


« 

€ 

t 

i 

i 

i 

i 

i 

i 


i 


i 



4 


i 

§ 


rtrrtrtrtryttTjrttffrtttttrrr 


DETERMINE IF Y IS AN INTEGER 


TTrttTtrrr y 


r r r t f 


4 


t 

« 

i 

< 

i 

€ 

f 

€ 

4 

4 

< 

< 

4 

I 


f 

C 

< 

t 


ROW 

PUSH 

IY 

; SAVE PTR TO Y 


PUSH 

IX 

f SAVE PTR TO X 


PUSH 

IY 

J PT TO Y 


POP 

IX 



CALL 

FIX A 

J CONVERT TO INTEGER 


POP 

IX 

i RESTORE IX 

* 

JR 

ZfPOWIO 

; Y WAS AN INTEGER 

; FOR 

NONINTEGER Y, TAKE EXP < Y*L.OG < X ) > 

r 

CALL 

LOG 

LOG OF X 


POP 

IY 

j RESTORE IY 


CALL 

MULT 

f MPY BY Y 


CALL 

EXP 

? EXP OF PRODUCT 


RET 



r 

; Y IS 

INTEGER, SET UP 

FOR LOOP> GET MAGNITUDE 

T 

F'OWIO 

XOR 

A 

5 CHECK FOR 0 


POP 

IY 

RESTORE IY 


CP 

< IY+7) 



JP 

Z, ONE 

5 X TO 0 POWER 


PUSH 

HL 

; SAVE Y 


BIT 

7 ; H 



JR 

Z.P0W15 

J Y IS POSITIVE 


EX 

DEfHL 

? NEGATE Y 


LD 

HL , 0 



AND 

A 



s BC 

HL , DE 



EX 

<SP) t HL 

SAVE NEGATED Y INSTEAD 


CALL 

REC 

? CALC 1/X 

P0W15 

CALL 

ST 04 

! SAVE COPY OF X 


POP 

HL 

r GET Y 


LD 

B, 16 

J SET COUNTER 

r 

; LOOK 
* 

FOR FIRST BIT 


T 

P0W20 

PUSH 

BC 

; SAVE COUNTER 


ADD 

HL, HL 

? SHIFT Y LEFT 1 


JR 

C , P0W40 

; 1 BIT FOUND 


POP 

BC 

; GET COUNTER 


D JN2 

P0W20 

; COUNTER NOT 0 


JP 

ONE 

J THIS JUMP SHOULD NOT HAPPEN 

T 

SQUARE FOR 

EACH ZERO 

BIT 

l SQUARE AND 

MULTIPLY 

FOR EACH l BIT 

y 

P0W30 

PUSH 

BC 

; SAVE COUNTER 


ADD 

HL 7 HL 

i SHIFT Y LEFT 1 


PUSH 

HL 

} SAVE Y 


PUSH 

AF 

J SAVE CY 


CALL 

ST 05 

J SAVE COPY OF ACCUM 


LD 

IY t TEMP5 



CALL 

MUL T 

? SQUARE X 


POP 

AF 

} GET CY 


I 

4 

I 

4 

i 

4 

4 

4 

4 

4 

4 

i 

4 

4 

4 

4 

4 

4 

4 


C 


i 



P0W35 

P0W40 


JR NC, ROW35 
CALL MULT4 
POP HL 
POP BC 
DJNZ P0W30 
RET 


; WAS A 0 BIT 
; MPY BY X 
; GET Y 
? GET COUNTER 
; COUNTER NOT 0 


; IPOW POWERS FOR INTEGER ARGUMENTS 
! USES REAL POW ROUTINE BY FLOATING ARCS 

! IN TEMPS, TEMP6, (UNUSED BY POW WHEN Y 

; IS INTEGER), CALLING POW, AND FIXING RESULT 


IPOW 


PUSH 

IX 

* 

7 

SAVE ORIGINAL PTR 

LO 

L, ( IY) 

t 

GET Y 

L D 

H,(IY+1) 



LD 

IX,TEMP6 



CALL 

FLOAT 1 

y 

FLOAT Y IN TEMF’6 

POP 

IX 

t 

7 

GET ORIGINAL PTR 

PUSH 

IX 



L.D 

L , (IX) 

t 

7 

GET X 

LD 

H,(IX + 1 ) 



LD 

IX,TEMPS 



CALL 

FLOAT 1 

7 

FLOAT X IN TEMPS 

LD 

IY , TEMP6 

T 

PT TO Y 

CALL 

POW 

7 

CALL REAL ROUTINE 

CALL 

FIX 

f 

FIX RESULT 

POP 

IX 

7 

GET ORIGINAL PTR 

LO 

( IX) ,L 

/ 

7 

STORE RESULT 

LD 

(IX+1),H 



RET 




:E ROOT 



t r r r 


J J r 

TTTrTTT77T7?7T777T 


77777777777777777777 


J SQUARE ROOT FUNCTION 


INPUT: IX PTS TO X 
OUTPUT: IX PTS TO SQR(X) 

ERRORS: X IS NEGATIVE, RETURNS SQR<ABS(X>) 
SPECIAL CASES: X IS 0 

ALGORITHM: INPUT IS SCALED TO ,1 <- X < 1 

BY REMOVING THE EXPONENT. HALF THE 
EXPONENT IS ADDED TO THE RESULT EXPONENT 
AND THE RESULT IS MULTIPLIED BY SQR(IO) 
IF THE ORIGINAL EXP WAS ODD. 

AN INITIAL APPROXIMATION IS GIVEN BY 
A 1ST ORDER POLYNOMIAL, AND IMPROVED TO 
REQUIRED ACCURACY IN 4 ITERATIONS OF 
NEWTON’S METHOD. 


! CHECK VALID ARGUMENTS 


SQR 


L..D A,(IX+7) 
OR A 
JP Z,ZERO 
BIT 7,(IX) 
CALL NZ,SQERR 


ZERO INPUT 
NEGATIVE INPUT 


: SCALE BY STRIPPING EXPONENT 



PUSH AF } 

LD (IX + 7>,80H‘ ? 


; CALC INITIAL APPROX 


CALL 

ST 04 

? 

PUSH 

IX 


LD 

IX,TEMP4 


CALL 

HALF 

7 

POP 

IX 


LD 

lYrF'SQOO 


CALL 

ADD 

T 

t 

; PERFORM 4 

ITERATIONS 


T 

LD 

B, 4 


S Q10 PUSH 

BC 

7 

CALL 

S T 0 5 

T 

LD 

HL,TEMP4 


CALL 

LD 

T 

LD 

IY,TEMPS 


CALL 

DIVI 

1 

PUSH 

IX 


LD 

IX,TEMPS 


CALL 

HALF 

’I 

7 

POP 

IX 


CALL 

ADD 

? 

POP 

BC 


DJNZ 

SQ10 

i 

T 

! CALC EXPONENT, SCALE 

RE 

T 

POP 

AF 

7 

AND 

A 


RRA 


% 

T 

RR 

B 

r 

ADD 

A, 40H 

7 

LD 

< IX + 7) , A 

7 

BIT 

7» B 

; 

RET 

Z 

r 

LD 

IY,SQRiO 


CALL 

MULT 

% 

7 

RET 



T 

i SQUARE ROOT OF NEGATIVE 

7 

SQERR LD 

HL , ERROR 


SET 

3 , ( Hl.) 

r 

RES 

7, (IX) 


RET 



*H EXP 




tfftttttrfr 

r ? 


SAVE ORIGINAL EXP 
SET ZERO EXP 


SAVE X 


CALC X/2 IN TEMP4 


INIT APPROX, 

SAVE COUNTER 

SAVE LAST APPROX IN TEMP5 
GET X/2 

DIVIDE BY LAST APPROX 

DIVIDE LAST APPROX BY 2 
CALC NEW APPROX 
COUNTER NOT 0 
5ULT 

GET ORIGINAL EXP 

DIVIDE BY TWO 
SAVE CY 

CORRECT THE BIAS 
STORE IT 

TEST ODD OR EVEN 
ORIGINAL EXP EVEN 

MPY BY SQR(iO) 

NUMBER 


FLAG ERROR 
HAKE POSITIVE 


J EXPONENTIATION ROUTINE 

? INPUT: IX PTS TO X 

t OUTPUT: IX PTS TO E+X 

; ERRORS: OVER/UNDERFLOW IN HANDLING EXPONENT PART 
? ALGORITHM: X : ~X*2*L.0G BASE 10 OF E 



SAVE J = I NT(X) 

K:-FRACTION PART(X ) / 2 
X;=10+X, WHERE 0 <= X < .5 
IF SAVE IS ODD THEN X*=X*SQR(10) 

X;=X*10+INT(SAVE/2) 

THE APPROXIMATION TO 10+X IS GIVEN BY ! 

( Q ( X + 2 ) + X*P(X+2>> / < Q ( X+2 ) - X*P(X+2>) 
WHERE Q IS 3RD ORDER POLYNOMIAL, P IS 2ND ORDER. 
REFERENCE IS "COMPUTER APPROXIMATIONS", HART, JOHN F. 
ET ALf FUNCTION NUMBER 1444 


5 CONVERT E + X TO 10 + X 

r 

EXP LD IY,L0GE2 ; MULT BY 2*L0G10(E) 

CALL MIJLT 

r 

; GET 8 BIT 2'S COMPL INTEGER PART IN A 

f 

CALL BRK 5 BREAK INTEGER AND FRACTION PARTS 

JP C,EXERR1 i INTEGER PART > 16 BITS 

L. D A , H 

AND 1 ? SIGN BIT IN AO 

ADD A,H ; 0 IFF H WAS FF OR 00 

JP NZ,EXERR2 ; INTEGER > 9 BITS 

PUSH HI. ? SAVE IT 

t 

? 

! SCALE FRACTION PART TO 0 <= X < .5 

f 

XOR A 
CP (IX+7) 

JR 2,EX10 ; FRACTION IS ZERO 

CALL HALF 

JP CfUNDER ? UNDERFLOW 

T 

? EVALUATE RATIONAL APPROX. TO 10+X 


L A L L 

V V o 

f \ f\ *L. 


STORE X, X+2 IN TEMP4, 5 

LD 

IY,PEX02 



LD 

B, 2 

* 

r 

2ND ORDER POLY 

C ALL- 

POLY 

T 

EVALUATE P(X+2> 

CALL 

MULT 4 

» 

y 

CALC. X*P(X + 2 > 

CALL 

ST 04 

r 

SAVE X*P IN TEMP4 

LD 

HL , TEMPS 



CALL 

LD 

t 

GET X+2 

LD 

IY,QEX02 



LD 

Bf2 

y 

3RD ORDER POLY(l) 

CALL 

P0LY1 

T 

EVALUATE Q(X+2) 

CALL 

ST 05 

y 

SAVE Q IN TEMPS 

LD 

IY f TEMF‘4 



CALL 

SUB 



CALL 

ST 06 

y 

SAVE Q-XP IN TEMP6 

LD 

HL,TEMPS 



CALL 

LD 

y 

GET Q 

LD 

IY,TEMP4 



CALL 

ADD 

y 

CALC. Q+XP 

LD 

IY,TEMP6 



CALL 

DIVI 

» 

y 

FINAL RESULT 


POSSIBLE SCALE BY SQR(IO) 


POP 

HL. 

; GET 9 BIT INTEGER 

SR A 

H 

i DIVIDE BY 2 

RR 

L. 


PUSH 

HI. 

J SAVE 8 BIT INTEGER 

LD 

IY,SQR10 


CALL 

C,MULT 

; INTEGER WAS ODD 


T 

; ADD INTEGER PART TO RESULT EXPONENT 


POP BC 
LO A, (IX + 7 ) 

SUB 80H 
ADD A , C 
JP PE,EKERR3 
ADD A 7 80H 
JP Z ,UNDER 
LD (IX+7),A 

RET 


PUT INTEGER IN C 
GET EXP 
REMOVE BIAS 
ADD INTEGER PART 
OVER/IJNDERFLOW 
ADD BIAS 
UNDERFLOW 
STORE EXP 


; EXP ERRORS 


EXERR1 

EX 

AF,AF ' 


LD 

H, A 

E XERR2 

BIT 

7, H 


JP 

NZ,UNDER 


JP 

OVER 

EXERR3 

JP 

P,UNDER 


JP 

OVER 


*H LOG 


GET ORIGINAL SIGN 

- SIGN 

+ SIGN AFTER OVERFLOW 

Vl%t******t%t***lr , r*»* 

777777777777777777777 


LOG ROUTINE 


INPUT; IX PTS TO X 

OUTPUT t IX PTS TO NATURAL LOG OF X 

ERRORS; X IS NEGATIVE, RETURN LOG(ABS(X)> 

X IS ZERO , RETURN NEGATIVE OVERFLOW 
SPECIAL CASES; X IS i, RETURN ZERO 
ALGORITHM; SAVE?=INT(X) 

X;= SQR(10)*FRACTI ON-PART(X > 

X;=LOG BASE E (X), WHERE 1/SQRdO) < X < 
X;=X - LOG BASE E < SQR(10)> + SAVE*LOG BASE 

LOG IS APPROXIMATED BY ; 

Z*P(Z+2) / QCZ+2) WHERE Z = <X~1) 
AND P IS A 3RD ORDER POLYNOMIAL, AND Q IS 
REFERENCE FUNCTION IS 2686. 


; CHECK FOR VALID ARGUMENTS 


LOG 


BIT 7,(IX) 

CALL NZ,LGERR1 ; NEGATIVE ARG 
LD IY,C1 

CALL CP 
JP Z,ZERO 


SQR(10 ) 
E (10) 


/ (X+l) 
4TH ORDER. 


LOG (1 ) 


LD 

A , ( IX + 7) 

i GET EXPONENT 

OR 

A 


UP¬ 

Z »LGERR2 

: ZERO ARG 

TAKE LOG OF 

EXPONENT 

PART 

SUB 

BOH 

? REMOVE BIAS 

L.D 

L r A 

? PUT IN HL 

LD 

H , 0 

5 SIGN EXTEND 

JR 

NCfLOIO 

? + SIGN 

LD 

H.OFFH 


010 PUSH 

HL 

i SAVE INTEGER 

SCALE TO 1/ 

SQR(10) < 

X < SQR(10) 

LD 

(IX+7),BOH 

; SCALE TO 0 < X < 1 

LD 

IY t SQR10 


CALL 

MULT 


APPROXIMATION FOR LOGARITHM 

CALL 

S T 0 4 

; SAVE X IN TEMP4 

PUSH 

IX 

i SAVE PTR 

LD 

IX , TEMP4 


LD 

IY , Cl 


CALL 

ADD 

J X+l IN TEMP4 

POP 

IX 


L. D 

IYfCI 


CALL 

SUB 

J X-l IN (IX) 

LD 

IY r TEMP4 


CALL 

DIVI 


CALL 

XX 2 

; STORE Z.Z+2 IN TEMP4, 

LD 

IY ,Ql.003 


LD 

B, 3 


CALL 

POL. Y1 

; EVALUATE Q(Z+2) 

CALL 

ST 06 

; Q IN TEMP6 

LD 

HL »TEMF'5 


CALL 

LD 

i LOAD Z+2 

L D 

IY,PL003 


LD 

B , 3 


CALL 

POLY 

; EVALUATE P(Z^2) 

CALL 

MIJLT4 

? CALC. Z*P ( Z t- 2 ) 

LD 

IY,TEMP6 


CALL 

DIVI 

J CALC ZP/Q 

L D 

IYrLNSQlO 


CALL 

SUB 

SUB LOGE ( SQR (10)) 

CALL 

ST 04 

? SAVE IN TEMP4 

CONVERT 16 

BIT EXPONENT PART TO FLOATING PT 

POP 

HL 

J GET INTEGER 

CALL 

FLOAT1 

? FLOAT IT 

LD 

IY , L. N10 


CALL 

MULT 

? MPY BY LOG BASE E(10) 

LD 

IY , T E M P 4 


CALL 

ADD 

) FINAL RESULT 


RET 

r 

1 ERRORS 


LGERRA 

LD 

HL,ERROR 



SET 

2 , ( HI. ) 

? FLAG LOG OF NEG 


RES 

RET 

7,(HL) 

} TAKE ABS VALUE 

LGERR2 

SET 

7,(IX) 

i SET NEGATIVE SIGN 

JP 

*H TANGENT 

OVER 

i LOG OF ZERO 

r y y y y y y 

r f t r j 

yyyyyyyyyy 

rtrryyiy?7?TtTrrTtryt 


; TANGENT ROUTINE 


INPUT? IK PTS.TO X (IN RADIANS) 
OUTPUT? IK PTS TO TAN(X) 

ERRORS? NONE DIRECTLY 
ALGORITHM? X s= X*4/PI 


7 

SAVE ?= 

INT(X ) LAND 

3 C 

SAVE 

LOW 2 BITS 

T 

X ?= FRACTION-PARK 

X) 



T 

IF SAVE 

IS 0,2 THEN 

X ? = 

1-X 


r 

X := T A N ( X * PI / 4 )» 

WHERE 

0 < = 

* X <> PI/4 

7 

IF SAVE 

IS 0,3 THEN 

V * «« 

t\ * — 

1/X 


r 

IF SAVE 

IS 0,1 THEN 

\/ t — 

s\ * — 

_x 



THE APPROXIMATION TO TAN(X*PI/4 ) IS GIVEN BY? 
X*P(X+2) / 0(X+2) 

WHERE P IS ORDER 3 POLYNOMIAL, 0 IS ORDER 4 
REFERENCE FUNCTION IS 4286 


r 

SCALE BY 4/PI 


TAN 

I..D 

IY,FOVRPI 




CALL 

MULT 




CALL 

BRK 

7 

GET INTEGER,FRACTION 


LD 

A r L 




PUSH 

AF 

* 

y 

SAVE FLAG BITS 


AND 

1 

y 

TEST FLAG 


JR 

Z , T A10 

y 

FLAG WAS 0 OR 2 


LD 

IY , C1 




CALL 

RSUB 

y 

CALC. 1-FRACTION PART 

y 

; TAN 

APPROX 




f 

TA10 

CALL 

XX2 

y 

SAVE X,X+2 IN TEMP4,5 


LD 

IY , 0 T A 0 3 




LD 

B , 3 




CALL 

POLY 1 

y 

EVALUATE Q(X+2) 


CALL 

S T 0 6 

y 

SAVE 0 


LD 

HL, TEMPS 




CALL 

LD 

y 

GET X-t-2 


LD 

IY , F'TA03 




L.D 

B, 3 




CALL 

POLY 

y 

EVALUATE P(X+2) 


CALL 

MULT 4 

y 

CALC. X*P(X+2) 


POP 

AF 

y 

GET 2 BIT FLAG 


PUSH 

AF 




AND 

3 




JP 

PE,TA20 

I 

FLAG WAS 0 OR 3 


CALL 

ST04 

y 

SAVE X*P 


LD 

HL , TEMP6 






CALL LD i GET Q 

L.D IY,TEMF‘4 ? SET FOR Q/XP 

JP TA30 

TA20 LD IY » TEMP6 J SET FOR XP/Q 

TA3Q CALL 01VI 


J CALC. SIGN OF RESULT 

f 

r 

POP AF j GET FLAG 

AND 2 

RET 7. j ON 0 OR 1 

SET 7, (IX) 

RET 

*H COSINE,SINE 

777777777777777777777777777777777777777777777777777777777777777777777777 


; COSINE ROUTINE 

J EFFECTIVELY ADDS PI/2 AND CALLS SINE 

r 

777777777777777777777777777 7 77777777777777777777777777777777777777777777 

COS LD IYfFOVRPI 

CALL MULT ; SCALE BY 4/PI 

CALL BRK t GET INTEGER, FRACTION 

L D A , L 

ADD A,2 ? ADD EFFECTIVE PI/2 

JP SI10 

777777777777777777777777777777777777777777777777777777777777777777777777 

* 

7 

} SINE ROUTINE 


? 

INPUT ; 

IX 

PTS TO X (IN RADIANS) 


7 

OUTPUT : 

IX 

PTS TO SIN(X) 



7 

ERRORS: 

NONE 




T 

ALGORITHM: 

X ?= X * 4/P I 



J 

SAVE 

; =s 

INT(X) LAND 7 



7 

X : = 

FRACTION-PART(X) 



7 

IF S 

AVE 

IS 1,3,5,7 THEN 

V * — 
f\ * 

i-X 

7 

IF S 

AVE 

IS 0,3,4,7 THEN 

v * — 

t\ * — 

SIN ( X * PI / 4 ) 

7 

IF S 

AVE 

IS 1,2,5,6 THEN 

X ? = 

COS(X*PI/4 ) 

? 

IF S 

AVE 

IS 4,5,6,7 THEN 

X : = 

-X 


THE SINE APPROX FOR 0 O X <= PI/4 IS? 
X*P(X+2) / Q < X t-2 ) 

WHERE P IS A 3RD ORDER POLYNOMIAL, Q IS 3RD 
REFERENCE FUNCTION IS 3063. 


ORDER. 


THE COSINE APPROX IS? 

P(X+2) / Q(X+2) 

WHERE P IS A 3RD ORDER POLY, 0 IS 3RD ORDER, 
REFERENCE FUNCTION IS 3843. 


7777777 7 777777777777777777777777777777777 7 777777777777777777777777777777 

SIN LD IYfFOVRPI 

CALL. MULT ; SCALE BY 4/PI 

CALL BRK ? GET INTEGER AND FRACTION 

L D A, L 

S110 PUSH AF ? SAVE 3 BIT FLAG 

AND 1 

JR 2,SI20 ; FLAG IS 0,2,4,6 

LD IY,Cl 


I 


i 


I 

« 

« 

t 

« 

< 

« 

€ 

« 

« 

1 

I 

S 

I 

€ 

< 

< 



CALL 

RSUB 

SI 20 

POP 

AF 


PUSH 

AF 


AND 

3 

* 

J P 

F'O, SI30 

J SIN 

APPROX 


t 

CALL 

XX2 


LD 

IY,QSI02 


LD 

B, 2 


CALL 

POL Y1 


CALL 

S T 0 6 


LD 

HL,TEMPS 


CALL 

LD 


LD 

IY,PSI03 


LD 

B, 3 


CALL 

POLY 


CALL 

MULT4 


L. D 

IY , T E M P 6 


CALL 

DIVI 


JP 

SI 40 

t 

; COS 

APPROX 


? 

SI 30 

CALL 

XX 2 


LD 

IY , QC002 


L D 

B, 2 


CALL 

P0LY1 


CALL 

ST06 


LD 

HL,TEMPS 


CALL 

LD 


LD 

IY, PC003 


LD 

B , 3 


l- A L. L 

POLY 


LD 

IY,TEMP6 


CALL 

DIVI 

f 

J CALC 

SIGN < 

OF RESULT 

SI 40 

POP 

AF 


AND 

4 


RET 

Z 


SET 

RET 

7, < IX ) 

*H ARC 

TANGENT 

f r r ? t t 


M M I M M 


; AR C T AN GE N T R0 U TIN E 


CALC. i-X 
GET FLAG 

1 , 2 , 5 , 6 

STORE X , X+2 IN TEMP4,5 

EVALUATE 0<X+2) 

SAVE Q 

GET X+2 

EVALUATE P<X+2> 

CALC. X* P (X + 2) 

CALC. XP/8 

STORE X , X +2 IN TEMP4,5 

EVALUATE Q(X + 2) 

GET X+2 

EVALUATE P(X+2) 

CALC. P/0 

GET FLAG 

0 , 1 , 2,3 


r ? 7 ? ? ? 


t r r i i 


i 

l 

€ 

i 

€ 

i 

i 

f 

< 

i 

i 

« 

« 

i 

« 

< 


m i INPUT: IX PTS TO X 

« ? OUTPUT: IX PTS TO ATN(X) (IN RADIANS) 

i ERRORS: NONE 
; ALGORITHM: SAVE := SGN(X) 

• J FLAG := 0 

; IF TAN ( 3*F'I / 6 ) <= X THEN X:=l/X, FLAG : = 1 

5 IF TANCPI/8) <= X THEN X:=(X-l)/(X+l), FLAG:=2 

• J X := ATN(X) , WHERE 0 <= X <= TAN(PI/8) 

? IF FLAG = 1 THEN X := X + PI/4 

^ ? IF FLAG = 2 THEN X := PI/2 - X 


i 

< 

i 

i 


f 


f 


! APPROXIMATION FOR ATN(X) IS? 

i X # P ( X + 2 ) / Q(X+2) 

} WHERE P IS 3RD ORDER POLYNOMIAL * 9 IS 4TH ORDER 

? REFERENCE FUNCTION IS 5075 

1 

tirtT7fiTrryttr?TTTTrrifrtttrttTtftrtfrtttrTrTrftrtttr7 

* 

} TAKE ABS VALUEr DETERMINE RANGE 


ATN 

L.D A, (IX) 

AMD 80H 

V 

T 

GET SIGN 


PUSH AF 

r 

SAVE SIGN 


RES 7,(IX) 

r 

MAKE POSITIVE 


LD IY,TAPIS 

CALL C P 

r 

COMPARE TO 1ST ! 


L.D B , 0 

JR C > AT20 

7 

< T A N ( PI / 8 ) 


LD IY i TA3PI8 

CALL CP 

1 

COMPARE TO 2ND 


JR C,AT1G 

7 

< TAN(3*PI/8) 

r 

TAN ( 3*F'I /8 ) <= X 



; TAKE 1/X NOW 5 SUB 

FROM PI72 AT END 


CALL REC 

7 

RECIPROCAL 


LD B 1 2 

T 

SET FLAG TO 2 


JP AT20 



7 

J TAN 

(F'I/8) <= X < TAN(3*PI/8) 

X T C 

ALC (X-l)/(X+i) NOW, 

ADD PI/4 AT END 

AT10 

CALL ST04 

PUSH IX 

L.D IX » TEMF‘4 

LD IY , C1 

7 

SAVE X IN TEMP4 


CALL ADD 

T 

X+l IN TEMP4 


POP IX 

LD IY , C1 

CALL SUB 

LD IY 7 TEMP4 

CALL DIVI 

7 

X-l IN ACCUM 


L D B 7 i 

7 

SET FLAG TO 1 

7 

; APPROX FOR ARCTAN 



7 

AT20 

PUSH BC 

7 

SAVE FLAG 


CALL XX2 

LD IY , QAT03 

T 

SAVE X 7 X * 2 


LD B 7 3 

CALL POLY! 

7 

CALC. Q(X+2) 


CALL ST06 

? 

SAVE Q 


LD HL,TEMPS 

CALL. LD 

LD IY,PAT03 

r 

GET X+2 


LD B, 3 

CALL POLY 

7 

CALC P(X+2) 


CALL MULT4 

LD IY , TEMF'6 

r 

X*P 





CALL 

DIVI 

i X*P/Q 



POP 

BC 

? GET FLAG 



BIT 

0 t B 




JR 

Z»AT30 

? FLAG IS 0 

OR 2 

; SECOND INTERVAL (FLAG 

IS 1) 



LD 

I Y r P I 4 




CALL 

ADD 

I ADD PI/4 



JP 

A T 4 0 



AT30 

BIT 

1 »B 




JR 

Z,AT40 

i FLAG IS 0 


; THIRD 

INTERVAL (FLAG 

IS 2) 



LD 

IYrPI2 




CALL 

RSUB 

; SUB FROM 

PI/2 

T 

} GIVE 

t 

RESULT 

SIGN OF ARGUMENT 


T 

AT 40 

POP 

AF 

? GET SIGN 



RET 

~r 

A. 

J SIGN WAS 

POSITIVE 


SET 

7,(IX) 

i MAKE NEGATIVE 


RET 





*H POLY 


} 

; POLYNOMIAL EVALUATOR 

? INPUT; IK PTS TO ACCUM 
f IY PTS TO HIGHEST COEFFICIENT 

; TEMPS CONTAINS X 

? B CONTAINS ORDER (ORDER-1 FOR POL.Y1 ) 

; output; ix pts to result 

; (POLYl ASSUMES FIRST COEFF, IS 1) 


POLY 

PUSH 

BC 

f 

SAVE COUNTER 


PUSH 

IY 

r 

SAVE COEFF. PTR 


JP 

POLY.10 



POLYl 

PUSH 

BC 

* 

T 

SAVE COUNTER 


PUSH 

IY 

J 

SAVE PTR 


CALL 

ADD 




LD 

IY f TEMPS 

7 

PT TO X 

POLY 10 

CALL 

MULT 




POP 

IY 




POP 

BC 




LD 

DE r - 8 

y 

ADJUST PTR 


ADD 

IY j DE 




D JNZ 

POLYl 

? 

COUNTER NOT 0 


CALL 

ADD 

* 

T 

FINAL ADD 


RET 




f 

POLYNOMIAL 

UTILITIES 



XX2 

CALL 

S T 0 4 

* 

STORE X IN TEMP4 


CALL 

MIJLT4 

r 

CALC. X+2 


CALL 

RET 

ST 05 



MULT 4 

LD 

IY 7 TEMP4 




CALL 

MULT 




RET 




*H CONS 

TANTS 




SQR10 

DEFW 

1603H 

V 

? 

3.1622 77660 168 


DEFW 

7 7 2 2 H 




DEFW 

01 6 6 H 






DEFW 

81.68 H 





LNSQiO 

DEFW 

1501H 

r 

1 . 1512 

92546 

497 


DEFW 

9212H 






DEFW 

6454H 






DEFW 

8197H 





LNIO 

DEFW 

3QQ2H 


6. « v? U J 

85092 

994 


DEFW 

8525H 






DEFW 

2909H 






DEFW 

8194H 





L0GE2 

DEFW 

6808H 

1 

7 

.86858 

89638 

065 


DEFW 

8958H 






DEFW 

8063H 






DEFW 

8065 H 





FOVRPI 

DEFW 

2701H 

7 

1.2732 

39544 

735 


DEFW 

3932H 






DEFW 

4754H 






DEFW 

813 5 H 





TAPIS 

DEFW 

i 4 0 4 H 

T 

.41421 

35623 

731 


DEFW 

352 1H 






DEFW 

3762H 






DEFW 

8031H 





T A3PI8 

DEFW 

4 i 0 2 H 

7 

2.4142 

13562 

373 


DEFW 

1342H 






DEFW 

2 3 5 6 H 






DEFW 

8173H 





PI4 

DEFW 

850 7 H 

t 

.78539 

81633 

974 


DEFW 

8139H 






DEFW 

3963H 






DEFW 

8074 H 





P12 

DEFW 

5701H 

7 

1.5707 

96326 

795 


DEFW 

9 6 0 7 H 






DEFW 

6732H 






DEFW 

819 5 H 





PSQOO 

DEFW 

710 i H 

7 

. 17157 

28752 

538 


DEFW 

2 8 5 7 H 






DEFW 

257 5 H 






DEFW 

8038 H 





PEXOO 

DEFW 

17 OSH 

% 

5178.0 

91 991 

516 


DEFW 

9180 H 






DEFW 

1599H 






DEFW 

8 416 H 





PEXOi 

DEFW 

31 OSH 

7 

831.40 

67212 

937 


DEFW 

6 7 4 0 H 






DEFW 

2921H 






DEFW 

8337 H 





PEX02 

DEFW 

8301H 

7 

18.312 

36015 

928 


DEFW 

3 612 H 






DEFW 

5901H 






DEFW 

8228 H 





QEXOO 

DEFW 

4 9 0 4 H 

7 

4497.6 

33557 

406 


DEFW 

3 3 7 6 H 






DEFW 

7455H 






DEFW 

8406 H 





QEXOl 

DEFW 

7002H 

7 

2709.3 

16940 

852 


DEFW 

1693H 






DEFW 

089 4 H 






DEFW 

8452H 





8EXQ2 

DEFW 

5901H 

? 

159 . 37 

41523 

603 


DEFW 

4137H 






DEFW 

3652H 










OEFW 

8303 H 

! QEX03 

IS 1 


PL000 

DEFW 

6502H 


OEFW 

2452H 


DEFW 

8590H 


DEFW 

8 316 H 

PL 001 

DEFW 

2 9 8 4 H 


DEFW 

3 4 4 8 H 


DEFW 

8682H 


DEFW 

8358 H 

PL 00 2 

DEFW 

9701H 


DEFW 

4 6 6 4 H 


DEFW 

7029H 


DEFW 

8335 H 

PL 00 3 

DEFW 

2782 H 


DEFW 

7 6 6 4 H 


DEFW 

7115H 


DEFW 

8215H 

Q LOGO 

DEFW 

3201 H 


DEFW 

1276H 


DEFW 

424 5H 


DEFW 

8359 H 

6 LOO 1 

DEFW 

5682H 


DEFW 

5499H 


DEFW 

9289H 


DEFW 

8300 H 

Q| 002 

DEFW 

5801H 


DEFW 

1860H 


DEFW 

2796H 


DEFW 

8327 H 

a LOO3 

DEFW 

1483H 


DEFW 

4 816 H 


DEFW 

824 4H 


DEFW 

8282H 

; QL004 

IS 1 


PTAOO 

DEFW 

0801H 


DEFW 

6088H 


DEFW 

3704H 


DEFW 

8828 H 

PTA01 

DEFW 

9588H 


DEFW 

8488H 


DEFW 

0640H 


DEFW 

8677 H 

PTA02 

DEFW 

4101H 


DEFW 

8589H 


DEFW 

5242H 


DEFW 

8576 H 

PI AO 3 

DEFW 

568 4 H 


DEFW 

3149H 


DEFW 

389 4 H 


DEFW 

8267 H 

UIA 0 0 

DEFW 

3801H 


DEFW 

7963H 


DEFW 

6366H 


DEFW 

8857H 

QTA01 

DEFW 

9983H 


DEFW 

0913H 


DEFW 

8051H 


DEFW 

8735H 

QTA02 

DEFW 

3501H 


; 265.52 24908 516 

J -429.48 34828 658 

: 197.64 46297 035 

; -22.764 76157 115 

? 132.76 12454 259 

; -258.99 54899 200 

J 158.60 18962 727 

? -31.416 48448 282 

10888 600.43 728 

? -89588 8.4400 677 

; 14189. 85425 276 

; -45649 31943 867 

; 13863 796.66 357 

? -39913 09.518 035 

; 13538 2.7128 051 





DEFW 

2 7 3 8 H 




DEFW 

8012H 




DEFW 

8651 H 



QTA03 

DEFW 

0181H 

7 

-1014.6 56190 253 


DEFW 

5 6 4 6 H 




DEFW 

0219H 




DEFW 

8453 H 



5 QTA04- 

IS 1 




PSIOO 

DEFW 

0602H 

i 

7 

20664 34.333 700 


DEFW 

3 4 6 4 H 




DEFW 

3733H 




DEFW 

8700 H 



PS 101 

DEFW 

818.1 H 

t 

-18160 3.9879 741 


DEFW 

3960H 




DEFW 

9787H 




DEFW 

8641 H 



P S102 

DEFW 

5903H 

7 

3599.9 30694 964 


DEFW 

3 0 9 9 H 




DEFW 

4969 H 




DEFW 

8464 H 



PS 103 

DEFW 

0182 H 

; 

-20.107 48329 459 


DEFW 

4807H 




DEFW 

9 4 3 2 H 




DEFW 

8259 H 



Cl SI 00 

DEFW 

6302H 

7 

26310 65.910 265 


DEFW 

651 OH 




DEFW 

0291H 




DEFW 

8765 H 



0 SI 01 

DEFW 

9 2 0 3 H 

? 

39270. 24277 465 


DEFW 

2470H 




DEFW 

7427H 




DEFW 

8565 H 



0 S1 0 2 

DEFW 

7802H 

7 

278.11 91948 108 


DEFW 

9111H 




DEFW 

8194H 




DEFW 

8308 H 



1 QSI03 

is i 




PC 000 

DEFW 

2901H 

r 

7 

12905 39.465 904 


DEFW 

3905H 




DEFW 

5946H 




DEFW 

8704 H 



P COO 1 

DEFW 

7 4 8 3 H 

7 

-37456 7.0391 572 


DEFW 

7 0 5 6 H 




DEFW 

15 3 9 H 




DEFW 

8672 H 



P COO 2 

DEFW 

3401 H 

7 

13432. 30098 654 


DEFW 

3 0 3 2 H 




DEFW 

860 9 H 




DEFW 

8554 H 



PC 00 3 

DEFW 

1281H 

\ 

7 

-112.31 45082 334 


DEFW 

4531H 




DEFW 

230 8 H 




DEFW 

8334 H 



& Co oo 

DEFW 

2901H 

? 

12905 39.465 904 


DEFW 

3905H 




DEFW 

5946H 




DEFW 

8704 H 



8 COOl 

DEFW 

3402H 


23467. 77310 725 


DEFW 

7 7 6 7 H 




DEFW 

0731H 




Q COO 2 

DEFW 

Q902H 

: 209.69 51819 673 



DEFW 

5169H 




DEFW 

968.IH 




DEFW 

8373 H 



; QC003 

IS 1 




P A TOO 

DEFW 

i 602H 

J 216.06 23078 972 



DEFW 

2306H 




DEFW 

8907H 




DEFW 

8372 H 



PAT 01 

DEFW 

2203H 

J 322.66 20700 133 



DEFW 

2066H 




DEFW 

0170H 




DEFW 

8333 H 



PAT02 

DEFW 

3201H 

5 132.70 23981 640 



DEFW 

2370H 




DEFW 

1698H 




DEFW 

8 3 4 0 H 



P AT03 

DEFW 

2801H 

5 12.888 38303 416 



DEFW 

3888H 




DEFW 

3430H 




DEFW 

8216H 



QATOG 

DEFW 

1602H 

J 216.06 23078 972 



DEFW 

DEFW 

2306H 

8907H 




DEFW 

8372H 



QATOl 

DEFW 

9403H 

; 394.68 28393 123 



DEFW 

2868H 




DEFW 

3139H 




DEFW 

8323H 



QAT02 

DEFW 

2102H 

J 221.05 08830 284 



DEFW 

0805H 




DEFW 

0283H 




DEFW 

8384H 



QAT03 

DEFW 

8503H 

? 38.501 48650 835 



DEFW 

480 IH 




DEFW 

0865H 




DEFW 

8235H 



J QAT04 

IS i 




7 

! REFERENCES 




T 

GLOBAL 

ROW , I F'OW , SQR , EXP , LOG , TAN , COS 

,SIN,ATN 


EXTERNAL- 

ERROR,OVER,UNDER,ZERO,ONE,Cl 



EXTERNAL 

ADD,SUB,RSUB,MULT,DIVI,REC 



EXTERNAL 

HALF 



EXTERNAL- 

ST 04,ST 05,ST 06 , TEMP4 , TEMPS,TEMP6 


EXTERNAL 

L.D , CP , FIX, FIXA , FLOAT , FLOAT 1 , 

BRK 


EXTERNAL 

SRBr NORMLZ 




! MISCELLANEOUS ROUTINES REQUIRED AS EXTERNALS FOR MATH PACKAGE 


* * ¥ Y * * 

F t t t f t 

LMARK 

L REST 
LBACK 
R C 2 


LD (POINT),HL 
RET 


RET 

RET 

RET 


77777777777777777777 


} STORE PTR FOR DTOB 


7777777777777 


l ENTRY FOR UNSUCCESSFUL DTOB (NO DIGITS FOUND) 
f DUMMY ROUTINE (BACK UP PTR FOR BASIC) 
t NORMAL RETURN» IF SOME KIND OF NUMBER FOUND 



LIN 

PUSH 

LD 

HL ; 

HL,(POINT) ; 

SAVE HL 

GET PTR 

LINIO 

LD 

INC 

CP 

JR 

U) 

POP 

RET 

A,(HL > } 

HL ; 

t r 

ZrLINiO ; 

(POINT)fHL f 
HL. ; 

GET NEXT CHAR 

ADVANCE PTR 

IGNORE BLANKS 

SAVE PTR 

RESTORE HL 

8 CP 

RET 

r 

STRING COMPARE 

POINT 

OEFS 

o » 

r 

TEMP PTR STORAGE 

ERROR 

DEFB 

o ; 

ERROR FLAGS FOR MATH PACKAGE 

ILLNUH 

EQU 

7 i 

ILLEGAL NUMBER (DTOB) 

OVRFLW 

EQU 

6 ; 

OVERFLOW 

UNDFLW 

EQU 

5 r 

UNDERFLOW 

DVDZRO 

EQU 

4 ; 

DIVISION BY ZERO 

SQRNGN 

EQU 

3 ; 

SQUARE ROOT OF NEGATIVE NUMBER 

LOGNGN 

EQU 

2 } 

LOG OF NEGATIVE NUMBER 


r 

i REFERENCES 


E R R 0 R , S C P , I... H ARK, I... R E S T , L B A C K , R C 2 , 


GLOBAL 



; CAL cm 

, ATOR 

PROGRAM US 

ING BCD MATH PACKAG 

t v r ? r J ? r 

V » » t 

7 T 7 7 

*********** 

77777777777 

**T**V*** , «**^*V**^* 

7777777777777777777 

IN IT 

LD 

(SSP) ,SP 

} SAVE STACK PTR 


CALL 

INOPN 

) OPEN TTY INPUT 


CALL 

OUT OPN 

OPEN TTY OUTPUT 


CALL 

RANDMZ 


INITIO 

LD 

HL, CI.RSCR 



CALL 

OUT ST R 

? CLEAR SCREEN 


LD 

HL,HEADNG 



CALL 

OUT ST R 

? PRINT HEADING 


CALL 

MENU 

? PRINT MENU 


CALL 

CL R ALL 

; CLEAR REGISTERS 


XOR 

A 



LD 

(ERROR),A 

? CLEAR ERROR FLAG 


CALL 

DP Y ALL 

} DISPLAY REGISTER 


CALL 

DPYMEM 

; DISPLAY MEMORY 

r 

; MAIN !. 

OOP 



? 

LOOP 

LD 

HL , CL RUN 



CALL 

OUT ST R 

? CLEAR INPUT LINE 


LD 

A, ( ERROR) 



AND 

A 

; CHECK FOR ERRORS 


LD 

HL , ERRMSG 



CALL 

NZ, OLITSTR 

; PRINT ERROR MESS 


CALL 

INCHR 

5 GET A CHAR 


LD 

HL,ERROR 



LD 

r p 

(HL) ,0 

t r 

; CLEAR ERRORS 


JR 

Z ,L00P05 

; ENTER NUMBER 


CP 

,Q, 



JR 

C,LOOP10 

) NOT A DIGIT 


CP 

7 9 r + 1 



JR 

NC.LOOP10 

; NOT A DIGIT 

LOOP05 

CALL 

NUM 

? A DIGIT 

LOOP 10 

CP 

' K ' 



JP 

Z>INITIO 

? RESTART 


LD 

HL,LOOZR 

; SET UP FOR ZERO 


LD 

IX,REGT 



PUSH 

HL 

; SET UP RET ADDR 


CP 

JP Z , RND 
CP 7 Z r 
JP Z,ZERO 
CP ’P’ 

JP Z,PI 
CP 

JP Z ,DUP ; DUPLICATE 

CP ' R' 

JP Z,MEHR 

LD IX,REG4 ; SET UP FOR UNARY OP 

LD HL,LOOUN 

EX (SP),HL ; SET NEW RET ADDR 

CP ' Q ' 

JP Z,SQR ? SQUARE ROOT 

CP ' E' 


4 


« 

i 

i 

< 

i 

< 

i 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

i 

i 

4 



i 


€ 

m 

t 

« 


t 

t 

« 

« 

t 

t 

« 


m 

m 

m 


LOOZR 
LOGBIN 


LOOUN 


CP 

' L' 



Jp 

Z , LOG 



CP 

r y r 



JP 

Z » TAN 



CP 

' C' 



Jp 

Z, COS 



CP 

'S' 



JP 

Z,SIN 



CP 

'A' 



JP 

Zi ATN 



CP 

'V' 



JP 

Zt REC 



CP 

t ft 



Jp 

ZfBRK 



CP 

' I' 



Jp 

Z, I NT 



CP 

7 B' 



Jp 

Zr ABS 



CP 

' N ’ 



Jp 

Z, NEC 



CP 

' X r 



JP 

Z»XCH 



CP 

'W r 



JP 

Z» MEHW 



ID 

IX,REG3 

* 

T 

SET UP FOR BINARY OPS 

L.D 

IY,REG4 



ID 

HL., L00B1N 



EX 

(SP),HL 

T 

NEW RET ADDR 

CP 




JP 

2, ADD 

r 

ADDITION 

CP 

r _ r 



JP 

2, SUB 



CP 




JP 

Z,MULT 



CP 

' l 7 



JP 

ZrDIVI 



CP 

r ,f. , 



JP 

Z.POW 



CP 

' D' 



RET 

*7 

t 

POP COMMAND 

LD 

HL,LOOP 

r 

SET UP FOR OTHERS 

EX 

(SP) ,HL. 

T 

NEW RET ADDR 

CP 

aDH 

t 

CR 

RET 

z 


CR IS ENTER 

LD 

A , OFFH 



I..D 

(ERROR) , A 

? 

FLAG COMMAND ERROR 

RET 


t 

T 

ALL UNDEFINED CHARS 

CALL 

PUSH 



JP 

LOOP 



L.D 

HL , REG4-1 

t 

t 

SRC 

LD 

DE,REGT"1 



LD 

BC , NUMSIZ*3 



L DDR 


r 

COPY REG1,2,3 INTO 2,3,4 

CALL 

DRYALL 

V 

r 

REDISPLAY REGS 

JP 

LOOP 



LD 

HL , DREG4 



LD 

IX , REG4 



CALL 

DPYR 

r 

REDISPLAY REG4 

CALL 

DPYROU 

r 

DISPLAY ROUNDED RESULT 




i 

< 

I 

« 

1 

I 

i 

1 

< 

« 

1 

< 

< 

< 

I 

« 

I 

I 

I 

I 


i 


J P l„ 0 0 P 

*H SPECIAL OPERATIONS 


; INPUT A NUMBER AND F'lJ 

SH 

IT 

? DIGITS , 

. , Er N, ARE 

ACCEPTED 

D AB.ORTS 

ENTER MODE 



I A DIGIT 

t 

HAS ALREADY 

BEEN TYPED 

t 

r 

} INITIALIZE 

BUFFER, FLAG 

, COUNT, PTR 

7 

NUM I..D 

(INBUF+1), 

A ? 

• 

SAVE FIRST DIGIT 

LD 

A , ’ + ' 



LD 

(INBUF),A 

7 

INIT SIGN 

LD 

HL,INBUF+2 



LD 

(INBPTR),HLJ 

INIT PTR 

LD 

B,INBLEN-2 



NUMIO LD 

<HL>,' ' 

T 

INITIALIZE BUFFER 

INC 

HL. 



D JNZ 

NUMIO 



LD 

A , 0 



L. D 

( NFI..G ) , A 

7 

CLEAR EXP FLAG 

LD 

A, ENTL.EN+2 



LD 

(ENTERL),A 

7 

INIT COUNT 

; CONVERT AND DISPLAY INBUF--REGT 

» SEND ERROR 

MESSAGES 



7 

NUM20 LD 

HL,ENTERL 



CALL 

OUT ST R 

7 

DISPLAY ENTER LINE 

CALL 

INCHR 

f 

GET A CHAR 

L D 

B, A 

7 

SAVE CHAR 

LD 

A,(NFLG) 



AND 

A 

7 

TEST FLAG 

LD 

A ,B 

7 

RESTORE CHAR 

JR 

NZ,NUM35 

V 

7 

IN EXPONENT PART 

r 

; PROCESS MANTISSA 



CP 

r O ? 



JR 

C, NUM25 

T 

NOT A DIGIT 

CP 

' 9 ' +1 



JR 

C, NUM60 

7 

A DIGIT 

NUM25 CP 

T T 

m 



JR 

2, NLIM60 

r 

A DECIMAL PT 

CP 

r E r 



JR 

Z , NUM65 

y 

SET UP FOR EXPONENT 

CP 

'N' 



JR 

t 

Z, NUM55 

y 

NEGATE 

7 

J END OF NUMBER INPUT, 

PU: 

SH RESULT 

t 

NUM30 CP 

T D ’ 



JR 

NZ,NUM32 

7 

NOT DELETE 

POP 

HL 

7 

KILL RET ADDR 

JP 

LOOP 

7 

IGNORE NUMBER 

NIJM32 PUSH 

AF 

7 

SAVE THIS CHAR 

LD 

HL ,INBUF 



LD 

IX,REGT 



CALL 

DT OB 

7 

CONVERT NUMBER IN REGT 



CALL PUSH i PUSH REGT ONTO STACK 

ID HL.CLRLIN 

CALL OUTSTR i CLEAR NUMBER DISPLAY 

POP AF 

RET 

¥ 

5 PROCESS EXPONENT PART 


NUM35 

CP 

'O’ 




JR 

Cr NUM37 

? 

NOT A DIGIT 


CP 

' 9 ' + 1 




JR 

C , NUM60 

r 

A DIGIT 

NUM37 

CP 

'N' 



* 

JR 

NZ»NUM30 

T 

NOT NEGATE 

r 

; NEGATE 

; EXPONENT 



r 

NUM40 

L.D 

HL,(ESGPTR) 

* 

T 

GET PTR TO SIGN 


LD 

A, (HL) 

r 

GET SIGN 


CP 

’ + ’ 




JR 

Z,NUM50 

r 

POSITIVE 


LD 

(HL),' + ' 




JP 

NUM20 



NIJM50 

LD 

(HL 




JP 

NUM20 



? 

; NEGATE 
* 

MANTISSA 



t 

NUM55 

LD 

r*' d 

A j (INBIIF) 

r « ? 

T 

GET SIGN 


L- r 

JR 

r 

Z t NLIM57 

r 

POSITIVE 


LD 

A; ' + ' 




LD 

(INBIIF) r A 




JP 

NUM20 



NUM57 

L.D 

A,’-’ 




LD 

( INBIJF) , A 




JP 

NUM20 



! STORE 

CHAR 

IN INBIIF 



f 

NUM60 

CALL 

NUM70 

r 

STORE CHAR 


JP 

NUM20 



r 

J SET EFLAG, 

ADJUST COUNT 

, PTR 

N LIMA 5 

CALL. 

NLIM70 

7 

STORE E IN INBUF 


LD 

A r OFFH 




LD 

( NFL.G ) , A 

» 

r 

TURN ON FLAG 


LD 

(ESGPTR),HL 

¥ 

SAVE PTR TO SIGN 


LD 

A; '+’ 


■ 


CALL 

NUM7Q 

7 

DEFAULT EXP SIGN 


JP 

NUM20 




; ROUTINE TO STORE INTO INBIIF 
? RETURNS PTR VALUE IN HL 


NUM70 

LD 

B. A 

; SAVE CHAR 


LD 

HL r(INBPTR) 

r GET PTR 


LD 

A,(ENTERL) 

; GET COUNT 


CP 

ENTLEN+INBL 

EN 


IGNORE 


RET 2 ; COUNT MAX, 

INC A i INC IT 

ID (ENTERL ) , A 
L.D (HL) , B ; STORE CHAR 

INC HL 

LD (IMBF'TR) , HL 

RET 


? DUP DUPLICATE TOP OF 

t 

r 

DUP LD HL,REG4 

LD DE,REG T 
L.D BC, NIJMSIZ 
L DIR 
RET 

T 

EXCHANGE REGISTERS 


XCH 

LD 

DE.REGT+7 


LD 

HL,REGT-1 


LD 

L DDR 

BC,NUMSIZ 


LD 

HL,REGT 


LD 

DE,REG3 


LD 

L DIR 

8C,NUMSIZ 


LD 

HL,DREG3 


LD 

IX,REG3 


CALL 

DPYR 


RET 


} EXIT 
* 

TO RIO 

t 

RIO 

LD 

SP,(SSP) 


RET 


7 

! PI ROUTINE 

RETURNS P 

PI 

PUSH 

IX 


POP 

DE 


LD 

HL , CPI 


LD 

LDIR 

RET 

8C,NUMSIZ 


T 

? MEMORY READ 


1 

ME MR 

CALL 

GETDIG 


JR 

C,MEMR99 


PUSH 

IX 


POP 

HL 


LD 

DE , REGT 


LD 

BC t NUMSI 


LDIR 



RET 


MEMR99 

POP 

HL. 


JP 

LOOP 


TACK 

COPY REG4 TO REGT 

DF.ST 

SRC 

COPY REG3,4 TO REG4,T 

COPY REGT TO REG3 

DISPLAY NEW REG3 

RESTORE SP 
RETURN TO RIO 

COPY IX TO DE 


GET MEM ADDR 
ERROR 

COPY REG PTR TO HL. 
PT TO TEMP REG 
COPY VALUE 
THROW AWAY RET ADDR 


J MEMORY WRITE 



ME MW 

CALL 

GETDIG 

; GET MEM ADDR 



P \Z J 

0 

j ERROR 



PUSH 

HL. 

; SAVE DISPLAY PTR 



PUSH 

IX 

? COPY REG PTR TO DE 



POP 

DE 




L.D 

HL , REG4 

! SRC IS REG4 



LD 

BC,NUMSIZ 




L DIR 


; COPY REG4 TO MEM 



POP 

HL 

GET DPY PTR 



L A L L 

DF'YR 

; REDISPLAY MEM REG 



RET 



* H 

SPECIAL FORMATTING 


DPYROU 

LD 

HL , REG4 

? SRC 



LD 

DE,REG T 




LD 

BC,NUMSIZ 




LDIR 


; COPY REG4 TO REGT 



LD 

IX, REGT 




LD 

IY,ROUCON 

; ROUNDING CONSTANT 



LD 

A, ( IX + NIJM8I 

Z-l ) 



LD 

(IY+NUMSIZ- 

1 ) , A 



LD 

A, (IX) 

; GET SIGN 



AND 

80H 




LD 

( IY) , A 




CALL 

ADD 

J ROUND TO 10 PLACES 

y 

; c 

:LEAR 

DISPLAY AREA 


r 


LD 

HL,DROU+ADR 

SI 2+1 



LD 

B, 18 


DFYRO5 

LD 

(HL),' ’ 




INC 

HL 




D JNZ 

DPYR05 




LD 

HL,DREG T 




PUSH 

HL 

? SAVE HL. 



CALL 

BTOD 

t CONVERT NUMBER 



EX 

AF.AF* 

J SAVE EXPONENT IN A' 

y 

? b 

ET UP HL 

F'TING TO INPUT, DE TO OUTPUT STRINGS 

y 


POP 

HL 

; GET PTR TO ASCII STRING 



LD 

DE,DROU+ADR 

SIZ+1J OUTPUT PTR 

? c 

HECK 

FOR 

ZERO INPUT 


y 


XOR 

A 

? CHECK FOR ZERO 



CP 

( IX+7) 


* 


JR 

Z.STR90 

; INPUT IS ZERO 

T 

! CALCUL 

ATE 

# SIGNIFICANT DIGITS 

y 


PUSH 

HL 

! SAVE INPUT PTR 



LD 

BC , 10 

l MAX COUNT 



ADD 

HL , BC 

; PT TO LAST DIGIT 



LD 

A, '0' 

; ASCII ZERO 

SIR 

10 

CP 

(HL) 

5 COMPARE DIGIT 



JR 

NZ.STR20 

1 DIGIT IS SIGNIFICANT 



DEC 

HL 

DECR PTR 



DEC 

c 

; DECR COUNTER 



JR 

NZ,STR10 

J COUNTER NOT ZERO 



POP 

HL 

,’ RESTORE PTR 





€ 

€ 

« 

« 

i 

I 

i 

t 

i 

i 

t 

l 

i 

i 

i 

€ 

t 

« 

t 

i 

€ 


JR 

STR 90 

7 

THIS JR SHOULDN’T HAPPEN 

3TR20 POP 

HL 

t 

RESTORE PTR 

J SET SIGN 




r 

LD 

A,(HL) 



LD 

(DE) f A 



INC 

HL 



INC 

DE 



7 

; CHECK IF X 

>= 1 



? 

EX 

AF,AF’ 

? 

GET EXPONENT 

LD 

B, A 

T 

t 

SAVE IN B 

AND 

A 



JR 

Z,STR40 

? 

EXPONENT IS ZERO 

JP 

M,STR60 

t 

EXPONENT IS < 0 

CP 

11 



JP 

NC , STR80 

7 

TOO BIG 

t 

) PRINT DIGIT 

‘S TO LEFT 

OF 

DECIMAL PT 

r 

STR30 LD 

A, (HL. ) 

r 

GET CHAR 

LD 

(DE),A 

7 

STORE IN STRING 

I NC 

HL 

t 

ADVANCE PTRS 

INC 

DE 



DEC 

C 

7 

DECR # SIGNIF DIGITS 

DJNZ 

STR30 

7 

EXP COUNTER NOT ZERO 

7 

! CHECK IF ANY SIGNIFICANT DIGITS LEFT 

f 

L.D 

A , C 

7 

GET # DIGITS 

AND 

A 



JP 

M, STR100 

7 

NONE LEFT 

J R 

Z,STR100 

J 

NONE LEFT 


J PRINT DECIMAL PI 

7 

STRAO LD A,’.’ 

LD (DE),A 
INC DE 

7 

J. PRINT SIGNIFICANT DIGITS TO RIGHT OF PT 

7 

STR50 LDIR ? (B IS ZERO, C IS COUNTER) 

JR STR100 


; CHECK IF NUMBER WILL FIT IN 10.PLACE FRACTION 


3TR60 

LD 

A, C 

? GET # 

SIG DIGITS 


SUB 

CP 

B 

11 

r ADD # 

LEADING ZEROS 


JR 

NC,STR80 

i WON'T 

FIT, USE E NOTATION 


7 

; OUTPUT DECIMAL PT AND LEADING ZEROS 


LD 

A, V 

LD 

(DE) , A 

INC 

DE 

L D 

A, ’0’ 


i 

c 

€ 

f 

t 

€ 

€ 

f 

i 

i 

< 

i 

i 

< 

I 

i 

i 

4 

I 

4 

4 


I 


ASCII 0 




SIR 70 

LD (DE) , A 

7 

STORE A 0 


INC DE 

x 

ADV PTR 


INC B 

T 

INC COUNTER 


JR NZ j STR70 

t 

COUNTER NOT 0 


JR STR50 

r 

GO PRINT SIG DIGITS 

r 

; USE E 

NOTATION FOR NUMBER 

r 

STR80 

LD BO,10 

1 

•DIGITS 


LDI 

i 

TRANSFER LEADING DIGIT 


L D A f ’ . ’ 




LD (DE)>A 

y 

DECIMAL PT 


INC DE 

L DIR 

* 

REST OF DIGITS 


INC HL 

INC HL 


SKIP 3 DIGITS 


INC HL 

INC DE 

r 

SKIP A SPACE 


LD BC 1 5 

LDIR 

r 

CHARS LEFT IN EXP 


JR STR100 



7 

; OUTPUT A ZERO 



T 

STR90 

LD A , ! 0 r 

INC DE 




LD ( D E ) , A 



8TR100 

LD HLfDROU 

CALL OUTSTR 

RET 

V 

7 

DISPLAY IT 

* H UTIL 

ITIES 



r 

? DISPLAY A REGISTER 



} HI. 

PTS TO DISPLAY 

BUFFER 

* T V 

t X /\ 

PTS TO REGISTER 

VALUE 

? 

DPYR 

PUSH HL 

J 

SAVE OUTPUT PTR 


PUSH IX 

POP HL 

r 

COPY IX TO HL 


LD DEfREGT 

t 

DEST IS TEMP REG 


LD BC,NUMSIZ 

LDIR 

7 

COPY REGISTER TO TEMP 


LD IX,REGT 

t 

PT IX TO IT 


LD HLrDREGT 

7 

PT TO TEMP DRY REG 


CALL BTOD 

7 

CONVERT TO ASCII 


POP HL 

PUSH HL 

7 

GET PTR 


LD BC,ADRSIZ+3 



ADD HL , BC 

EX DE , HL 

7 

PT TO SIGN BYTE 


LD HLrDREGT 

7 

PT TO BTOD OUTPUT 


LD BC , NUMDIG+1 



LDI 

7 

COPY SIGN 


LDI 

L D A , ’ . ' 

7 

COPY FIRST DIGIT 


LD {DE),A 

7 

DECIMAL PT 


INC DE 

LDIR 

7 

COPY REST OF DIGITS 


LD A,' r 

LD (DE) , A 

7 

BLANK 



INC HL 


ID BC,NUMEX P 


LDIR 

; EXPONENT 

POP HL 


CALL OUTSTR 

; PRINT RESULT 

RET 



r 

; PUSH, DP Y ALL. 

: VALUE IN TEMP REG IS PUSHED ONTO STACK 

i ALL REGISTERS ARE REDISPLAYED 


PUSH 

L..D 

HL,REG2 J SOURCE 


LD 

DE,REGi J DEST 


L D 

BC,NUM8IZ*4 


LDIR 

? MOVE ALL VALUES 

DPYALL 

LD 

HL,DREG! 


LD 

IX,REG1 


CALL 

DPYR 


L D 

HL , DREG2 


LD 

IX,REGS 


CALL 

DPYR 


LD 

HL,DREG3 


L.D 

IX,REG3 


CALL 

DPYR 


LD 

HL , DREG4 


LD 

IX,REG4 


CALL 

DPYR 


CALL 

RET 

DPYROU ; DISPLAY ROUNDED 

; CL.. R ALL 

CLEAR ALL REGISTERS 

t 

CL RAIL. 

LD 

HL., REG1 


LD 

B,NUMSIZ*15 

CL RAID 

LD 

(HL),0 


INC 

Hl. 


DJNZ 

RET 

CLRAiO 


i DISPLAY MEMORIES 


7 

DPYMEM 

LD 

HL,DMEMO 

7 

FIRST MEM REG 


LD 

IX,MEMO 




!.. D 

B , 10 

7 

COUNT 

DPYM10 

PUSH 

BC 

r 

SAVE COUNT 


PUSH 

HL 

7 

SAVE REG PTE 


PUSH 

IX 

r 

SAVE DPY PTR 


CALL 

DPYR 




POP 

IX 




L D 

BC,HUMSIZ 

7 

DISPLAY SIZE 


ADD 

IX , BC 

t 

PT TO NEXT DIS 


POP 

HL 




LD 

BCfOSIZE+1 




ADD 

HL , BC 

7 

PT TO NEXT REG 


POP 

BC 




DJNZ 

DPYM10 

7 

COUNT NOT 0 

? 

RET 





; DISPLAY MENU 



MENU 


LD HL » MEN 111 

CALL OUTSTR 
LD HL,MENU2 
CALL OUTSTR 
LD HL , MENU3 
CALL OUTSTR 
LD HL., MENU4 

CALL OUTSTR 
LD HL , MENUS 
CALL OUTSTR 
LD HL ,MENU6 
CALL OUTSTR 
RET 


J GETDIG GETS ADDRESSES FOR MEMORY OPERATIONS 
J INPUTS A CHAR, MULTIPLIES BY NIJMSIZ AND DSIZE, 
! ADDS TO REG AND DPY F'TRS 
; RETURNS CY IF NOT SUCCESSFUL 


GETDIG 

CALL 

INC HR 


CP 

' 0 7 


JR 

C.GETD99 


CP 

' 9' +1 


JR 

NC»GETD99 


AND 

OFH 


LD 

C, A 

GETDIO 

XOR 

A 


LD 

B, NUMSIZ 

GETD20 

ADD 

A , C 


DJNZ 

GETD20 


LD 

IX,MEMO 


L.D 

E, A 


LD 

D , 0 


ADD 

IX, DE 


EX 

DE»HL 


LD 

E, C 


LD 

D , 0 


LD 

B,DSIZE+i 

GETD30 

ADD 

HL, DE 


DJNZ 

GETR30 


LD 

DE,DMEMO 


ADD 

HL, DE 


AND 

RET 

A 

GETD99 

CP 

c* r* rr 

' D 7 


•:* t* r 

RET 

Z 


LD 

A,OFFH 


LD 

RET 

(ERROR),A 

*H I/O 

UTILITIES 

T 

r INPUT 

FILE 

OPEN 

r 

INOPN 

LD 

IY,INVECT 


LD 

(IY+1),4 


LD 

(IY + 4) ,0 


CALL 

RET 

SYSTEM 


; GET A CHAR 

; NOT A DIGIT 

5 NOT A DIGIT 
f CONVERT TO BINARY 
i SAVE IN C 
} CLEAR A 


I MULTIPLY BY NUMSIZ 

; PT TO FIRST REG 
} PUT RESULT IN DE 


J PT TO ADDRESSED MEM 
} PUT DISPL. IN HL 
? PUT ORIGINAL VAL IN DE 

NUMSIZ 

; MULTIPLY BY DSIZ 

; PT TO FIRST DPY 
} PT TO ADDRESSED MEM DPY 
; CY OFF 


J TURN ON CY FLAG 
1 DELETE COMMAND 

5 FLAG ERROR 


? PT TO VECTOR 
; OPEN 
; LEN 






OUTOPN 


LD IY,OTVECT ; F'T TO VECTOR 

LD (IY+1),4 ? OPEN 

I D (IY+4)» 0 ; LEN 

CALL SYSTEM 
RET 

t 

: INCHR INPUT 1 CHAR 


RETURNED 

IN A 




\R LD 

IY,INVECT 




LD 

(IY+1),QCH 

r 

READ 

LINE 

LD 

(IY+4)ii 

? 

DATA 

LENGTH 

LD 

HL , IMCBIJF 




LD 

(IY+2),L 

7 

DATA 

ADDR 

LD 

(IY + 3), H 




CALL 

SYSTEM 




L D 

(IY+1),4 0H 

7 

READ 

STATUS 

L.D 

{IY + 4) ,1 

1 

DATA 

LENGTH 

LD 

HL -CSTAT 


GET 

CSAT ADDR 

LD 

( IY + 2) r L. 




L D 

(IY + 3)> H 




CALL 

SYSTEM 




L.D 

A, (CSTAT) 

7 

GET 

FLAGS 

BIT 

5 r A 




JP 

2, RIO 

t 

7 

ESC 

PENDING 

LD 

A,(INCBUF) 

t 

GET 

CHAR 

BIT 

6 > A 




RET 

Z 

? 

NOT 

UPPER CAS 

RES 

5 , A 

t 

T 

MAKE 

UPPER CA 


RET 

t 

f 

? OUTSTR OUTPUT A STRING 
; HL PTS TO COUNT FOLLOWED BY STRING 


OUTSTR 


LD A>(HL > ? GET COUNT 

LD IY i OTVECT 

LD (IY+1 ) > 1 OH WRITE ASCII 

LD (IY + 4)»A 

INC HL 

L.D ( IY + 2) ,L 

LD (IY+3);H 

CALL SYSTEM 

RET 


; VECTORS 


INVECT DEFB 1 
DEFB 4 
DEFW 0 
DEFW 0 
DEFW 0 
DEFW 0 
DEFB 0 
DEFW 0 
OTVECT DEFB 2 
DEFB 4 
DEFW 0 
DEFW 0 
DEFW 0 
DEFW 0 


DEFB 0 
DEFW 0 

*H CONSTANTS, WORK AREAS 


; EQUATES 


r 

ASCESC 

EQU 

1BH 

SYSTEM 

EQU 

14 0 3 H 

NUMSIZ 

EQU 

8 ', 

NIJMDIG 

EQU 

13 ; 

NIJMEXP 

EQU 

4 ; 

r 

; cons 

TANTS 


T 

CPI 

DEFW 

1403H 


DEFW 

9215H 


DEFW 

3565H 


DEFW 

819 0 H 

ROUCON 

DEFW 

0 


DEFW 

0 


DEFW 

50 OH 


DEFW 

0 

T 

; MESS 

AGES, £ 

STRINGS 

r 

HEADNG 

DEFB 

27+ADRSIZ 


AD DR 

1,20 


DEFM 

[7. 8 0 C A 

? 

; MENU 

LINES 


T 

MENU! 

DEFB 

60+ADRSIZ 


AD DR 

4,8 


DEFM 

’+ Y + X 


DEFM 

'Q SQR(X) 


DEFM 

'T TAN(X) 


DEFM 

'V 1/X 


DEFM 

'X EXCH X,Y 

MENU 2 

DEFB 

60+ADRSIZ 


AD DR 

5,8 


DEFM 

7 _ v — V 

1 A 


DEFM 

’E EXP(X) 


DEFM 

'S SIN(X) 


DEFM 

'N CHG SGN 


DEFM 

'= DUP X 

MENUS 

DEFB 

60+ADRSIZ 


A DDR 

6,8 


DEFM 

’* Y*X 


DEFM 

'L LOG(X) 


DEFM 

’C COS(X) 


DEFM 

' B ABS(X) 


DEFM 

’D DEL/POP 

MENIJ4 

DEFB 

60+ADRSIZ 


AD DR 

7,8 


DEFM 

'/ Y/X 


DEFM 

’# RANDOM 


DEFM 

7 A ARCTAN(X) 


DEFM 

r F FRACTION 


DEFM 

' R < N > READ 

MENUS 

DEFB 

60+ADRSIZ 


ESCAPE CHAR 

SIZE OF DATA 

NUMBER OF SIGNIF. DIGIT 

NUMBER OF EXPON DIGITS 


3.141592653590 


L C U L A T 0 R’ 


r 

? 


r 

r 


f 

r 


r 

r 

r 







AD DR 

8 , 8 



DEFM 

'+ Y + X 

y 


DEFM 

'Z ZERO 

7 


DEFM 

'P PI 

T 


DEFM 

'I INTEGER 

T 


DEFM 

’ W < N> WRITE 

t 

MENIJA 

DEFB 

60+ADRSIZ 



ADDR 

9,8 



DEFM 

7 

7 


DEFM 

7 

7 


DEFM 

r 

7 


DEFM 

’ <ESC> EXIT 

7 


DEFM 

'K CLEAR ALL 

7 

ERRMSG 

DEFB 

5+ADRSIZ 



ADDR 

C i~~ r X 



DEFM 

'ERROR' 


NOP 

DEFB 

1 



DEFB 

r r 


CLRLIN 

DEFB 

39+2*ADRSIZ 



ADDR 

n O i * 

JL *1 f X r 

ROW 


DEFM 

7 



DEFM 

1 



ADDR 

22,1 i 

ROW 

PRHPTA 

DEFB 

2+ADRSIZ 



ADDR 

c. a. r X r 

ROW 


DEFB 

7 > r 


FRSTDG 

DEFS 

1 


PROMPT 

DEFB 

l 



DEFB 

T 7 


CLRSCR 

DEFB 

CLRSIZ 



CLEAR 



7 

! DISPLAY BUFFERS 


7 

DSIZE 

ECU.) 

22+ADRSIZ 

DREG! 

DEFB 

DSIZE 


ADDR 

13,8 


DEFM 

7 7 


DEFS 

20 

DREG 2 

DEFB 

DSIZE 


ADDR 

14,8 


DEFM 

7 7 


DEFS 

20 

DR EG 3 

DEFB 

DSIZE 


ADDR 

15,8 


DEFM 

' Y ' 


DEFS 

20 

DR EG 4 

DEFB 

DSIZE 


ADDR 

16 , 8 


DEFM 

'X ' 


DEFS 

20 

DREGT 

DEFS 

NUMDIG+NUMEXP+ 2 

DMEMO 

DEFB 

DSIZE 


ADDR 

13,40 


DEFM 

'0 ’ 


DEFS 

20 

DMEM1 

DEFB 

DSIZE 


ADDR 

14,40 


DEFM 

' 1 ’ 


DEFS 

20 





DMEM2 

DEFB 

DSIZE 


ADDR 

15,4Q 


DEFM 

? r 


DEFS 

20 

DMEM3 

DEFB 

DSIZE 


ADDR 

16,40 


DEFM 

'3 ' 


DEFS 

20 

DMEM4 

DEFB 

DSIZE 


ADDR 

17,40 


DEFM 

'4 ' 


DEFS 

20 

DMEM5 

DEFB 

DSIZE 


ADDR 

18,40 


DEFM 

, 5 , 


DEFS 

20 

DMEM 6 

DEFB 

DSIZE 


ADDR 

19,4 0 


DEFM 

'6 ’ 


DEFS 

20 

DMEM7 

DEFB 

DSIZE 


ADDR 

20,4 0 


DEFM 

'7 ’ 


DEFS 

20 

DMEM 8 

DEFB 

DSIZE 


ADDR 

21,40 


DEFM 

’8 ’ 


DEFS 

20 

DMEM9 

DEFB 

DSIZE 


ADDR 

22,40 


DEFM 

> 9 » 


DEFS 

20 

DROU 

DEFB 

18+ADRSI 


ADDR 

19,10 


DEFS 

18 

y 

5 REGIS 

TERS 


r 

REGl 

DEFS 

NUMSIZ 

REG 2 

DEFS 

NIJMSIZ 

REG3 

DEFS 

NUMSIZ 

REG4 

DEFS 

NUMSIZ 

REGT 

DEFS 

NUMSIZ 

MEMO 

DEFS 

NUMSIZ 

MEM 1 

DEFS 

NUMSIZ 

MEM2 

DEFS 

NUMSIZ 

ME M3 

DEFS 

NUMSIZ 

MEM4 

DEFS 

NUMSIZ 

MEM5 

DEFS 

NUMSIZ 

MEM 6 

DEFS 

NUMSIZ 

MEM7 

DEFS 

NUMSIZ 

MEMS 

DEFS 

NUMSIZ 

MEM9 

DEFS 

NUMSIZ 

T 

? BUFFERS 


f 

ENTLEN 

EQU 

7+ADRSIZ 

INBLEN 

EQU 

30 

ENTER!... 

DEFS 

1 


ADDR 

22,1 


i SIZE OF ENTER MSG 
J SIZE OF INBUF 
J MSG LENGTH 
} ENTER MESSAGE 



INBUF 

DEEM ’ENTER; 
DEFS INBl.EN 

NFLG 

DEFS 1 

DEFS 1 

INBPTR 

DEFS 2 

ESGPTR 

DEFS 2 

I NCBUF 

DEFS 1 

CSTAT 

DEFS 1 

SSP 

DEFS 2 

! REFERENCES 


EXP SIGN HAY GO HERE 
NUMBER FLAG 
PTR IN INBUF 
PTR TO EXP SIGN 
INPUT CHAR 
CONSOLE STATUS 
SAVED STACK PTR 


EXTERNAL 

EXTERNAL 

EXTERNAL 

EXTERNAL 

EXTERNAL 

EXTERNAL 


DTOB,BTOD 

ADD,SUB,MULT,DIVI,ROW 

INT » ABS » SGN > BRK » RND 

REC, NEG ,EANDMZ,ZERO 

SQR,EXP,LOG r TAN jCOS,SIN,ATN 

ERROR 


ASM #.1. CAL (M NOL. 0 = CAL. OBJ) 

LINK $-4400 CAL. DMATHM DMATHO OMATHJ. DMATH2 (NOM ) 


*H MULTIPLICATION 


8 BIT BY 8 BIT UNSIGNED MULTIPLY 306 CYCLES AVG. 

MULTIPLICAND IN E 
MULTIPLIER IN C 
PRODUCT IN AE 


MP88 XOR 

A 

l.D 

B t 8 

RR 

E 

MP88A JR 

NC r $+3 

ADD 

A, C 

RRA 


RR 

E 

DJNZ. 

MP88A 

RET 


*EJECT 


T 

! 16 BIT BY 

16 BIT IJNS 

T 

! MULTIPLICAND IN DE 

; MULTIPLIER IN AC 

; PRODUCT 

IN HL.AC 

t 

M P16 LD 

HLrO 

LD 

8,16 

RRA 


RR 

c 

MF'161 JR 

NC.MP162 

ADD 

HI., DE 

MP162 RR 

H 

RR 

L 

RRA 

r« r« 


K k 

DJNZ 

MP161 

RET 


*EJECT 


t 

r 24 BIT BY 

24 BIT UNS 

t 

MULTIPL ICAND IN CD! 

; MULTIPLIER IN H' L' 

; PRODUCT 

IN AHLH'L’I 

T 

MP24 XOR 

A 

LD 

Hr A 

LD 

L, A 

L D 

B»24 

EXX 


RR 

H 

RR 

L 

RR 

D 

EXX 


MP241 JR 

NC,MP242 

ADD 

HI.»DE 

ADC 

A r C 

MP242 RRA 


RR 

H 


CLEAR ACCUM 

SET COUNTER 

1ST MPLIER BIT IN CY 

BIT WAS A 0 

ADD MULTIPLICAND 

ROTATE ACC-MULT I PL. IER INTO CY 

COUNTER NOT 0 


CLEAR ACCUM 
SET COUNTER 
1ST MPLIER BIT IN CY 

MPLIER BIT IS 0 

ADD MULTIPLICAND 

ROT ACC-MPLIER INTO CY 


COUNTER NOT 0 


2022 CYCLES AVG, 


CLEAR ACCUM 


SET COUNTER 
ALT BANK 

ROT 1ST MPLIER BIT INTO CY 


MAIN BANK 
MPLIER BIT WAS 0 
ADD MPLICAND 

ROT ACC-MPLIER INTO CY 




EKK 

P P 

LI 

! ALT BANK 

K ft. 

RR 

a p 

n 

L 

n 


l\ ft. 
p V 

u 

5 MAIN BANK 

OJNZ 

RET 

MPZ41 

l COUNTER NOT 0 


*EJf=CT 

j 

; 32 BIT BY 32 BIT UNSIGNED MULTIPLY 36.1.9 CYCLES AVG. 

f 

i MULTIPLICAND IN D’E’DE 
J MULTIPLIER IN B'C'CA 


; PRODUCT 

IN H'L'HLB' 

C ’ C A 

? 

HP32 I..D 

HL , D 

5 CLEAR ACCUM 

EKK 


; ALT BANK 

L D 

H L , 0 


EKK 


; MAIN BANK 

LD 

B, 33 

1 SET COUNTER 

JP 

MP323 


HP321 JR 

NC,MP322 

; MPLIER BIT IS 0 

ADD 

HL, OE 

} ADD MULTIPLICAND 

EKK 


! ALT BANK 

ADC 

HL , DE 


EKK 


J MAIN BANK 

MP322 EKK 


? ALT BANK 

RR 

H 

; ROTATE ACCUM INTO CY 

RR 

L 


EKK 


; MAIN BANK 

RR 

D P 

H 

I 


ft ft. 

HP323 EKK 

L 

; ALT BANK 

RR 

B 

? ROTATE MPLIER INTO CY 

RR 

C 


EKK 


! MAIN BANK 

RR 

c 


RRA 



D JN2 

MP321 

; COUNTER NOT 0 

RET 



*H DIVISION 




’ 8 BIT BY 8 BIT UNSIGNED DIVIDE 346 CYCLES AVG. 
i USES RESTORING DIVISION, PRODUCING COMPLEMENTED QUOTIENT 


i 

E 

IS DIVIDEND 

r 

c 

IS DIVISOR 

? 

A 

IS QUOTIENT 

T 

B 

IS REMAINDE 


DV88 

KOR 

A 

V 

CLEAR ACCUM 


LD 

B, 8 

} 

LOOP COUNTER 

D V B A 

RL 

E 

1 

ROTATE CY INTO ACC-DIVIDEND 


RLA 


T 

CY WILL BE OFF 


SUB 

C 

r 

TRIAL SUB DIVISOR 


JR 

NC, 

7 

SUBTRACT OK 


ADD 

A , C 

\ 

RESTORE ACCUM, SET CY 


DJNZ 

DV8A 




L.D 

B, A 


PUT REMAINDER IN B 


LD 

A, E 

; 

GET QUOTIENT 


RLA 


* 

SHIFT IN LAST RESULT BIT 



CPL i COMPLEMENT BITS 

RET 

♦ EJECT 


I 16 BIT BY 8 BIT UNSIGNED DIVIDE 473 CYCLES AVG. 
; RESTORING DIVISION, FAKING A 9 BIT ACCUM 

f 

I DIVIDEND IN AE 

J DIVISOR IN C 

J QUOTIENT IN E 

? REMAINDER IN A 

; . CY ON FOR OVERFLOW (A >= C ON INPUT) 


D V16 8 

CP 

CCF 

c 

J 

CHECK FOR OVERFLOW 


RET 

c 

r 

ft >~ c 


LD 

B , 8 

• 

LOOP COUNTER 


JP 

DV 163 



DV161 

SUB 

c 

T 

SUB DIVISOR 

D V16 2 

SCF 


? 

FORCE RESULT TO 1 


DEC 

B 




JR 

Z,DV165 

} 

COUNTER IS 0 

D V16 3 

RL 

E 

r 

ROTATE CY INTO ACC-DIVIDEND 


RLA 





JR 

C , DV161 

7 

NEED 9 BIT ACCUM 


SUB 

c 

? 

TRY SUB DIVISOR 


JR 

NC,DV162 

? 

SUB OK, OUTPUT A 1 


ADD 

A , C 

r 

RESTORE ACCUM 


AND 

A 

r 

MAKE CY 0 


DJNZ 

DV163 

r 

COUNTER NOT 0 

D V16 5 

RL 

E 

r 

GET LAST RESULT BIT, 0 CY 


RET 




♦EJECT 






? 

i 16 BIT BY 15 BIT UNSIGNED DIVIDE 1164 CYCLES AVG. 
f NON--RESTORING DIVISION ALGORITHM 


DIVIDEND IN AC 
DIVISOR IN DE (15 BITS) 
QUOTIENT IN AC 
REMAINDER IN HL 


D V15 
D V151 


D Vi 52 
D V15 3 


D V15 4 


LD HL , 0 
LD B,16 

RL C 
RLA 

ADC HL , HL. 

SBC HL,DE 
CCF 

JR NC,DV155 

DJNZ DV151 

RL C 

RLA 

RET 

RL C 
RLA 

ADC Hl., HL 
AND A 


; CLEAR ACCUM 
} SET COUNTER 
; 0 <~ HL < DE 
? ROTATE ACC-RESULT LEFT 

J NO CY POSSIBLE 
i SUBTRACT DIVISOR 
} CALC RESULT BIT 
r ACCUM WENT NEGATIVE 
i COUNTER NOT 0 
! SHIFT IN LAST RESULT BIT 


; -DE <= HL < 0 
r ROTATE ACC-RE8ULT LEFT 


t 


TURN OFF CY 








